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FOREWORD 

This  technical  report  is  submitted  to  the  Georgia  Institute  of 
Technology  to  comply  with  the  report  requirements  of  contract  l-A-2550, 
which  is  a  subcontract  under  United  States  Navy  contract  N-00039-80-C- 
0032.  This  report  is  published  in  two  volumes,  and  each  volume  consists 
of  two  parts. 
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REPORT  SUMMARY 


Auburn  University,  under  contracts  N66314-73-C-0565,  N66314-74-C-1352, 
N66314-74-C-1634,  N00228-75-C-2080,  N00228-76-C-2069,  and  N00228-78-C-2233 
with  the  United  States  Navy,  and  has  investigated  various ^spects  of  the 
/Marine  Air  Traffic  Control  and  Landing  System  (MATCALSj.  This  report 
contains  the  results  of  the  continuation  of  ^hese^nvestigations^under 
contract  1-A-255G  with  the  Georgia  Institute  of  Technology  .jf The  report 

is  organized  into  three  main  sections,  namely  Part  Two,  Part  Three,  and 

/ 

Part  Four.  Part  Two  presents  a  method  of  estimating  the  centroid 
location  of  a  target  utilizing  a  scan  return  amplitude  versus  angle 
information.  Part  Three  contains  the  results  of  an  investigation  into 
replacing  the/a-j£  filter  in  the  MATCAL  digital  controller  with  an  observer, 
in  order  to  reduce  the  effects  of  radar  noise.  Part  Four  presents  the 
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results  of  an  investigation  into  replacing  the  same  a-e  filter  with  a 


tri -state  adaptive  filter,  in  order  to  reduce  the  effects  of  radar  noise. 
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Centroid  Estimation 


A 
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Essential  to  the  performance  of  any  tracking  radar  is  an  effective 
target  centroid  estimator.  The  purpose  of  the  work  reported  in  Part  Two 
is  to  examine  the  accuracy  of  several  target  centroid  estimators  in  a 
comparative  fashion,  and  to  introduce  a  non-thresholding  algorithm 
developed  as  part  of  this  research.  This  analysis  was  conducted  using  a 
software  simulation  of  a  landing  system  radar  tracking  a  passive  target. 
The  algorithm  developed  in  Part  Two  is  a  method  of  estimating  the  centroid 
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location  of  a  target  utilizing  scan  return  amplitude  versus  angle  infor¬ 
mation.  The  method  is  compared  to  three  thresholding  estimators  and  a 
first  moment  estimator  in  a  computer-simulated  automatic  landing  system. 

It  was  found  that  the  method  introduced  was  the  most  robust  and  ac¬ 
curate  of  the  estimators  in  noise,  due  to  its  unique  scan  rejection  capa¬ 
bility.  In  periods  of  high  signal-to-noise  ratio  the  method  had  less 
error  than  the  thresholding  methods,  and  was  similar  in  ability  to  the 
first  moment  estimator.  Further,  the  pulse  transmissions  required  to 
obtain  a  desired  level  of  performance  is  much  reduced  from  the  thresholding 
methods  employed  in  this  simulation. 

Observer  Design 

Presently  a  problem  exists  in  the  closed- loop  control  of  the  MATCALS 
system  due  to  the  noise  generated  in  AN/TPN  22  radar.  An  a-B  filter  in 
the  flight  dynamic  and  control  module  is  employed  to  reduce  the  noise 
effects  while  estimating  the  position  and  the  velocity  of  the  aircraft. 

An  observer  may  also  be  used  to  estimate  the  status  of  the  aircraft.  Part 
Three  of  this  report  presents  the  results  of  an  investigation  of  the 
replacement  of  the  a-B  filter  with  an  observer. 

The  F4J  aircraft  lateral  control  system  is  employed  as  an  example 
in  this  investigation.  Several  different  controllers  are  utilized  to 
determine  which  yield  the  best  radar-noise  response  and  which  yield  the 
best  wind  response. 

The  proposed  MATCALS  system  contain  an  a-B  filter  in  the  controller. 


Alternative  controllers  are  constructed  by  replacing  the  a-B  filter  with  an 
observer. 
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In  geniral  the  observer  control  systems  exhibit  significantly  less 
radar-noise  response  than  do  the  a-s  systems,  but  exhibit  somewhat  more 
wind  response.  The^e  studies  indicate  that  the  observer  controllers  im¬ 
prove  the  MATCALS  system's  operation  when  compared  to  the  a-e  controllers, 
and  that  the  observer  systems  should  be  considered  further. 

Tri-State  Adpative  Filters 

A  tri-state  adaptive  tracking  filter  was  designed  for  use  in  the 
F4J  aircraft  lateral  control  system  in  an  automatic  landing  configuration.  1 

The  system  presently  uses  an  a-6  tracking  filter  to  estimate  the  aircraft's 
lateral  position  and  velocity.  The  tri-state  adaptive  filter  is  designed 
to  replace  the  a-e  filter. 

Three  digital  tracking  filters,  each  based  upon  a  different  air¬ 
craft  dynamic  model,  were  combined  to  form  the  tri-state  adaptive  tracking 
filter.  The  selection  of  the  appropriate  filter  output  was  determined  by 
the  variance  of  the  filters'  smoothed  position  estimates.  The  tri -state 
adaptive  filter  was  implemented  in  the  simulation  of  the  F4J  lateral 
control  system.  The  results  given  in  Part  Four  suggest  that  the  performance 
of  the  F4J  lateral  control  system  may  be  improved  through  the  use  of  a 
tri-state  adaptive  tracking  filter.  Since  the  F4J  longitudinal  control 
system  is  structurally  identical  to  the  lateral  control  system,  the  tri¬ 
state  adpative  filter  may,  in  a  similar  manner,  provide  an  improvement  in 
the  performance  of  the  longitudinal  control  system. 

The  overall  performance  of  the  tri-state  adaptive  tracking  filter 
may  be  enhanced  by  selecting  the  parameters  of  each  of  the  three  compo¬ 
nent  filters  in  such  a  manner  as  to  achieve  a  more  complementary  filter 
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response.  Another  modification  which  might  improve  the  performance  of 
the  tri-state  adaptive  filter  is  the  adjustment  of  the  variance  thresh¬ 
olds  of  the  alpha  and  alpha-beta  filters.  As  was  shown  by  the  results 
of  the  F4J  lateral  control  system  simulation,  the  frequency  response  of 
the  tri -state  adaptive  filter  may  be  altered  by  the  selection  of  the  ap¬ 
propriate  variance  thresholds. 
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ABSTRACT 

A  method  of  estimating  the  centroid  location  of  a  target  utilizing 
scan  return  amplitude  versus  angle  information  is  introduced.  The  method 
is  compared  to  three  thresholding  estimators  and  a  first  moment  estimator 
in  a  computer-simulated  automatic  landing  system. 

It  was  found  that  the  method  introduced  was  the  most  robust  and  ac¬ 
curate  of  the  estimators  in  noise,  due  to  its  unique  scan  rejection  capa¬ 
bility.  In  periods  of  high  signal-to-noise  ratio  the  method  had  less 
error  than  the  thresholding  methods,  and  was  similar  in  ability  to  the 
first  moment  estimator.  Further,  the  pulse  transmissions  required  to  ob¬ 
tain  a  desired  level  of  performance  is  much  reduced  from  the  thresholding 
methods  employed  in  this  simulation. 
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I.  INTRODUCTION 


Essential  to  the  performance  of  any  tracking  radar  is  an  effective 
target  centroid  estimator.  The  purpose  of  this  work  is  to  examine  the  ac¬ 
curacy  of  several  target  centroid  estimators  in  a  comparative  fashion, 
and  to  introduce  a  non-thresholding  algorithm  developed  as  part  of  this 
research.  This  analysis  was  conducted  using  a  software  simulation  of  a 
landing  system  radar  tracking  a  passive  target.  Gaiccari  and  Nucci  [1], 
Shradar  [2],  Mueke  [3],  and  Gilbert  [4],  provide  an  excellant  discussion 
of  air  traffic  control  radars.  The  results  of  this  work  are  most  appli¬ 
cable  to  sequential-lobing  tracking  radars. 


II.  OVERVIEW  OF  THE  SIMULATION 


The  computer  simulation  used  in  this  work  describes  a  large  jet 
fighter  aircraft  in  a  normal  ground  controlled  approach  (GCA)  with  the 
radar  antenna  located  500  meters  from  the  runway  touchdown  point,  as 
shown  in  Figure  2-1.  The  simulation  initiates  the  flight  with  the  tar¬ 
get  3.72  nmi  downrange  from  the  runway  touchdown  point,  or  4.0  nmi  down- 
range  from  the  radar  antenna.  The  target  model  is  allowed  to  approach 

the  runway  at  a  constant  148.6  mph  on  a  3.5  degree  glideslope,  which  is 

»  » 

a  typical  approach  for  the  jet  fighter  being  modelled  [5].  the  radar  is 
a  phased-array  3-0  pencil  beam  radar  utilizing  a  null-to-null  cross-type 
scan,  which  scans  the  target  as  it  moves.  Since  the  tracking  mode  of  an 
operating  radar  attempts  to  find  the  target  within  a  small  area  of  space 
designated  by  the  search  mode,  this  simulation  varies  the  location  of 
the  target  in  the  scanning  window  by  use  of  a  uniform  random  number  gen¬ 
erator  before  the  start  of  each  scan.  The  scanning  window  is  always 
wide  enough  to  fully  scan  the  target. 

The  simulation  executes  a  single  scan  on  the  moving  target  and  then 
Increments  time  to  allow  the  modelled  radar  to  perform  its  other  search 
and  track  duties,  and  to  move  the  target  down  the  gli depath.  The  simula¬ 
tion  aborts  when  the  target  is  within  90  meters  of  the  runway  touchdown 
point. 

The  target  model  used  is  an  ensemble  of  three  anisotropic  scatter¬ 
ing  complexes  representing  the  left  wing,  right  wing,  and  fuselage, 
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Siting  of  the  precision  approach  radar  and  final  approach 
glideslope  used  in  the  computer  simulation. 


slightly  modified  from  the  model  of  Loomis  and  Graf  [6].  The  location  of 
the  scattering  complexes  in  the  target  coordinate  system  is  shown  in  Fig¬ 
ure  2-2(a),  and  the  arrangement  of  the  scattering  points  in  a  scattering 
complex  are  shown  in  Figure  2- 2 ( b ) .  The  equations  describing  the  scatter¬ 
ing  complexes  is  given  in  Table  2-1.  In  this  work,  the  angles  %  and  9 
are  not  the  typical  spherical  phi  and  theta,  but  rather  relative  angles 
measured  from  the  nose  axis  of  the  target  coordinate  system.  Phi  de¬ 
scribes  the  angle  in  azimuth,  and  theta  describes  the  angle  in  elevation. 
Figures  2-3,  2-4,  and  2-5  are  plots  of  the  rade-  cross  section  (RCS)  in 
azimuth  of  the  fuselage,  right  wing,  and  left  wing,  respectively.  The 
composite  cross  sections  of  the  target  model  in  azimuth,  Figure  2-6,  and 
in  elevation,  Figure  2-7,  are  not  used  by  the  simulation,  and  are  pre¬ 
sented  here  for  completeness.  The  radar  cross  sections  in  polar  form  of 
the  fuselage,  right  wing,  and  left  wing,  are  shown  in  Figures  2-8,  2-9, 
and  2-10,  respecti vely.  The  built-in  shadowing  effect  of  the  fuselage 
on  the  wings  is  especially  evident  in  Figures  2-9  and  2-10.  The  compos¬ 
ite  cross  sections  in  azimuth,  Figure  2-11,  and  elevation.  Figure  2-12, 
are  again  shown  for  completeness.  All  figures  are  for  a  wavelength  of 
3.3  cm. 

The  individual  returns  from  each  of  the  scattering  complexes  are 
weighted  by  the  antenna  pattern  before  being  summed  on  a  power  basis. 

This  process  is  repeated  for  every  simulated  transmission  of  a  pulse  from 
the  radar.  Although  only  one  pulse  is  transmitted  at  each  beam  pointing 
location,  time  is  incremented  as  though  t!v.<se  pulses  are  transmitted. 

When  the  simulation  noise  option  is  enabled,  random  gaussian  noise  is 
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Table  2-1.  Radar  cross  section  equations  for  the 
target  model  scattering  complexes 


RCS  equation  for  all  points: 

cr(e,$)  =  A(9,$)|Ax(a)  +  Ay(s)  +  A  2  ( 5 )  |  (m2) 


where: 


kdx 

A  (a)  =  COS  (—  COS  a) 
kd 

A  (5)  =  cos  (-^  cos  6) 
kd 

Az(s)  =  cos  (—  cos  s) 


a,  5,  6  are  assumed 

the  same  for  each  scatterer 

and  are  defined  in  Figure  2-2. 
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Figure  2-3.  Radar  cross  section  of  fuselage  scattering  complex  in 
azimuth,  with  the  azimuth  angle  measured  from  the  nose 
axis  of  the  coordinate  system. 
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Figure  2-5.  Radar  cross  section  of  left  wing  scattering  complex  i 
azimuth,  with  the  azimuth  angle  measured  from  the  nos 
axis  of  the  coordinate  systan. 


Composite  cross  section  in  azimuth,  with  the  azimuth  a 
measured  from  the  nose  axis  of  the  coordinate  system. 
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Figure  2-7.  Composite  cross  section  in  elevation,  with  the  elevation 

angle  measured  from  the  nose  axis  of  the  coordinate  system. 


Figure  2-8.  Radar  cross  section  of  fuselage  scattering  complex  In  azimuth, 
with  the  azimuth  angle  measured  from  the  nose  axis  of  the 
coordinate  system.  Amplitudes  are  in  dB  down  from  maximum. 


Figure  2-10.  Radar  cross  section  of  left  wing  scattering  complex  in 
azimuth,  with  the  azimuth  angle  measured  from  the  nose 
axis  of  the  coordinate.  Amplitudes  are  in  dB  down  from 
maximum. 


16 


180' 


Figure 


90* 


12.  Composite  cross  section  in  elevation,  with  the  elevation 
angle  measured  from  the  nose  axis  of  the  coordinate  sys¬ 
tem.  Amplitudes  are  in  dB  down  from  maximum. 


17 

added  to  the  resultant  return  on  a  power  basis.  The  magnitude  of  the 
noise  power  is  such  that  the  variance  of  the  noise  is  15  dB  down  from  a 
relative  maximum  scan  (without  noise)  at  far  range. 

To  simulate  turbulence,  the  target  coordinate  system  is  allowed 
roll,  pitch,  and  yaw,  with  the  origin  of  the  target  coordinate  system 

? 

locked  on  the  3.5  degree  glideslcpe.  To  simulate  calm  air,  the  target 
model  maintains  a  "wings  level"  attitude  for  the  duration  of  the  flight. 

The  target  returns  are  calculated  with  the  simplified  form  of  the 
radar  equation  [7],  and  are  output  to  the  centroid  estimators.  The  basic 
system  parameters  are  listed  in  Table  2-2. 
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Table  2-2. 


Parameters  for  landing  system  simulation. 


Frequency 

Pulse  repetition  frequency 
Target  initial  elevation 
Target  model  initial  range 
from  touchdown 
Target  model  speed 
Turbulence  rates 


Signal-to-noise  ratio  at  a 
far  range 

Antenna  beamwidth  (null-to-null) 
Azimuth 
Elevation 

Simulation  duration 


9.1  GHz 
6  KHz 

56.6  mrad 

6890  meters 

148.6  mph 

10  deg/s  roll 
5  deg/s  pitch 
5  deg/s  yaw 

15  dB 

1.83° 

1 .77° 

103  scans 


III.  SIGNAL  PROCESSING 


The  computer  simulation  just  described  creates  a  sequence  of  scan 
returns  from  the  target.  In  order  to  neglect  the  effects  of  multipath, 
this  work  will  address  itself  solely  to  that  data  generated  by  the  scan 
in  azimuth.  The  target  centroid  is  calculated  from  the  returns  as  fol¬ 
lows.  A  threshold  determined  from  the  scan  returns  is  applied  to  the 
scan.  Moving  in  from  the  edges  of  the  scan,  the  first  occurrence  of  two 
consecutive  return  voltages  exceeding  the  threshold  is  located.  The  out¬ 
ermost  of  those  return  voltages  are  tagged  as  the  edge-points  of  the  tar¬ 
get.  Since  the  angle  to  the  returns  are  known,  the  centroid  of  the  tar¬ 
get  is  judged  to  be  midway  between  the  edge-points. 

Three  methods  of  setting  the  threshold  are  used  in  this  work.  Two 
are  the  mean,  and  median,  post-determined  thresholds.  That  is,  the  tar¬ 
get  is  scanned  and  the  returns  are  recorded.  The  mean  of  the  scan  returns 
is  calculated,  and  a  threshold  is  set  at  that  level.  Likewise,  the  medi¬ 
an  scan  return  is  found  and  a  threshold  is  set  at  that  level. 

A  third  method  is  a  pre-determi ned  thresholding  method.  The  an¬ 
tenna  beam  is  placed  in  the  center  of  the  scanning  window  to  measure  the 
anticipated  maximum  return  from  that  scan.  The  threshold  is  set  12  dB 
down  from  that  return  level,  which  was  empirically  determined  as  optimal 
with  regard  to  certain  system  model  parameters.  When  two  consecutive  re¬ 
turns  are  above  the  12  dB  threshold,  the  edge  is  marked  and  the  scanning 


translates  to  the  other  side  to  determine  the  other  edge-point.  The 
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requirement  that  the  target  be  fully  scanned  no  longer  exists  for  this 
method,  so  that  fewer  pulses  are  needed  to  locate  the  target  with  little 
or  no  increase  in  estimating  error. 

A  fourth  method  used  is  a  non-thresholding  technique  called  the 
radar  centroid  (RADARCG),  similar  to  that  used  by  Gordon  and  Casowitz 
[8].  This  estimator  weights  each  antenna  pointing  angle  in  the  scanning 
window  by  the  return  from  that  angle,  and  divides  the  sum  of  the  weighted 
angles  by  the  sum  of  the  weights  (returns).  The  result  is  the  angle  to 
the  center  of  gravity  of  the  body  of  the  return.  Since  it  requires  that 
the  window  be  fully  scanned,  all  available  pulses  are  used. 

The  above  methods  form  the  basis  for  comparison  with  the  centroid 
algorithm  based  upon  return  ampl i tude-versus-angl e  signature  introduced 


next. 


IV.  THE  TARGET  CENTROID  ESTIMATING  ALGORITHM 


Introduction 

Since  all  target  centroid  estimators  are  based  on  scan  returns,  it 
is  instructive  to  examine  the  flight  scan-return  history  of  a  target. 
Figure  4-1  is  the  scan  return  history  of  the  model  in  still  air  without 
noise  added,  which  shall  now  be  referred  to  as  a  baseline  flight.  This 
plot  was  made  with  the  target  in  the  center  of  the  scanning  window.  The 
first  and  last  beam  pointing  locations  have  negligible  return  amplitudes 
since  a  null-to-null  cross  track  is  employed;  the  first  null  in  the  an¬ 
tenna  pattern  is  placed  on  the  target  at  those  beam  locations.  As  is  to 
be  expected,  the  maximum  return  occurs  in  the  center  of  the  scan.  It  is 
readily  seen  that  the  scan  returns  over  the  flight  are  modulated,  spec¬ 
ifically  by  the  scintillation  of  the  target  model  radar  cross  section. 

In  particular,  note  scan  number  90.  At  this  scan,  the  antenna  is  clearly 
in  a  null  of  the  target  RCS.  We  can  also  pick  out  scans  78,  and  with 
greater  difficulty,  scan  58,  as  being  in  nulls  of  the  target  model  cross 
section.  It  is  in  these  scans,  with  poor  target  returns,  that  we  would 
expect  the  target  location  error  of  the  estimators  to  increase. 

A  flight  with  noise  is  shown  in  Figure  4-2.  The  two  large  bodies 
of  return  between  scans  58  and  90  are  still  clearly  seen,  but  the  effect 
of  noise  is  pronounced  on  the  rest  of  the  flight.  Beam  pointing  loca¬ 
tions  1  and  49  are  no  longer  at  zero  amplitude,  but  vary  with  noise.  It 

is  clearly  seen  from  observation  of  scans  90,  78,  and  58  that  an  accurate 
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rget  return  history  for  a  baseline  flight  with  noise  added,  49  beam  pointing  locations 
the  scanning  window,  target  in  the  center  of  the  cross  track,  azimuth  scan.  Shown  is 
portion  of  the  flight  from  the  30th  to  the  90th  scan,  inclusive. 
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determination  of  the  presence  of  a  target  at  those  scans  would  be  very 
difficult  and  prone  to  error,  whereas  the  detection  of  the  target  with 
a  good  signal  return,  even  in  the  presence  of  noise,  is  less  prone  to 
error. 

Figure  4-3  is  of  a  baseline  flight  with  turbulence.  The  many  nulls 
in  this  plot  are  the  result  of  the  modulation  of  the  target  model  radar 
cross  section  on  the  target  returns  as  the  model  rotates  on  its  axis  in 
simulation  of  turbulent  wind  conditions.  Again,  beam  locations  1  and  49 
exhibit  negligible  returns  as  the  null  in  the  antenna  pattern  is  on  the 
target. 

Addition  of  noise  to  the  flight  with  turbulence  is  shown  in  Figure 
4-4.  The  many  returns  that  were  of  low  signal  level  are  now  filled  in 
with  noise.  Only  those  scans  whose  signal  level  rises  above  the  noise 
are  suitable  for  target  detection. 

It  is  in  this  light  that  the  work  to  develop  a  new  centroid  algo¬ 
rithm  was  conducted.  The  goal  was  to  produce  an  algorithm  which  would  be 
able  to  determine  which  scans  are  suitable  for  target  detection  and  loca¬ 
tion  -  and  to  discard  all  others. 

Before  proceeding,  a  determination  of  the  expected  scan  signal-to- 
noise  ratios  in  a  typical  flight  is  in  order.  Three  baseline  flights 
with  noise  were  made  with  the  target  at  the  initial  point  3.72  nmi  from 
touchdown  to  the  release  point.  Figure  4-5  is  the  flight  with  a  granu¬ 
larity  of  9  beam  pointing  locations  in  the  scanning  window,  Figure  4-6 
has  a  granularity  of  29,  and  Figure  4-7  has  a  granularity  of  49  beam 
pointing  locations.  It  is  seen  that  the  SNR  amplitude  over  the  flights 
have  the  same  envelope  for  all  granularities.  That  is,  all  are  around  14 


Figure  4-4.  Target  return  history  for  a  baseline  flight  with  noise  and  turbulence  added,  49  beam 
pointing  locations  in  the  scanning  window,  target  in  the  center  of  the  cross  track, 
azimuth  scan.  Shown  is  a  portion  of  the  flight  from  the  30th  to  the  90th  scan,  inclusive 
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fiction  of  scan  number. 
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function  of  scan  number. 


function  of  scan  number. 


30 


dB  until  the  vicinity  of  the  fiftieth  scan  where  the  SNR  drops  a  bit,  due 
to  a  relative  null  in  the  RCS  of  the  target.  From  there  the  SNRs  general¬ 
ly  increase  significantly  as  the  target  draws  nearer  to  the  release  point. 
Note  that  in  all  three  plots  there  are  nulls  at  scan  numbers  58,  78,  and 
90,  which  verify  the  observations  made  earlier  on  Figures  4-1  and  4-2. 
There  are  also  nulls  at  scans  94,  101,  and  103  on  the  plots.  The  fact 
that  the  nulls  appears  on  all  three  plots,  that  "is,  independent  of  granu¬ 
larity,  at  those  scans  is  due  to  a  peculiarity  of  this  simulation.  In 
order  to  reduce  computer  execution  time,  the  target  is  not  scanned  10 
times  a  second  as  in  the  actual  system,  but  rather  the  entire  flight  is 
broken  into  equal  time  units  of  such  a  length  as  to  provide  a  large  num¬ 
ber  of  scans  for  a  statistical  analysis  while  keeping  the  execution  time 
down.  This  was  done  by  scanning  the  target,  moving  it  down  the  glide- 
path,  scanning  it  again,  etc.,  until  the  release  point  was  reached, 
yielding  a  large  dead  time  between  scans.  In  order  to  provide  an  equal 
number  of  statistical  data  points  for  all  flights  regardless  of  granular¬ 
ity  (and  therefore  independent  of  the  number  of  pulses  transmitted),  the 
dead  time  between  scans  is  variable.  It  is  the  greatest  when  the  gran¬ 
ularity  is  9,  and  the  least  when  there  are  49  pulses  to  be  transmitted. 
Each  scan  begins  at  the  same  range  from  touchdown.  Therefore,  range  and 
scan  number  are  related,  and  will  be  used  synonymously  in  this  work.  So 
the  fact  that  nulls  occur  at  scans  58,  78,  90,  94,  101,  and  103,  regard¬ 
less  of  granularity,  is  because  the  target  is  at  the  same  point  in  space 
at  the  beginning  of  the  scan.  The  target  position  at  the  end  of  the  scan 
will  vary  according  to  the  number  of  pulses  that  need  to  be  transmitted 
in  the  scanning  window. 
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Returning  to  Figures  4-5,  4-6,  and  4-7,  it  is  seen  that  all  the 
target  centroid  estimators  must  work  with  scan  signal-to-noise  ratios  of 
around  14  dB  for  at  least  half  of  the  flight. 

Figures  4-8,  4-9,  and  4-10  are  of  a  typical  flight  with  noise  and 
turbulence  for  granularities  of  9,  29,  and  49  beam  pointing  locations, 
respectively.  All  plots  show  a  general  degradation  in  the  SNR  due  to 
the  fluctuating  target  model  RCS  in  turbulence.  Half  the  flight  is  now 
down  to  between  10  and  12  dB,  a  loss  in  signal  strength  of  half  from  the 
flights  without  turbulence. 

The  Scan  Return  Amp! itude-Versus- 
Angle  Signature  Algorithm 

It  is  observed  in  Figure  4-1,  which  graphically  depicts  the  scan 
history  of  a  baseline  flight,  that  all  scan  envelopes  have  a  high  degree 
of  symnetry.  That  is,  as  the  antenna  beam  illuminates  the  target  first 
with  the  pattern  null,  then  increasing  the  illumination  as  the  main  lobe 
moves  onto  the  target,  reaching  the  maximum  when  the  beam  is  centered  or. 
the  target,  then  diminishing  as  the  target  is  placed  in  the  pattern  null, 
the  overall  scan  envelope  takes  on  a  bell  shape  due  to  the  modulation  of 
the  antenna  beam.  Since  the  return  envelopes  are  of  this  shape,  each 
side  of  the  bell  shape  has  a  unique  point,  the  point  of  maximum  slope. 
Returning  to  Figure  4-2  it  is  observed  that  the  maximum  slope  of  a  scan 
with  a  low  SNR  (such  as  scans  58,  78,  and  90)  is  relatively  small,  and 
those  scans  with  large  SNR's  have  a  relatively  large  maximum  slope. 

This,  then,  is  the  chosen  criteria: 

Find  the  point  of  maximum  slope; 


versus  range  for  a  flight  with  noise  and  turbulence 
a  scanning  window  containing  9  beam  pointing  locatii 
ted  as  a  function  of  scan  number. 


IR  versus  range  for  a  flight  with  noise  and  turbulence 
id  a  scanning  window  containing  29  beam  pointing  locat 
otted  as  a  function  of  scan  number. 


Figure  4-10.  SNR  versus  range  for  a  flight  with  noise  and  turbulence 
and  a  scanning  window  containing  49  beam  pointing  loca¬ 
tions,  plotted  as  a  function  of  scan  number. 


Compare  the  slope  at  that  point  to  a  minimum  acceptable 
value  based  on  the  characteristics  of  the  receiver; 

Accept  or  reject  the  scan; 

Determine  target  location  if  the  scan  was  accepted. 

The  method  used  to  find  the  point  of  maximum  slope  is  based  on  the 
scan  shape.  Referring  to  Figure  4-11,  let  us  assume  that  we  are  using  a 
cross  track  with  a  granularity  of  7  beam  pointing  locations.  The  rela¬ 
tive  amplitudes  of  the  expected  returns  are  marked  by  the  lettered  X's 
on  the  drawing.  Moving  from  left  to  right,  the  first  three  returns  have 
a  positive  second  derivative,  since  the  slope  BC  is  greater  than  slope 
AB.  Points  B,  C,  and  D  have  a  negative  second  derivative,  since  slope 
CD  is  less  than  slope  BC.  Since  the  point  of  maximum  slope  is  where  the 
second  derivative  is  zero,  that  is,  where  the  second  derivative  changes 
sign,  the  maximum  slope  must  have  occurred  between  points  B  and  C. 

Having  found  the  maximum  slope,  we  check  to  ensure  that  its  magnitude  is 
greater  than  the  minimum  acceptable  slope.  If  it  is,  the  target  edge  is 
marked  as  being  midway  between  points  B  and  C,  and  scanning  translates  to 
the  other  side  of  the  scan.  The  process  is  then  repeated  for  returns  G, 
F,  E,  and  D.  When  the  two  target  edges  are  found,  the  centroid  is  placed 
midway  between  the  edge  points.  Since  the  target  is  located  by  calculat¬ 
ing  second  derivatives,  this  method  shall  be  referred  to  in  this  work  as 
the  second  derivative  method  or  SDRV. 

A  method  used  to  integrate  the  scan  returns  with  the  second  deri¬ 
vative  method  Is  as  follows.  The  first  half  of  the  scan  is  broken  into 
four  equal  parts,  or  windows,  as  shown  in  Figure  4-12.  For  the  remainder 
of  this  work,  scanning  window  is  an  area  identified  by  the  search  mode 


VOLTAGE  LEVEL 


TARGET  CENTROID  LOCATION 


BEAM  POINTING  LOCATIONS 


Figure  4-11.  Illustration  of  the  method  employed  to  determine  the  tar 
get  centroid  location  based  on  the  shape  of  the  scan 
envelope.  The  signal  returns  are  marked  by  X. 
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Figure  4-12.  Illustration  of  the  pulse  integration  technique  employed 
with  the  second  derivative  algorithm. 
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of  the  radar  which  contains  a  target.  A  section  of  the  scanning  window 
segmented  by  SDRV  shall  be  referred  to  as  a  window  in  this  work.  The 
number  of  pulses  in  each  window  is  determined  by  the  granularity  and  is 
easily  calculated.  The  returns  from  the  beam  pointing  oositions  inside 
the  windows  are  averaged  together,  and  the  location  of  the  averaged  re¬ 
turn  is  placed  in  the  center  of  the  window.  That  is,  the  average  of  the 
five  returns  in  window  1  is  placed  in  the  center  of  the  window  (in  this 
case,  it  has  the  same  angular  location  as  the  third  return)  as  shown  by 
the  average  return  labelled  A.  The  amplitude  of  A  is  not  necessarily 
equal  to  the  center  return  in  the  window.  When  the  four  average  re¬ 
turns  A,  B,  C,  and  D,  are  calculated,  the  determination  of  a  target  edge 
proceeds  as  has  been  previously  described  in  Figure  4-11.  Since  the 
noise  is  gaussian  with  zero  mean,  the  effect  of  noise  will  decrease  as 
more  returns  are  placed  in  the  window  and  averaged  together.  In  the 
event  that  the  averaged  returns  A,  B,  C,  and  D,  do  not  satisfy  the  cri¬ 
terion,  the  windows  are  shifted  one  beam  location  to  the  right,  and  A, 

B,  C,  and  D,  recalculated.  This  process  is  continued  until  a  target 
edge  is  found,  whereby  the  scanning  then  translates  to  the  other  side  of 
the  target.  Otherwise  the  scanning  continues  until  the  fourth  window 
has  shifted  to  within  two  window  widths  of  the  right  side  of  the  scanning 
window,  in  which  case  the  scan  is  rejected  and  further  pulse  transmis¬ 
sions  are  aborted.  In  the  event  that  a  target  edge  was  found,  scanning 
will  proceed  from  right  to  left  until  the  last  beam  location  in  the 
fourth  window  reaches  the  position  calculated  as  the  first  target  edge. 

If  this  occurs,  the  scan  is  rejected  and  further  transmissions  aborted. 
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but  the  area  of  the  last  two  windows  in  the  first  target  edge  was  re¬ 
scanned,  both  using  more  pulses  and  being  doubly  sure  that  no  target  was 
present.  No  centroid  location  decision  is  determined  from  a  rejected 
scan. 

In  this  simulation,  the  minimum  acceptable  rise  of  the  maximum 
slope  in  the  scan  is  related  to  the  amount  of  noise  introduced  into  the 
flight  by  1.5  times  the  standard  deviation  of  the  noise.  This  value  was 
chosen  after  making  many  runs  of  the  simulation  and  observing  the  effect 
of  the  minimum  slope  on  both  the  number  of  rejected  scans  and  on  the  ac¬ 
curacy  of  the  estimator.  The  accuracy  of  ar,  accepted  scan  is  affected 
since,  with  scans  of  poor  SNR,  there  may  be  more  than  one  set  of  returns 
for  which  the  criteria  are  satisfied,  due  to  the  effect  of  noise. 

A  third  and  final  criterion  is  implemented  in  the  algorithm.  With 
no  or  little  signal  present,  it  is  possible  for  return  B  to  be  below  both 
returns  A  and  C,  because  the  probability  of  obtaining  a  negative  value  at 
any  time  from  a  zero  mean  gaussian  process  is  one  half.  If  B  is  a  nega¬ 
tive  quantity,  subtracting  B  from  C  is  a  large  number,  sometimes  greater 
than  the  minimum  slope  criterion.  To  this  end,  return  B  is  first  com¬ 
pared  to  a  voltage  reference.  If  it  is  below  the  reference,  the  windows 
are  shifted  and  B  recalculated.  The  voltage  reference  is  initially  set 
to  zero  at  the  beginning  of  the  flight,  and  is  then  modified  as  follows. 
Each  time  a  set  of  returns  is  rejected,  a  beam  pointing  location  drops 
out  of  window  1  in  Figure  4-12  as  the  windows  shift  to  the  right.  The 
return  from  that  beam  location  is  averaged  into  the  existing  value  of 
the  voltage  reference.  That  is,  the  return  is  added  to  a  register  which 
contains  the  sum  of  all  past  returns  dropped  out  of  the  first  window.  A 


second  register  is  incremented  by  one  to  record  the  number  of  returns  in 
the  sum.  The  average  is  recalculated  for  every  shift  of  the  windows,  for 
the  duration  of  the  flight.  A  voltage  reference  such  as  just  described 
can  be  set  at  a  constant  level  according  to  Ward  [9]  for  a  practical 
radar  receiver. 

In  summary,  the  second  derivative  algorithm  uses  three  criteria. 

In  the  order  the  criteria  are  implemented,  they  are 

1)  (Averaged)  return  B  must  be  greater  than  a  specified 
voltage  reference, 

2)  (Averaged)  returns  A,  B,  and  C  must  have  a  positive 
second  derivative  and  (averaged)  returns  B,  C,  and  D 
must  have  a  negative  second  derivative,  and 

3)  the  maximum  slope  ¥c  must  be  greater  than  the  minimum 
acceptable  slope. 

If  the  scan  returns  satisfy  the  criteria,  a  target  edge  is  assigned  to 
be  midway  between  B  and  C,  and  scanning  translates  to  find  the  second 
target  edge. 

The  first  responsibil ity  of  the  algorithm  is  to  reject  scans  for 
which  an  estimate  of  the  target  position  is  subject  to  severe  error. 

Figure  4-13  is  a  plot  which  illustrates  this  capability.  A  flight  with 
noise  was  flown  with  a  granularity  of  9  beam  locations  in  the  scanning 
window,  and  with  the  antenna  on  the  runway  centerline.  The  antenna  was 
then  moved  to  25  meters  from  the  runway  centerline  and  the  target  reflown. 
The  antenna  position  was  moved  again  by  25  meters,  and  so  on,  ending  with 
the  antenna  located  250  meters  from  the  runway  centerline  for  a  total  of 
11  flights.  In  all  flights,  and  for  all  data  in  the  remainder  of  this 
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Figure  4-13.  Rejection/acceptance  performance  of  the  second  derivative 
algorithm  for  a  granularity  of  9  beam  pointing  locations, 
averaged  over  11  flights  with  noise.  Location  of  the  tar¬ 
get  in  the  scanning  window  is  randomly  varied. 


42 


work,  the  target  location  in  each  scanning  window  was  varied  randomly. 
Figure  4-13  is  the  result  of  those  11  flights.  It  clearly  shows  that  the 
average  number  of  rejected  scans  closely  follows  the  total  number  of 
scans  for  low  scan  SNRs,  and  decreases  as  the  SNR  increases.  There  the 
number  of  accepted  scans  increases  and  matches  the  total  number  of  scans. 
The  crossover  point  is  between  12  and  15  dB. 

Flights  with  noise  and  turbulence  were  flown  as  before  and  results 
plotted  in  Figure  4-14.  There  is  a  higher  average  number  of  scans  with 
both  low  and  high  SNRs  due  to  the  effects  of  turbulence.  These  plots 
are  truncated;  no  scan  SNR  less  than  -9  dB  or  greater  than  50  dB  were 
taken  into  account.  Figure  4-14  is  flatter  than  Figure  4-13,  and  the 
crossover  point  between  the  number  of  accepted  and  rejected  scans  is 
still  about  12-15  dB. 

The  cumulative  number  of  accepted  and  rejected  scans  for  the  flight 
with  noise  is  shown  in  Figure  4-15.  The  rejected  scans  closely  follow 
the  total  number  of  scans  until  around  10  dB,  when  it  begins  to  level 
off  and  the  number  of  accepted  scans  increases.  On  the  average,  59  scans 
are  accepted  and  about  44  scans  rejected  out  of  a  total  of  103  scans, 
which  verify  our  observations  of  Figure  4-5  that  about  half  the  scans 
were  in  low  SNR. 

Figure  4-16  depicts  an  average  flight  with  noise  and  turbulence. 
Here  there  are  fewer  than  sixty  scans  accepted  with  100  scans  evaluated 
(three  of  the  scans  were  below  -9  dB  or  greater  than  50  dB,  and  do  not 
appear  in  the  plot).  It  is  again  evident  that  there  are  more  poor  scans 
with  turbulence  than  without  turbulence,  as  was  observed  in  Figure  4-8. 


PLIGHT  WITH  NOISE  AND  TURBULENCE 
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Figure  4-14.  Rejection/acceptance  performance  of  the  second  derivative 
algorithm  for  a  granularity  of  9  beam  pointing  locations, 
averaged  over  11  flights  with  noise  and  turbulence.  Loca¬ 
tion  of  the  target  in  the  scanning  window  is  randomly  varied. 
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Figure  4-15.  Cumulative  number  of  accepted,  rejected,  and  total  scans 
for  a  flight  with  noise,  granularity  of  9  beam  pointing 
locations. 


45 


Figure  4-16.  Cumulative  number  of  accepted,  rejected,  and  total  scans 
for  a  flight  with  noise  and  turbulence,  granularity  of  9 
beam  pointing  locations. 
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For  a  granularity  of  9  beam  pointing  locations  in  a  scanning  win¬ 
dow  where  the  location  of  the  target  is  varied  randomly,  there  is  only 
one  beam  location  in  each  window.  For  a  granularity  of  29  beam  pointing 
locations,  there  are  two  beam  locations  in  each  window  and  integration 
can  occur. 

Figure  4-17  was  generated  in  a  similar  manner  as  was  Figure  4-13, 
but  with  29  beam  locations.  The  effect  of  pulse  integration,  even  for 
only  two  pulses,  is  already  apparent  by  noticing  that  the  crossover 
point  is  moved  further  down,  to  8  or  9  dB.  Also,  the  average  number  of 
rejected  and  accepted  scans  follow  the  total  number  of  scans  exactly 
below  2  dB  and  above  17  dB,  respectively.  The  flight  with  noise  and 
turbulence,  Figure  4-18,  has  more  scans  at  low  SNR's  cs  does  the  flight 
with  noise  only,  but  the  crossover  remains  at  8-9  dB.  The  advantage  of 
the  integration  is  clearly  seen  in  the  cumulative  number  of  scans  ac¬ 
cepted  and  rejected,  Figures  4-19  and  4-20.  By  comparing  these  figures 
with  Figures  4-16  and  4-17  the  increase  in  the  number  of  accepted  scans, 
and  aecrease  in  the  number  of  rejected  scans,  is  readily  apparent.  But 
reference  to  Figures  4-6  and  4-9  shows  no  real  increase  in  the  number  of 
scans  with  a  high  SNR.  Therefore,  the  increase  in  the  number  of  accepted 
scans  was  due  solely  to  the  integration  process. 

A  granularity  of  49  beam  locations  places  4  beam  pointing  posi¬ 
tions  in  each  window,  permitting  4  returns  to  be  averaged.  Figure  4-21, 
an  average  of  flights  with  noise,  shows  the  crossover  point  moved  back 
to  6  dB.  Virtually  every  scan  above  12  dB  is  accepted,  which  was  the 
crossover  point  for  the  single  return  integration  of  Figure  4-13.  Also 
note  that  tne  crossover  region  becomes  narrower.  In  Figure  4-13,  a  gap 
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algorithm  for  a  granularity  of  29  beam  pointing  locations, 
averaged  over  11  flights  with  noise  and  turbulence.  Loca¬ 
tion  of  the  target  in  the  scanning  window  is  randomly  varied. 
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Figure  4-19.  Cumulative  number  of  accepted,  rejected,  and  total  scans 
for  a  flight  with  noise,  granularity  of  29  beam  pointing 
locations. 
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Figure  4-20.  Cumulative  number  of  accepted,  rejected,  and  total 
for  a  flight  with  noise  and  turbulence,  granularit 
beam  pointing  locations. 


Figure  4-21.  Rejection/acceptance  performance  of  the  second  derivative 
algorithm  for  a  granularity  of  49  beam  pointing  locations 
averaged  over  11  flights  with  noise.  Location  of  the  tar 
get  in  the  scanning  window  is  randomly  varied. 
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of  27  dB  occurs  between  the  points  where  all  scans  are  either  accepted 
or  rejected,  centered  at  13  dB.  In  Figure  4-17,  with  two  returns  aver¬ 
aged  in  a  window,  this  gap  is  reduced  to  12  dB,  centered  at  8  dB.  With 
four  returns  integrated  in  each  window,  the  gap  is  down  to  11  dB  centered 
on  6-7  dB.  In  Figure  4-22,  the  average  of  flights  with  turbulence,  the 
same  observations  can  be  made.  The  crossover  is  lower  than  the  other 
granularities,  now  at  6  dB,  and  virtually  all  scans  above  12  dB  are  ac¬ 
cepted.  As  is  to  be  expected,  the  average  cumulative  number  of  accepted 
scans  for  flights  with  noise,  Figure  4-23,  has  increased,  again  due  to 
the  integration  of  pulses  in  windows.  The  number  of  rejected  scans  also 
is  reduced.  For  flights  with  noise  and  turbulence,  Figure  4-24,  the 
same  is  true.  The  number  of  accepted  scans  has  increased  by  over  20  on 
the  average,  and  the  number  of  rejected  scans  dropped  by  a  like  amount 
in  comparison  to  the  implementation  of  the  second  derivative  algorithm 
without  pulse  integration. 

The  algorithm  having  been  introduced  and  its  rejection  ability 
verified,  the  comparison  of  the  centroid  estimation  accuracy  of  the 
second  derivative  method  to  the  other  estimators  may  now  proceed. 
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Figure  4-23.  Cumulative  number  of  accepted,  rejected,  and  total  scans 
for  a  flight  with  noise,  granularity  of  49  beam  pointing 
locations. 
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V.  COMPARISON  OF  THE  TECHNIQUES 


Each  of  the  five  centroid  estimators  operate  on  the  same  returns 
In  a  different  fashion.  It  would  be  appropriate  to  examine  the  action 
of  each  on  scans  of  different  signal-to-noise  ratios  by  way  of  introduc¬ 
tion  to  a  statistical  analysis  of  the  estimating  error. 

A  plot  of  the  scan  returns  for  a  scan  with  a  SNR  of  21  dB  is  shown 
in  Figure  5-1  with  the  antenna  located  125  m  from  the  runway  centerline. 
The  mean  algorithm  set  the  highest  threshold,  followed  by  the  12  dB  and 
median  thresholds,  all  set  low  on  the  curve.  The  average  returns  inside 
the  windows  of  SDRV  cause  the  windows  to  closely  follow  the  shape  of  the 
curve.  The  set  of  windows  shown  by  the  triangles  which  calculated  the 
first  edge,  shown  by  the  small  left-most  arrow,  are  well  up  on  the  body 
of  the  return  where  the  least  error  should  occur.  The  second  target 
edge  calculated  by  SDRV  occurred  when  the  windows  shifted  to  the  posi¬ 
tions  and  amplitudes  shown  by  the  diamonds.  The  actual  target  location 
was  calculated  during  the  transmission  from  the  22nd  beam  location,  and 
is  denoted  by  the  large  arrow.  The  target  location  calculated  by  the 
radar  center  of  gravity  (RADARCG)  method  is  shown  by  the  large  X.  In 
this  scan,  the  mean,  median,  and  RADARCG  methods  used  all  41  pulses,  12 
dB  used  28,  and  SDRV  used  38.  The  error  ror  this  scan  for  all  five  es¬ 
timators  is  shown  on  the  plot. 

The  second  derivative  method  used  one  fourth  more  of  the  scan  than 

did  the  12  dB  estimator  to  be  sure  that  it  was  on  an  actual  target 
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return.  Since  the  thresholds  are  all  set  relatively  low  on  the  scan,  it 
would  be  easy  for  the  thresholds  to  be  prematurely  tripped  by  noise. 
Comparison  of  Figures  5-2  with  5-3  show  this  to  be  the  case.  Both  the 
median  and  12  dB  estimators  were  satisfied  by  the  noise  outliers  on  the 
left  part  of  the  scan  in  Figure  5-2.  These  were  rejected  by  SDRV  which 
shifted  the  windows  until  it  moved  onto  the  main  body  of  the  return.  Due 
to  the  uniformity  of  the  noise  in  Figure  5-3,  the  12  dB  estimator  proper¬ 
ly  set  the  edges  for  the  most  accurate  estimate  of  the  target  location. 
SDRV  did  move  well  onto  the  return,  but  misjudged  the  location  of  the 
target  by  0.6  mrad  to  the  left.  The  mean  and  median  estimators  made  the 
same  estimate  as  SDRV,  and  the  radar  center  of  gravity  estimator  posi¬ 
tioned  the  target  as  shown  by  the  X  due  to  the  noise  of  the  right  part 
of  the  scan. 

A  scan  of  low  SNR  is  shown  in  Figure  5-4.  Note  the  signal  return 
magnitudes  on  this  plot  in  comparison  to  Figures  5-1,  5-2,  and  5-3.  The 
thresholding  methods  all  set  the  thresholds  low  on  the  scan,  in  the  noise, 
and  came  out  well  due  to  the  uniformity  in  the  noise.  Although  a  main 
body  of  return  appears  obvious  in  this  plot,  it  is  of  such  low  amplitude 
that  it  was  dismissed  as  noise  by  SDRV  and  the  scan  rejected. 

The  comparison  of  the  techniques  is  in  two  basic  parts.  In  the 
first  part,  the  output  of  the  estimators  is  plotted  as  a  function  of  an¬ 
tenna  offset  for  three  different  scanning  granularities.  In  the  second 
part,  the  output  of  the  estimators  is  plotted  as  a  function  of  granular¬ 
ity  for  a  given  antenna  offset.  Both  parts  are  composed  of  the  results 
of  a  baseline  flight,  flight  with  noise,  flight  with  turbulence,  and 
flight  with  noise  and  turbulence,  for  each  data  point,  in  that  order. 


ERROR  RRDRRCG  -0.96  MR  -6.85 
MEAN  0.58  MR  4.13  M  12DB  -8.40  MR  -59.8 
MEDIAN  -6.71  MR  -47.86  M  SDRV  0.02  MR  0.13  M 
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Figure  5-2.  Target  returns  in  the  scanning  window  for  a  granularity  of  41  beam  pointing  locations. 

The  distance  from  target  to  antenna  is  7128  M,  scan  SNR=12.0  dB.  The  location  of  the 
centroid  calculated  by  RADARCG  is  shown  by  the  symbol  X.  The  triangles  and  diamonds 
are  the  SDRV  window  averaged  returns  when  the  target  edges  (small  arrows)  were  found. 
The  large  arrow  is  the  actual  location  of  the  target.  The  error  for  the  scan  is  shown. 


ERROR  RADARCG  1.29  MR  9.14  M 

MEAN  -0.59  MR  -4.20  M  12DB  -0.03  MR  -0.24  M 

MEOIAN  -0.59  MR  -4.20  M  SDRV  -0.59  MR  -4.20  M 
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The  distance  from  target  to  antenna  Is  7062  M,  scan  SNR=12.4  dB.  The  location  of  the 
centroid  calculated  by  RADARCG  is  shown  by  the  symbol  X.  The  tri.  -’es  and  diamonds 
are  the  SDRV  window  averaged  returns  when  the  target  edges  (small  arrows)  were  found. 
The  large  arrow  is  the  actual  location  of  the  target.  The  error  for  the  scan  is  shown. 
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Figure  5-4.  Target  returns  in  the  scanning  window  for  a  granularity  of  41  beam  pointing  locations. 

The  distance  from  target  to  antenna  is  4538  M,  scan  SNR=5.3  dB.  The  location  of  the 
centroid  calculated  by  RADARCG  is  shown  by  the  symbol  X.  No  target  edges  were  found 
by  the  SDRV  algorithm.  The  large  arrow  is  the  actual  location  of  the  target.  The 
error  for  the  scan  is  shown. 
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The  output  of  the  estimators  requires  eight  plots  for  each  set  of  flights, 
one  plot  each  for  the  mean  error,  mean  square  error,  standard  deviation 
of  the  error,  and  the  variance  of  the  error,  in  mill i radians  and  meters. 
Each  data  point  on  each  plot  is  the  result  of  103  scans  (one  flight)  on 
the  target.  The  equations  used  are: 

X  *  (Estimated  value  -  Actual  value) 

Mean  Error  =  jjj-  S  X 

1  2 

Mean  Square  Error  =  g  r  )! 

Standard  Deviation  =  ^  (e*)2) 

Variance  =  jtpj-  (EX^  -  ^  (EX)^) 

where  N  =  Number  of  samples 

The  data  input  to  the  estimators  in  the  baseline  flights  and  flights 
with  noise  are  somewhat  correlated.  After  the  estimators  used  the  scan 
returns  of  the  baseline  flight,  zero  mean  random  gaussian  noise  was  ad¬ 
ded  to  each  scan  return,  and  the  estimators  were  called  again.  That 
data  is  plotted  as  the  flight  with  noise.  In  the  same  way,  the  data  in¬ 
put  to  the  estimators  for  the  flights  with  turbulence  and  noise  were 
made  from  the  scan  returns  of  the  flights  with  turbulence. 

In  the  next  thirty  two  plots,  all  scans  were  used  In  the  error 
analysis  irregardless  of  the  scan  SNR.  A  scanning  granularity  of  9  beam 
pointing  locations  was  employed.  The  plots  of  the  baseline  flights  ver¬ 
sus  antenna  location.  Figures  5-5  through  5-12,  all  show  the  RADARCG 
estimator  to  have  the  least  error,  with  the  SDRV  estimator  almost  as 
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Figure  5-7.  Standard  deviation  of  error  of  estimators  In  milliradians 
for  a  granularity  of  9  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Variance  of  error  of  estimators  in  milllradians  for  a  grai 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-9.  Mean  error  of  estimators  in  meters  for  a  granularity  of  9 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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Figure  5-10.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-11.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  9  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-12.  Variance  of  error  of  estimators  in  meters  for  a  granular 
of  9  beam  pointing  locations.  Each  data  point  is  the  re 
suit  of  one  flight,  all  scans  used. 
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accurate.  The  three  thresholding  methods  are  very  nearly  equal,  but  com 
parison  of  the  plots  will  show  that  the  12  dB  estimator  has  the  least 
error,  then  the  mean  estimator,  and  finally  the  median  estimator.  Under 
these  no-noise,  non-turbulent  conditions,  there  does  not  appear  to  be 
any  appreciable  error  introduced  by  siting  the  antenna  250  meters  from 
the  runway  centerline  as  compared  to  being  on  the  runway  itself. 

The  flights  with  noise,  Figures  5-13  through  5-20,  begin  to  show 
the  relative  merits  of  the  estimators.  Here  the  second  derivative  algo¬ 
rithm  clearly  has  the  least  error  in  noise  alone,  less  than  1.5  mrad  or 
around. 4  meters  on  the  average  according  to  Figures  5-13  and  5-17.  It 
is  also  the  most  stable,  as  shown  by  Figures  5-15,  5-16,  5-19,  and  5-20. 
Whereas  RADARCG  was  excellent  without  noise,  the  figures  indicate  that 
it  it  significantly  degraded  in  noise.  The  thresholding  estimators  are 
still  very  close  to  each  other,  with  the  same  approximate  order  of  error 
as  in  the  baseline  flights.  While  all  estimators  were  degraded  with  the 
introduction  of  noise,  there  is  still  no  apparent  effect  due  to  antenna 
offset. 

The  flights  with  turbulence  only,  Figures  5-21  through  5-28,  show 
little  error  difference  in  comparison  to  the  corresponding  plots  of  the 
baseline  flights,  Figures  5-5  through  5-12.  This  illustrates  that  for  a 
low  order  granularity,  the  effect  of  turbulence  on  the  estimators  is 
smal 1 . 

The  flights  with  noise  and  turbulence,  Figures  5-29  through  5-36, 
again  show  the  second  derivative  method  to  be  the  most  accurate  and  most 
robust  estimator  for  a  scan  with  nine  beam  locations  in  noise  and  turbu¬ 
lence.  The  thresholding  methods  still  appear  to  be  about  equal  in 
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Figure  5-14.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 


uw  ‘yoyyj  jo  nouhiajq  oybONHi§ 


00  *8  00*9  00  *h  00*2  00*Od 


<o 

</>  -4-J 
c  ra 
<0  T3 

•F" 

■a  .r: 

<0  u 

W  ro 

•p*  LU  • 
p-  T3 

*—  a> 

•t“  •  SA 

E  w  3 

C  o  cn 
*p»  *f-  c 
4->  r© 
(/)  (TJ  U 
L.  U  V) 
O  O 

■*->  r—  »— 

tO  •— 

E  cn  fo 


w  -p  -p 

oj  c  -c 

•p  cr» 
4-  o  *»“ 
O  Q.*— 

O  H  03 

I-  a)  c 

s_  .3  o 

03 

CTi  4- 
4-  O 
O  >+- 
O  +■» 
C  i— 
O  >>  3 

■r  +)  1/1 

■M  t-  <U 

ig  L  I. 

•r-  nj 
>  r—  <u 

a)  3  -e 
■O  C  -t-> 
IT} 

13  !.  U1 
L.  U)*r- 
<0 

c  c 
t-  •*- 
**->  o  o 

(/)  4-  fl. 


I 

in 


a> 


75 


Figure  5-16.  Variance  of  error  of  estimators  in  mllliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-17.  Mean  error  of  estimators  in  meters  for  a  granularity  of  9 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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Figure  5-18.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-19.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  9  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-20.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-22.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-23.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  9  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-24.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-25.  Mean  error  of  estimators  in  meters  for  a  granularity  of  9 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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Figure  5-26.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-27.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  9  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-28.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-29.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-30.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-32.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 


Figure  5-33.  Mean  error  of  estimators  in  meters  for 
beam  pointing  locations.  Each  data  po 
one  flight,  all  scans  used. 
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Figure  5-35.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  9  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-36.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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quality,  with  a  slight  preference  to  12  dB.  RADARCG  has  now  deteriorated 
to  the  least  accurate  and  least  robust  estimator,  due  to  the  low  data 
rate  which  does  not  allow  the  algorithm  to  properly  decouple  the  noise. 

While  the  above  graphs  included  all  scans,  it  would  be  worthwile 
to  examine  the  performance  of  the  estimators  with  signal -to-noise  ratios 
of  13  dB  or  greater.  Only  those  scans  were  selected  from  the  data  and 
plotted  in  the  next  series  of  figures. 

The  flights  with  noise,  all  scans  above  13  dB,  Figures  5-37  through 
5-44,  indicate  that  SDRV  and  RADARCG  are  both  approximately  equal  in  per¬ 
formance.  SDRV  is  slightly  favored  since  it  appears  to  be  more  robust 
at  the  various  antenna  offsets.  The  thresholding  estimators  are  of  equal 
qual ity. 

With  turbulence.  Figures  5-45  through  5-52,  the  same  is  true. 

While  it  would  appear  that  SDRV  is  slightly  less  accurate  than  RADARCG 
by  viewing  the  plots  with  the  error  in  milliradians.  Figures  5-45  through 
5-48,  SDRV  was  more  accurate  in  the  actual  meter  error  from  the  target, 
Figures  5-48  through  5-52.  This  error  occurred  when  the  target  was  at 
close  range,  when  a  large  error  in  milliradians  is  a  small  error  in 
meters.  The  thresholding  techniques  are  comparable  to  each  other  and 
are  not  as  accurate  as  SDRV  or  RADARCG. 

Scans  which  had  signal-to-noise  ratios  of  10  dB  or  less  reflect 
the  ability  of  the  estimators  to  find  the  target  in  noise.  The  flights 
with  noise,  Figures  5-53  through  5-60,  show  the  second  derivative  method 
to  be  the  most  accurate  with  the  least  deviation  or  variance  in  noise. 
RADARCG  has  the  greatest  error  of  the  estimators.  The  thresholding 
methods  are,  again,  of  similar  quality. 
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Figure  5-37.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
used. 
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Figure  5-39.  Standard  deviation  of  error  of  estimators  in  mi  1 1 i rad ia ns 
for  a  granularity  of  9  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  13 
dB  or  greater  used. 
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Figure  5-40.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  poincing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-41.  Mean  error  of  estimators  in  meters  for  a  granularity  of  9 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  13  dB  or  greater  used. 
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Figure  5-42.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater- 
used. 
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Figure  5-43.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  9  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-44.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
used. 
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Figure  5-45.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
used. 
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Figure  5-46.  Mean  square  error  of  estimators  in  mill iradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-47.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  9  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  13 
dB  or  greater  used. 
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Figure  5-48.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-49.  Mean  error  of  estimators  in  meters  for  a  granularity  of  9 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  13  dB  or  greater  used. 
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Figure  5-51.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  9  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-52.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
used. 
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Figure  5-53.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used . 
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Figure  5-55.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  9  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  at 
or  below  10  dB  used. 
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Figure  5-57.  Mean  error  of  estimators  in  meters  for  a  granularity  of  9 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  at  or  below  10  dB  used. 
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Figure  5-58.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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Figure  5-59.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  9  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  at  or 
below  10  dB  used. 
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Figure  5-60.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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The  flights  with  turbulence  and  noise,  Figures  5-61  through  5-68, 
reflect  the  same  observations  made  earlier.  It  is  interesting  to  note 
that  with  turbulence  and  noise,  Figure  5-61,  the  second  derivative  esti¬ 
mator  is  as  accurate  as  the  thresholding  methods  were  with  an  infinite 
SNR  and  no  turbulence,  Figure  5-5.  This  clearly  shows  the  increased 
ability  of  SDRV  to  track  targets  over  the  thresholding  methods. 

A  granularity  of  9  beam  pointing  locations  forced  the  second  deriv¬ 
ative  method  to  have  only  one  beam  location  in  each  window.  No  pulse 
integration  was  possible.  With  29  beam  locations  in  the  scanning  window, 
each  window  contains  two  beam  locations,  and  some  noise  rejection  action 
will  occur  due  to  the  averaging  of  the  returns  in  each  window. 

The  results  of  a  baseline  flight  versus  antenna  offset  for  a  scan¬ 
ning  window  with  29  beam  locations  are  shown  in  Figures  5-69  through 
5-76.  The  second  derivative  method  has  the  least  error,  followed  closely 
by  RADARCG ,  followed  by  the  thresholding  estimators.  Of  the  thresholding 
methods,  the  12  dB  method  has  the  least  deviation  and  variance,  shown  in 
Figures  5-71,  5-72,  5-75,  and  5-76,  and  thus  distinguishes  itself  from 
the  mean  and  median  methods.  Comparison  of  these  plots  with  the  base¬ 
line  flight  plot  of  9  beam  pointing  locations,  specifically  comparing 
Figures  5-69  and  5-73  with  Figures  5-5  and  5-9,  shows  that  whereas  the 
thresholding  methods  and  RADARCG  have  no  real  change  in  their  estimating 
ability  due  to  the  increased  amount  of  data,  the  error  of  the  SDRV  algo¬ 
rithm  was  cut  by  approximately  two-thirds. 

The  flights  with  noise,  Figures  5-77  through  5-84,  show  the  SDRV 
method  to  be  the  most  accurate  as  well  as  the  most  robust  estimator. 
RADARCG  improved  significantly  with  the  increased  rate,  and  the  thresh¬ 
olding  algorithms  are  still  comparable,  with  the  12  dB  estimator  leading. 
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Figure  5-61 .  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used . 
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Figure  5-62.  Mean  square  error  of  estimators  in  mill iradians  for  a  gran¬ 
ularity  of  S  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
10  dB  used. 
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Figure  5-64.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  9  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
10  dB  used. 
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Figure  5-65.  Mean  error  of  estimators  in  meters  for  a  granularity  of  9 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  at  or  below  10  dB  used. 
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Figure  5-66.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used . 
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Figure  5-68.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  9  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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Figure  5-69.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-70.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-72.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-73.  Mean  error  of  estimators  in  meters  for  a  granulatity  of  29 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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Figure  5-74.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-75.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  29  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-76.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re- 
sule  of  one  flight,  all  scans  used. 
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Figure  5-77.  Mean  error  of  estimators  in  mi  1 1 i radians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-79.  Standard  deviation  of  error  of  estimators  in  milliraoians 
for  a  granularity  of  29  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-80.  Variance  of  error  of  estimators  in  mill iradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 


BPSELINE  FLIGHT  WITH  NOISE 


Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  29  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-84.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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The  results  of  the  flights  with  turbulence,  Figures  5-85  through 
5-92,  show  results  similar  to  the  baseline  flights,  just  as  the  flights 
with  turbulence  were  similar  to  the  baseline  flights  with  a  scanning 
granularity  of  9  beam  locations.  This  is  due  to  the  lack  of  noise  in 
the  scan  returns. 

The  flights  with  noise  and  turbulence  for  all  scans,  shown  in  Fig¬ 
ures  5-93  through  5-100,  indicates  a  slight  degradation  of  estimating 
quality  in  the  thresholding  methods  from  the  baseline  flights.  The  sec¬ 
ond  derivative  method  was  degraded,  but  clearly  remains  the  most  accu¬ 
rate  of  the  estimators.  The  mean  error  in  milliradians  of  the  radar 
centroid  algorithm  indicates  that  the  mean  error  is  less  than  that  of 
the  thresholding  methods,  but  the  mean  square  error  in  milliradians  and 
meters,  Figures  5-94  and  5-98,  and  the  mean  error  in  meters.  Figure  5-97, 
are  in  excess  of  the  thresholding  methods. 

Selecting  only  those  scans  which  are  equal  to,  or  greater  than  13 
dB,  the  estimating  ability  of  the  methods  was  examined.  Figures  5-101 
through  5-108  are  the  results  of  the  flights  with  noise,  scans  13  dB 
and  greater  only.  As  was  to  be  expected,  RADARCG  is  the  best  estimator 
in  periods  of  high  SNR,  and  the  thresholding  estimators  are  of  equal 
quality.  With  turbulence  effects,  Figures  5-109  through  5-116,  all  esti¬ 
mators  are  degraded  moderately.  As  in  the  flights  with  noise,  RADARCG 
is  the  most  accurate  estimator,  followed  by  SDRV,  with  the  thresholding 
techniques  of  lower  quality. 

The  accuracy  of  the  estimators  using  scans  with  SNRs  of  10  dB  or 
less  in  flights  with  noise  are  depicted  in  Figures  5-117  through  124. 

The  second  derivative  appears  to  be  the  best  of  the  estimators  in  noise, 
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Figure  5-85.  Mean  error  of  estimators  in  inilliradians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-87.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  29  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-88.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-89.  Mean  error  of  estimators  in  meters  for  a  granularity  of  29 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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figure  5-90.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-91.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  29  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-92.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 


Figure  5-93.  Mean  error  of  estimators  in  mil  1 iradians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-95.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  29  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-96.  Variance  of  error  of  estimators  In  mllllradlans  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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quare  error  of  estimators  in  meters  for  a  granularity 
beam  pointing  locations.  Each  data  point  is  the  re- 
if  one  flight,  all  scans  used. 
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Figure  5-99.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  29  beam  pointing  locations.  Each  data  point 
is  the  resist  of  one  flight,  all  scans  used. 
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Figure  100.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-101.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  5NR  13  dB  or  greater 
used. 
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Figure  5-103.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  29  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  13 
dB  or  greater  used. 


Figure  5-104.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-105.  Mean  error  of  estimators  in  meters  for  a  granularity  of  29 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  13  dB  or  greater  used. 
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Figure  5-106.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
used. 
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Figure  5-107.  Standard  deviation  of  error  of  estimators  In  meters  for  a 
granularity  of  29  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-108.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  of  13  dB  or  greater 
used. 
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Figure  5-109.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
used. 
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Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-111.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  29  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  13 
dB  or  greater  used. 
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Figure  5-113.  Mean  error  of  estimators  in  meters  for  a  granularity  of  29 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  13  dB  or  greater  used. 
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Figure  5-114.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  dB  or  greater 
used. 


FLIGHT  WITH  TURBULENCE  RND  NOISE 
O  MEAN  +  RAOARCG  O  SDRV 

A  MEDIAN  X  12DB 


177 


Figure  5-115.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  29  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-117.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10 
dB  used. 
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Figure  5-119.  Standard  deviation  of  error  of  estimators  in  mi  1 1 iradians 
for  a  granularity  of  29  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  at 
or  below  10  dB  used. 
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Figure  5-121.  Mean  error  of  estimators  in  meters  for  a  granularity  of 

29  beam  pointing  locations.  Each  data  point  is  the  result 
of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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granularity  of  29  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  at  or 
below  10  dB  used. 
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Figure  5-124.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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Figure  5-124.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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due  to  the  method  of  rejecting  scans  which  do  not  satisfy  the  criteria. 
RADARCG  has  a  lower  mean  and  mean  square  error  than  the  thresholding 
methods,  but  it  is  not  robust  in  noise.  Comparison  to  Figures  5-53 
through  5-60  discloses  a  general  improvement  in  all  estimators,  due  to 
the  increased  data  provided  by  a  scan  of  finer  granularity. 

The  error  due  to  scans  with  signal -to-noise  ratios  of  10  dB  or  low¬ 
er  in  the  flights  with  turbulence  and  noise,  Figures  5-125  to  5-132,  is 
on  the  order  of  the  flights  without  turbulence.  The  thresholding  esti¬ 
mators  show  no  additional  degradation  due  to  turbulence.  However,  both 
RADARCG  and  SDRV  experience  an  increase  in  error,  particularly  with  the 
antenna  at  a  large  distance  from  the  runway.  SDRV  is  still  the  most  ac¬ 
curate  of  the  estimators,  but  RADARCG  is  greatly  affected  by  the  noise 
and  produces  estimates  of  the  target  location  with  greater  error  than 
the  thresholding  methods. 

Having  discussed  the  error  produced  by  the  estimators  with  scanning 
granularities  which  permit  the  averaging  of  1  and  2  beam  location  returns 
in  its  windows,  the  discussion  of  a  still  finer  scanning  granularity  is 
in  order.  With  49  beam  locations  in  the  scanning  window,  the  returns 
from  four  beam  pointing  locations  are  averaged  in  each  window  of  the  sec¬ 
ond  derivative  algorithms. 

The  baseline  flights,  Figures  5-133  through  5-140,  show  little  im¬ 
provement  in  the  mean  or  mean  square  error  for  the  thresholding  methods 
or  RADARCG.  The  standard  deviation  and  variance  of  the  thresholding 
methods  did  improve,  causing  their  plots  to  closely  coincide.  The  sec¬ 
ond  derivative  method  displays  an  improvement  in  estimating  capability, 
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Figure  5-125.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  29  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10 
dB  used. 
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Figure  5-126.  Mean  square  error  of  estimators  in  mill iradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
10  dB  used. 
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Figure  5-127.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  29  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  at 
or  below  10  dB  used. 
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Figure  5-128.  Variance  of  error  of  estimators  in  mi  1 1 iradians  for  a  gran¬ 
ularity  of  29  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
10  dB  used. 
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Figure  5-129.  Mean  error  of  estimators  in  meters  for  a  granularity  of  29 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  at  or  below  10  dB  used. 
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FLIGHT  WITH  TURBULENCE  AND  NOISE 


Standard  deviation  of  error  of  estimators  in  meters  for 
granularity  of  29  beam  pointing  locations.  Each  data  pi 
is  the  result  of  one  flight,  all  scans  with  a  SNR  at  or 
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Figure  5-132.  Variance  of  error  of  estimators  in  meters 
of  29  beam  pointing  locations.  Each  data 
suit  of  one  flight,  all  scans  with  a  SNR  , 
used. 
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rigure  5-133.  Mean  error  of  estimators  in  milliradians  for  a  granularit 
of  49  beam  pointing  locations.  Each  data  point  is  the  re 
suit  of  one  flight,  all  scans  used. 
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square  erro»  of  estimators  in  milliradians  for  a  gran- 
ty  of  49  beam  pointing  locations.  Each  data  point  is 
esult  of  one  flight,  all  scans  used. 
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Figure  5-135.  Standard  deviation  of  error  of  estimators  in  milliradia 
for  a  granularity  of  49  beam  pointing  locations.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-137.  Mean  error  of  estimators  in  meters  for  a  granularity  of  49 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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Figure  5-138.  Mean  square  error  of  estimators  in  meters  for  a  gramT  ritj 
of  49  beam  pointing  locations.  Each  data  point  is  the  r* 
suit  of  one  flight,  all  scans  used. 


Figure  5-140.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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and  is  the  most  robust  of  the  estimators  with  regard  to  antenna  offset 
location  versus  actual  meter  error. 

The  flights  with  noise,  Figures  5-141  through  5-148,  show  an  im¬ 
provement  in  all  estimators  over  the  29  beam  location  scans,  RADARCG  in 
particular.  The  standard  deviation  and  variance  of  the  errors  are  more 
stable  than  in  the  coarser  granularities,  due  to  the  increase  in  data 
points  upon  which  a  location  decision  can  be  made. 

The  baseline  flights  with  turbulence,  Figures  5-149  through  5-156, 
show  little  change  over  those  baseline  flights  without  turbulence  (Fig¬ 
ures  5-133  to  5-140).  There  is  a  curious  jump  in  the  standard  deviation 
and  variance  of  the  SDRV  milliradian  error  at  antenna  locations  of  225 
and  250  meters  from  the  runway,  Figures  5-151  and  5-152,  but  there  is 
only  the  slightest  deviation  from  a  straight  line  at  those  antenna  loca¬ 
tions  in  the  plots  of  the  standard  deviation  and  variance  of  the  error 
in  meters.  Figures  5-155  and  5-156.  The  error  probably  occurred  at  close 
range. 

The  flights  with  turbulence  and  noise,  Figures  5-157  through  5-164, 
indicate  that  SDRV  is  again  the  estimator  with  the  least  error  in  turbu¬ 
lence  and  noise.  There  is  an  improvement  in  all  estimators  over  the 
coarser  granularities,  particularly  in  the  case  of  RADARCG.  T'ne  estima¬ 
tor  is  much  more  stable  in  the  higher  data  rates,  but  the  standard  devia¬ 
tion  and  variance  are  still  in  excess  of  those  of  the  thresholding  meth¬ 
ods. 

As  before,  it  would  be  advantageous  to  examine  the  performance  of 
the  estimators  at  the  increased  data  rate  in  both  high  and  low  scan 
SNR's. 
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Figure  5-141.  Mean  error  of  estimators  in  mi  11 iradians 
of  49  beam  pointing  locations.  Each  da t. 
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Figure  5-142.  Mean  square  error  of  estimators  in  mi  1 1 iradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-143.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  49  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-144.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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lean  error  of  estimators  in  meters  for  a  granularity 
ieam  pointing  locations.  Each  data  point  is  the  resu 
me  flight,  all  scans  used. 
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Figure  5-147.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  49  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-148.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-149.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-150.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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Figure  5-151.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  49  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-153.  Mean  error  of  estimators  in  meters  for  a  granularity  of  49 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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Figure  5-154.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-155.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  49  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-156.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figure  5-157.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 


Figure  5-158.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 
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igure  5-159.  Standard  deviation  of  error  of  estimators  in  mill iradians 
for  a  granularity  of  49  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-161.  Mean  error  of  estimators  in  meters  for  a  granularity  of  49 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  used. 
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Mgure  5-164.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  used. 
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Figures  5-165  through  5-172  are  the  results  of  the  calculation  of 
the  mean  error,  mean  square  error,  and  the  standard  deviation  and  vari¬ 
ance  of  the  error  in  milliradians  and  meters  using  only  the  data  from 
those  scans  of  a  flight  with  a  SNR  of  at  least  13  dB.  The  figures  show 
the  RADARCG  estimator  to  be  the  most  accurate,  in  contrast  the  SDRV  was 
most  accurate  when  the  scan  SNR  was  infinite  (baseline  flights).  The 
three  thresholding  methods  are  all  comparable,  with  a  slight  edge  given 
to  the  mean  estimator. 

The  addition  of  turbulence,  Figures  5-173  through  5-180,  shows  lit¬ 
tle  change  from  the  flights  without  turbulence.  This  again  indicates 
that  little  decorrelation  occurs  from  pulse  to  pulse  due  to  turbulence 
for  this  simulation. 

The  scans  less  than  or  equal  to  a  SNR  of  10  dB  for  a  flight  with 
noise  were  selected  and  the  data  plotted  in  Figures  5-181  through  5-188. 
Most  noticeable  is  the  improvement  in  RADARCG.  The  algorithm,  now  having 
more  data  points  to  assimilate,  is  working  well  in  reducing  the  effects 
of  noise.  The  most  accurate  and  robust  estimator  is  still  SDRV,  although 
the  choise  is  not  as  obvious  as  it  was  in  the  coarser  granularities.  It 
is  interesting  to  note,  by  comparing  Figure  5-53  with  Figure  5-181  for 
example,  that  the  mean  error  in  milliradians  of  SDRV  for  nine  beam  loca¬ 
tions,  scans  less  than  10  dB,  is  at  least  as  accurate  as  are  the  thresh¬ 
olding  estimators  under  those  same  conditions,  but  using  49  beam  loca¬ 
tions. 

The  flights  with  turbulence  and  noise,  scans  at  or  below  10  dB,  are 
plotted  in  Figures  5-189  through  5-196.  All  estimators  are  improved  over 
similar  conditions  with  the  coarser  granularities.  Note  that  the  second 
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Figure  5-166.  Mean  square  error  of  estimators  in  mil  1 iradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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variance  or  error  or  estimators  in  miliiradians  for  a  gran 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-169.  Mean  error  of  estimators  in  meters  for  a  granularity  of  49 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  13  dB  or  greater  used. 
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Figure  5-170.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
used. 
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Figure  5-171.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  49  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-174.  Mean  square  error  of  estimators  in  mi  1 1 iradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-175.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  49  beam  pointing  locations.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  13  dB  or  greater  used. 
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Figure  5-176.  Variance  of  error  of  estimators  in  milliradians  for  a  gran 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-177.  Mean  error  of  estimators  in  meters  for  a  granularity  of  49 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  13  dB  or  greater  used. 
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Figure  5-178.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
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Figure  5-180.  Variance  of  error  of  estimators  in  meters  for  a  grar.ilari 
of  49  beam  pointing  locations.  Each  data  point  is  the  re 
suit  of  one  flight,  all  scans  with  a  SNR  13  dB  or  greater 
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Figure  5-182.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
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Figure  5-183.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  49  beam  pointing  locations.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  at  or  below  10  dB  used. 


BASELINE  FLIGHT  WITH  NOISE 
0  MEAN  +  RADARCG  <!>  SDRV 

a  MEDIAN  X  1 2DB 


249 


Figure  5-184.  Variance  of  error  of  estimators  in  mi  1 1 iradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
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Figure  5-185.  Mean  error  of  estimators  in  meters  for  a  granularity  of  49 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SHR  at  or  below  10  dB  used. 
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Figure  5-186.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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Figure  5-187.  Standard  deviation  of  error  of  estimators  in  meters  for  a 
granularity  of  49  beam  pointing  locations.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  at 
or  below  10  dB  used. 
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Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 


Figure  5-189.  Mean  error  of  estimators  in  milliradians  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
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Figure  5-190.  Mean  square  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
10  dB  used. 
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Figure  5-191.  Standard  deviation  of  error  of  estimators  in  milliradians 
for  a  granularity  of  49  beam  pointing  locations.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  at  or  below  10  dB  used. 


FLIGHT  WITH  TURBULENCE  RND  NOISE 
O  MEAN  +  RRDflRCG  <!>  SDRV 

A  MEDIAN  X  12DB 


UW  'UQUU3  30  33NUIUUA  ° 
00*091  00*021  00*09  00*0h  00*0= 


UW  U0UU3  30  30NHIUUA 


Figure  5-192.  Variance  of  error  of  estimators  in  milliradians  for  a  gran¬ 
ularity  of  49  beam  pointing  locations.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  at  or  below 
10  dB  used. 
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Figure  5-193.  Mean  error  of  estimators  in  meters  for  a  granularity  of  49 
beam  pointing  locations.  Each  data  point  is  the  result  of 
one  flight,  all  scans  with  a  SNR  at  or  below  10  dB  used. 
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Figure  5-194.  Mean  square  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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Figure  5-195.  Standard  deviation  of  error  of  estimators  in  meters  for  a 

granularity  of  49  beam  pointing  locations.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  at  or  be¬ 
low  10  dB  used. 
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Figure  5-196.  Variance  of  error  of  estimators  in  meters  for  a  granularity 
of  49  beam  pointing  locations.  Each  data  point  is  the  re¬ 
sult  of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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derivative  estimator,  using  only  those  scans  at  or  below  10  dB,  is  still 
more  accurate  than  the  thresholding  methods  in  baseline  flights,  Figures 
5-133  through  5-140, 

Since  there  was  little  or  no  effect  of  antenna  location  on  the  er¬ 
ror  of  the  estimators,  we  shall  consider  the  effect  of  scan  granularity 
on  error  for  only  one  antenna  location  which  shall  be  125  meters  from 
the  runway  centerline.  For  the  next  32  figures,  all  scans  will  be  in¬ 
cluded  in  the  analysis  regardless  of  scan  SNR. 

The  errors  plotted  as  a  function  of  granularity  for  baseline  flights 
are  shown  in  Figures  5-197  through  5-204.  The  most  surprising  aspect  of 
these  plots  is  that  there  is  no  decrease  in  error  with  the  increase  in 
scanning  granularity,  except  for  SDRV  which  generally  improves  with  the 
finer  granularities.  It  would  then  appear  that,  on  the  average,  the  lim¬ 
it  of  accuracy  for  the  thresholding  methods,  RADARCG,  and  SDRV  is  3  mr, 

0.4  mrad,  and  0  mrad,  respectively,  as  shown  by  Figure  5-197.  The  stan¬ 
dard  deviation  and  variance  also  show  little  change  with  the  exception 
of  SDRV,  which  decreases. 

Figures  5-205  through  5-212  are  the  results  of  the  baseline  flights 
with  noise.  The  thresholding  estimators  do  not  increase  in  accuracy  with 
an  increase  in  granularity,  but  rather  converge  on  a  value  of  error.  Both 
the  second  derivative  and  RADARCG  improve  with  finer  granularities.  The 
second  derivative  has  an  increase  in  error  between  the  granularities  of 
9  and  23,  when  there  is  only  one  beam  location  in  each  window.  Two  beam 
locations  are  in  each  window  with  granularities  of  23  through  29,  three 
beam  locations  are  used  from  31  to  45,  and  granularities  of  47  and  49  use 
four  beam  locations  in  each  window.  Note  that  the  error  is  a  relative 
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Mean  error  of  estimators  in  milliradians  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-198.  Mean  square  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-199.  Standard  deviation  of  error  of  estimators  in  milliradians 
with  the  antenna  located  125  M  from  the  runway  centerline 
Each  data  point  is  the  result  of  one  flight,  all  scans 


Figure  5-200.  Variance  of  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  th*.?  result  of  one  flight,  all  scans  used. 
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Figure  5-201.  Mean  error  of  estimators  in  meters  with  the  antenna  located 
125  M  from  the  runway  centerline.  Each  data  point  is  the 
result  of  one  flight,  all  scans  used. 


Figure  5-202.  Mean  square  error  of  estimators  in  meters  v/ith  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-204.  Variance  of  error  of  estimators  In  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-207.  Standard  deviation  of  error  of  estimators  in  milliradians 
with  the  antenna  located  125  M  from  the  runway  centerline. 
Each  data  point  is  the  result  of  one  flight,  all  scans 
used. 
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Figure  5-208.  Variance  of  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-209.  Mean  error  of  estimators  in  meters  with  the  antenna  loca¬ 
ted  125  M  from  the  runway  centerline.  Each  data  point  is 
the  result  of  one  flight,  all  scans  used. 


BASELINE  FLIGHT  WITH  NOISE 
©  MEAN  +  RADARCG  O  SDRV 

A  MEDIAN  X  1 2DB 


276 


tO 

00 '08 


ix  sy3i3w  4yoyy3  3ybnos 

00 ‘09  00  ’Oh  00*02 

*  1  l  i  1  l  i  1  -i  1  .1— *  1  1  i  1  I  1  1— i.  i  1  1  t  1  1  t  i  Ufa  1  .A  I 


NB^J 

00*0 


CO 

lO 


Figure  5-210.  Mean  square  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-211.  Standard  deviation  of  error  of  estimators  in  meters  with 
the  antenna  located  125  M  from  the  runway  centerline. 
Each  data  point  is  the  result  of  one  flight,  all  scans 
used. 
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Figure  5-212.  Variance  of  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 


minimum  at  granularities  of  9,  23, tff,  and  47.  This  implies  that  the 
coarsest  granularity  for  the  windows  to  hold  a  desired  number  of  beam 
locations  should  be  used  to  minimize  error.  This  result  is  to  be  expect¬ 
ed.  The  second  derivative  method  is  scan  shape  oriented,  and  works  best 
when  the  four  windows  just  fit  into  half  of  the  scan  return.  Shrinking 
the  size  of  the  windows  in  angle  makes  it  more  difficult  for  the  algo¬ 
rithm  to  detect  the  maximum  slope,  and  easier  to  trigger  on  noise. 
Spreading  out  the  windows  forces  any  major  change  in  signal  strength  to 
the  target  related,  not  noise  related.  Therefore,  the  coarser  the  gran¬ 
ularity  the  better  the  algorithm  works  for  a  given  number  of  beam  loca¬ 
tions  in  each  window.  The  most  favorable  granularity  is  the  antenna 
null-to-null  beamwidth  divided  by  the  number  of  beam  pointing  locations 
minus  one. 

The  flights  with  turbulence,  Figures  5-213  through  5-220,  again 
show  little  change  from  the  baseline  flights. 

The  flights  with  turbulence  and  noise.  Figures  5-221  through  5-228, 
again  show  the  second  derivative  to  be  the  best  estimator,  especially 
with  the  finer  granularities.  The  thresholding  methods  are  still  close 
to  each  other,  but  the  12  dB  method  appears  best,  followed  by  the  mean. 
RADARCG  is  a  bit  unusual,  as  the  mean  error  in  both  milliradians  and 
meters.  Figures  5-221  and  5-225,  show  it  to  be  quite  good,  yet  the  devia¬ 
tion  and  variance  of  the  estimator  is  quite  poor,  especially  in  actual 
meter  error,  Figures  5-227  and  5-228.  The  difference  between  the  mean 
error  in  meters.  Figure  5-225,  and  the  mean-square  error,  Figure  5-226, 
is  quite  large,  indicating  that  when  RADARCG  missed  the  target,  it  missed 
by  a  sizeable  amount. 
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Figure  5-213.  Mean  error  of  estimators  In  milllradians  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-214.  Mean  square  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-215.  Standard  deviation  of  error  of  estimators  in  milliradians 
with  the  antenna  located  125  M  from  the  runway  centerline. 
Each  data  point  is  the  result  of  one  flight,  all  scans 
used. 
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Figure  5-216.  Variance  of  error  of  estimators  in  mill iradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 


284 


00*91 


o/T2r3I3W  J#yy3  Ny3w  g 

^PP..V..  ,  00,9  00*h  t  00  vg 

1 8 


00 *2!  00*8  on  ’ 

SU313W  V0UU3  NH3M 


I  <✓> 
<T} 

o 

O 

r—  c 


•xs  o 

C  CL 

“ 

c 

o 

CD  rQ 

o 

*W  +-> 
C  * 

.co 

*TS  -O 

.a* 

0)  -C 

x;  u 

■*->  (t3 

,c--g 

■w  V) 

So 

„  =  </> 
</»  -p-  c 
L.  i—  ig 

<U  S-  u 
+->  0)  in 

SJ  4J 
E  C  r- 
(V  r— 
c  U  IO 


£ 

3  §.? 
I  ^ 

*r— 

*>r  o 
>p  *->  c 

E  ° 

<♦-0  4- 
O  S-  O 
H- 

S-  4J 

2£  = 
t-  W  W 
<U  cm  qj 

C  *" 

lO  t5  <u 
'■  <U  JC 


00  "u0 


£ 


Csi 

I 

IT) 

0> 

3 

cn 


I 


.1 


:i 

3 

I 

I 


BASELINE  FLIGHT  WITH  TURBULENCE 
O  MEAN  -I-  RADARCG  <!>  SDRV 

A  MEDIAN  X  1 2DB 


285 


Figure  5-218.  Mean  square  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-219.  Standard  deviation  of  error  of  estimators  in  meters  with 

the  antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-220.  Variance  of  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-221.  Mean  error  of  estimators  in  mill iradians  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-222.  Mean  square  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-223.  Standard  deviation  of  error  of  estimators  in  milliradians 
with  the  antenna  located  125  M  from  the  runway  centerline. 
Each  data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-224.  Variance  of  error  of  estimators  in  mil  1 iradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-225.  Mean  error  of  estimators  in  meters  with  the  antenna  located 
125  M  from  the  runway  centerline.  Each  data  point  is  the 
result  of  one  flight,  all  scans  used. 
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Figure  5-226.  Mean  square  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 


Figure  5-227.  Standard  deviation  of  error  of  estimators  in  meters  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each  data 
point  is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-228.  Variance  of  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  used. 
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Figure  5-229.  Mean  error  of  estimators  in  milliradians  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB 
or  greater  used. 
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Figure  5-230.  Mean  square  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  13  dB  or  greater  used. 
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Figure  5-232.  Variance  of  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  13  dB  or  greater  used. 
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Figure  5-233.  Mean  error  of  estimators  In  meters  with  the  antenna  loca¬ 
ted  125  M  from  the  runway  centerline.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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gure  5-234.  Mean  square  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-235.  Standard  deviation  of  error  of  estimators  in  meters  with 

the  antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  13  dB  or  greater  used. 
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Figure  5-236.  Variance  of  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 


304 


Proceeding  to  the  scans  with  a  SNR  of  at  least  13  dB,  Figures  5-229 
through  5-236,  the  flight  with  noise  only  shows  the  RADARCG  estimator  to 
be  the  most  accurate  and  reliable  with  respect  to  granularity.  The  sec¬ 
ond  derivative  is  generally  better  than  the  thresholding  methods,  deci¬ 
sively  so  in  the  finer  granularities.  The  thresholding  methods  are  of 
equal  quality  in  high  SNRs. 

The  flights  with  turbulence  and  noise,  Figures  5-237  through  5-224, 
again  show  RADARCG  to  be  the  best  estimator  in  periods  of  high  signal 
strength.  The  second  derivative  is  the  next  best  method,  degrading  to 
the  level  of  the  thresholding  estimators  for  the  granulari ties  between 
15  and  21  beam  locations  in  the  scan.  The  mean  estimator  appears  the 
best  in  accuracy  of  the  thresholding  methods,  followed  by  12  dB  and  the 
median  estimators,  in  that  order. 

The  scans  with  SNR's  of  10  dB  or  less  for  the  flights  with  noise 
are  plotted  in  Figures  5-245  through  5-252.  The  degradation  of  all  esti¬ 
mators  is  evident.  RADARCG  and  SDRV  are  the  two  best  estimators  particu¬ 
larly  in  the  higher  granularities.  The  12  dB  thresholding  method  again 
distinguishes  itself  as  the  best  of  the  thresholding  methods  in  noise, 
although  the  mean  and  median  methods  are  similar  in  ability. 

The  flights  with  turbulence  and  noise,  Figures  5-253  through  5-260, 
now  show  SDRV  to  be  the  most  accurate,  and  RADARCG  to  be  the  least  accu¬ 
rate.  The  thresholding  methods  do  not  appear  to  be  affected  by  the  tur¬ 
bulence.  It  is  shown  that  there  is  no  general  improvement  in  the  ac¬ 
curacy  of  the  thresholding  estimators  with  an  increasing  number  of  beam 
pointing  locations  in  the  scan.  If  the  target  returns  are  below  the 
noise  floor,  then  no  amount  of  data  will  display  the  target  without  some 
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Figure  5-237.  Mean  error  of  estimators  in  milliradians  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-238.  Mean  square  error  of  estimators  in  milliradians  with  the 
antenna  located  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  dB  or  greater  used. 
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Figure  5-239.  Standard  deviation  of  error  of  estimators  in  mill iradians 
with  the  antenna  located  125  M  from  the  runway  centerline. 
Each  data  point  is  the  result  of  one  flight,  all  scans 
with  a  SNR  13  dB  or  greater  used. 
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Figure  5-241.  Mean  error  of  estimators  in  meters  with  the  antenna  loca¬ 
ted  125  M  from  the  runway  centerline.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-242.  Mean  square  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-243.  Standard  deviation  of  error  of  estimators  in  meters  with 
the  antenna  located  125  M  from  the  runway  centerline. 
Each  data  point  Is  the  result  of  one  flight,  all  scans 
with  a  SNR  13  dB  or  greater  used. 


Figure  5-244.  Variance  of  error  of  estimators  in  meters  with  the  antenna 
located  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  5-246.  Mean  square  error  of  estimators  in  milliradians  with  the 
antenna  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  at  or 
below  10  dB  used. 
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Figure  5-247.  Standard  deviation  of  error  of  estimators  in  mili'iradians 
with  the  antenna  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  at  or  below  10  dB  used. 
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Figure  5-248.  Variance  of  error  of  estimators  in  milliradians  with  the 
antenna  125  M  from  the  runway  centerline.  Each  data  point 
is  the  result  of  one  flight,  all  scans  with  a  SNR  at  or 
below  10  dB  used . 
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Figure  5-249.  Mean  error  of  estimators  in  meters  with  the  antenna  125  M 
from  the  runway  centerline.  Each  data  point  is  the  result 
of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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Figure  5-250.  Mean  square  error  of  estimators  in  meters  with  the  antenna 
125  M  from  the  runway  centerline.  Each  data  point  is  the 
result  of  one  flight,  all  scans  with  a  SNR  at  or  below  10 
dB  used. 


FLIGHT  WITH  NOISE 

O  MEAN  +  RADARCG  <!>  SDRV 

a  MEDIAN  X  12DB 


329 


w  'yoyy3  yo  noiihia30  ayuaNyisg 

00*91  00'hf  00  *0  1  00  *9  00  *S* 


Figure  5-251.  Standard  deviation  of  error  of  estimators  in  meters  with 
the  antenna  125  M  from  the  runway  centerline.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  at 
or  below  10  dB  used. 
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Figure  5-252.  Variance  of  error  of  estimators  in  meters  with  the  antenna 
125  M  from  the  runway  centerline.  Each  data  point  is  the 
result  of  one  flight,  all  scans  with  a  SNR  at  or  below  10 
dB  used. 
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Figure  5-254.  Mean  square  error  of  estimators  in  milliradians  with  the 
antenna  125  M  from  the  runway  centerline.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR 
at  or  below  10  dB  used. 
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sort  of  filtering  action.  The  thresholding  methods  have  none,  so  their 
accuracy  is  limited  by  the  amount  of  signal  rising  out  of  the  noise 
floor.  With  enough  beam  locations  in  each  window,  SDRV  can,  in  theory, 
be  able  to  average  out  the  noise  in  the  target.  When  averaging  is  not 
done,  as  in  the  granularities  at  and  below  21  beam  locations,  SDRV  is 
only  marginally  better  than  the  thresholding  methods, as  shown  by  the 
standard  deviation  and  variance  in  mil  1 iradians ,  Figures  5-255  and  5-256. 
When  averaging  occurs  in  the  higher  granularities,  the  method  stands  by 
itself. 
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Figure  5-255.  Standard  deviation  of  error  of  estimators  in  milliradian 
with  the  antenna  125  M  from  the  runway  centerline.  Each 
data  point  is  the  result  of  one  flight,  all  scans  with  a 
SNR  at  or  below  10  dB  used. 
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Figure  5-256.  Variance  of  error  of  estimators  in  milliradians  with  the 
antenna  125  tl  from  the  runway  centerline.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR  at 
or  below  10  dB  used. 
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Figure  5-257.  Mean  square  of  estimators  in  meters  with  the  antenna  125  M 
from  the  runway  centerline.  Each  data  point  is  the  result 
of  one  flight,  all  scans  with  a  SNR  at  or  below  10  dB 
used. 
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Figure  5-258.  Mean  square  error  of  estimators  in  meters  with  the  antenna 
125  M  from  the  runway  centerline.  Each  data  point  is  the 
result  of  one  flight,  all  scans  with  a  SNR  at  or  below  10 
dB  used. 


Figure  5-259.  Standard  deviation  of  error  of  estimators  in  meters  with 
the  antenna  125  M  from  the  runway  centerline.  Each  data 
point  is  the  result  of  one  flight,  all  scans  with  a  SNR 
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Figure  5-260.  Variance  of  error  of  estimators  in  meters  with  the  antenna 
125  M  from  the  runway  centerline.  Each  data  point  is  the 
result  of  one  flight,  all  scans  with  a  SNR  at  or  below  10 
dB  used. 
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VI,  SUMMARY 

The  subject  of  this  work  was  to  introduce  a  new  non-thresholding 
algorithm,  and  to  determine  the  estimating  ability  of  it  with  a  second 
non-thresholding  algorithm  (radar  center  of  gravity)  and  three  thresh¬ 
olding  techniques,  two  post-determined  (mean  and  median)  and  one  pre¬ 
determined  (12  dB).  It  was  observed  that  during  periods  of  high  signal  - 
to-noise  ratio  the  radar  center  of  gravity  was  superior  to  the  other 
methods,  with  the  second  derivative  method  a  very  close  second.  With  the 
same  signal-to-noise  ratios,  the  thresholding  methods  were  of  equal  es¬ 
timating  ability,  with  a  slight  favoring  of  the  12  dB  threshold.  The 
addition  of  turbulence  alone  into  the  simulation  did  not  appreciably  af¬ 
fect  the  error  produced  by  the  estimators  since  there  was  no  noise  floor 
present  in  the  computer  simulation.  Therefore,  the  signal-to-noise  ra¬ 
tio  was  infinite,  regardless  of  the  actual  return  amplitude  of  the  scan, 
and  the  flights  with  only  turbulence  were  really  just  baseline  flights 
revisited. 

When  15  dB  of  zero-mean  random  gaussian  noise  was  added  to  the 
baseline  flight  or  flight  with  turbulence  scan  returns,  the  second  deri¬ 
vative  method  emerged  as  the  most  accurate  estimator,  especially  in  pe¬ 
riods  of  low  signal-to-noise  ratio.  The  radar  center  of  gravity  is  de¬ 
graded  because  it  always  makes  a  location  estimate,  regardless  of  whether 

or  not  a  target  is  present.  With  enough  beam  locations,  the  noise  will 
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average  out,  and  the  estimator  will  improve  as  shown  by  the  plots  of 
error  versus  granularity.  The  thresholding  methods  also  make  an  estimate 
of  target  location  regardless  of  the  presence  or  absence  of  a  target, 
except  in  the  rare  instance  when  no  two  consecutive  returns  are  above 
the  scan  threshold.  In  this  simulation,  if  there  is  no  target,  but 
only  noise,  the  mean  estimator  will  calculate  a  threshold  close  to  zero. 
Since  gaussian  noise  is  randomly  positive  and  negative,  the  probability 
that  all  returns  in  the  scan  be  alternately  positive  and  negative  is  low. 
Therefore,  the  probability  that  the  there  is  at  least  one  instance  where 
two  consecutive  returns  are  positive  is  high,  and  a  "target"  is  found, 
since  two  positive  returns  are  above  the  threshold.  Again  in  the  noise 
only  case,  the  median  estimator  will  set  the  threshold  at  the  median 
return  amplitude,  which  will  be  approximately  zero,  and  the  above  argu¬ 
ment  holds  for  this  estimator  also.  In  the  case  of  the  12  dB  estimator, 
12  dB  down  from  a  noise  voltage  places  the  threshold  in  noise,  and  the 
above  argument  can  be  reapplied.  Therefore,  none  of  the  above  estimators 
are  truely  capable  of  target  detection,  since  the  estimators  will  do 
their  best  to  find  a  "target"  regardless  of  whether  or  not  one  is  present 
Only  the  second  derivative  method  is  able  to  make  a  decision  as  to  the 
absence  or  presence  of  a  target  by  operating  on  the  shape  of  the  scan 
return.  If  no  target  is  present,  the  shape  of  the  return  is  flat  with 
zero  slope  and  a  constant  second  derivative.  If  a  target  is  present, 
the  slope  will  change  and  the  second  derivatives  of  each  set  of  estimator 
windows  will  change  sign,  locating  the  point  of  maximum  slope  and  there¬ 
fore  a  target  edge.  The  plots  of  scan  SNR  versus  average  number  of  scans 
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Figures  4 - 1 3 ,  14,  17,  18,  21,  and  22,  show  the  second  derivative  method 
to  '-eject  those  scans  of  low  SNR.  The  plots  of  error  for  low  SNR,  in 
particular  Figures  5-253  through  5-260,  show  that  when  a  scan  is  accepted 
and  a  decision  of  target  location  is  made,  that  the  decision  is  accurate, 
because  even  in  noise  the  basic  scan  shape  was  suitable  enough  to  detect 
and  locate  a  target. 

The  figures  of  error  versus  antenna  offset  show  no  apparent  in¬ 
crease  in  error  due  to  antenna  offset.  This  is  primarily  due  to  the 
high  PRF  which  "stops"  the  target  in  the  flight  path  during  the  scan. 

Even  with  a  large  granularity,  the  target  will  move  only  on  the  order  of 
its  own  length  during  the  scan.  So  the  effect  due  to  translation  is  small. 
Since  the  target  movement  is  small,  in  translation  as  well  as  in  rotation 
due  to  turbulence,  the  modulation  of  the  radar  cross  section  of  the  tar¬ 
get  is  small,  and  that  effect  is  also  negligible.  Therefore,  the  results 
of  this  computer  simulation  is  accurate  in  those  regards.  Yet,  in  actual 
practice,  error  is  expected  when  the  antenna  is  delocated  from  the  runway. 
This  error  is  introduced  by  the  ranging  accuracy  of  the  radar.  If  the 
antenna  is  a  distance  to  the  left  of  the  runway,  and  the  range  estimate 
is  short  of  the  actual  target  location,  then  the  estimate  of  position 
will  place  the  target  to  the  left  of  the  runway  independent  of  the  error 
to  the  target  centroid  in  milliradians. 

The  plots  of  error  versus  granularity  display  little  decrease  in 
error  with  an  increase  in  granularity  for  the  thresholding  methods,  but 
dees  show  an  increase  in  robustness  with  granularity.  Both  non-thresh¬ 
olding  methods  improve  in  both  average  error  and  robustness  with  granula¬ 
rity. 
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Overall,  The  second  derivative  algorithm  is  clearly  the  estimator 
least  prone  to  error  and  most  robust  in  estimating  ability  for  virtually 
all  granularities  in  the  presence  of  noise.  The  radar  center  of  gravity 
estimator  is  most  accurate  in  periods  of  high  SNR  and  with  fine  granular¬ 
ities.  It  suffers  degradation  in  noise,  and  in  large  noise  becomes  unus¬ 
able  as  an  estimator  due  to  noise  outliers  in  the  scan  returns.  The 
thresholding  methods  are  rather  robust  but  less  accurate  than  the  second 
derivative  method  and  RADARCG  with  high  SNRs. 

Pulse  economy  is  of  great  importance.  That  estimator  which  uses 
the  fewest  pulse  transmissions  without  an  increase  in  target  location  will 
in  general  be  the  estimator  employed  in  practice.  The  12  dB  thresholding 
method  is  especially  suited  to  pulse  economy.  When  a  target  edge  is  found, 
scanning  translates  to  the  other  side  of  the  scanning  window,  savina  Dulse 
transmissions  by  removing  the  requirement  to  fully  scan  the  target.  The 
second  derivative  method  does  require  that  the  target  be  scanned  more 
fully  than  the  12  dB  method,  as  shown  in  Figures  5-1  to  5-4.  While  this 
reduces  the  probability  of  false  alarm,  it  requires  the  expenditure  of 
additional  se  transmissions  in  the  scanning  window.  However,  virtually 
every  data  point  on  the  256  figures  in  the  previous  chapter  have  shown  the 
mean  error,  mean  square  error,  and  standard  deviation  and  variance  of  the 
error  of  SDRV  to  be  less  than  the  12  dB  estimator.  How  do  these  methods 
compare  in  terms  of  target  location  accuracy  versus  pulse  budget? 

A  plot  of  average  pulses  used  per  scan  v  rsus  granularity  is  shown 
in  Figure  6-1.  It  should  be  pointed  out  that  in  this  simulation,  time  was 
incremented  as  though  3  pulses  were  transmitted  at  each  beam  location, 
but  only  one  pulse  was  used.  Since  there  is  a  one  to  one  correspondence 
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Figure  6-1.  Average  number  of  pulses  used  per  scan.  Each  data  point  is 
the  result  of  103  scans  on  the  target. 
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between  pulses  per  scan  and  the  number  of  beam  pointing  locations  in  the 
scan,  the  mean,  median  and  RADARCG  estimators  use  a  number  of  pulses 
equal  to  the  scanning  granularity.  The  12  dB  method  transmits  one  pulse 
in  the  direction  of  the  target  to  set  the  threshold,  and  thus  uses  one 
pulse  before  the  start  of  the  scan.  The  second  derivative  estimator  has 
the  capability  of  rescanning  part  of  the  target  twice  due  to  the  nature 
of  the  algorithm.  As  such,  it  can  use  more  pulses  than  the  magnitude  of 
the  scanning  granularity.  Referring  tc  Figure  6-1,  the  12  dB  estimator 
is  shown  to  use  about  half  of  the  pulses  budgeted  for  the  scan.  The 
second  derivative  stays  very  close  to  the  pulse  budget,  sometimes  re¬ 
scanning  the  target  more  frequently  at  certain  granularities.  Those  are 
the  granularities  where  the  number  of  beam  locations  in  each  window  of 
the  estimator  increases  by  one.  The  highest  granularities  that  do  not 
increase  the  number  of  beam  locations  in  the  windows,  namely  the  gran¬ 
ularities  21,  29,  and  45,  have  the  least  pulse  usage  but  the  highest 
error. 

In  periods  of  high  SNR,  what  should  be  expected  with  regards  to 
pulse  usage?  Figure  6-2  is  a  plot  of  the  ayerage  number  of  pulses  used 
per  scan,  only  those  scans  of  at  least  13  dB  selected.  There  is  almost 
no  change  in  the  number  of  pulses  used  by  SDRV.  The  12  dB  method  used 
more  pulses,  because  it  had  set  the  threshold  higher  due  to  the  increase 
in  return  amplitude. 

In  periods  of  low  SNR,  as  shown  by  Figure  6-3,  the  12  dB  method 
used  fewer  pulses,  because  '.t  was  able  to  trigger  on  noise  at  the  ends  of 
the  scanning  window,  rather  than  triggering  on  the  target  itself.  Also 
as  expected,  there  was  a  slight  increase  in  the  number  of  pulsed  used  by 


336 


NdOS  d3d  Q3SH  S3S~ind  30dU3AB 


Figure  6-2.  Average  number  of  pulses  used  per  scan.  Each  data  point  is 
the  result  of  one  flight,  all  scans  with  a  SNR  13  dB  or 
greater  used. 
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Figure  6-3.  Average  number  of  pulses  used  per  scan.  Each  data  point  is 
the  result  of  one  flight,  scan  SNR  at  or  below  10  dB  used. 
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SDRV  as  it  was  forced  to  search  further  in  the  scan  to  locate  a  target, 
the  noise  alone  being  unable  to  satisfy  the  algorithm. 

But  we  have  not  yet  addressed  the  posed  question  -  what  is  the  re¬ 
lationship  between  pulses  used  and  error  induced?  Returning  to  Figure 
5-221,  let  us  choose  the  granularity  where  the  12  dB  estimator  is  the  most 
accurate.  This  would  be  a  granularity  of  45  beam  locations  in  the  scan¬ 
ning  window.  From  Figure  6-1,  the  12  dB  method  used  approximately  25 
pulses  per  scan,  on  the  average,  at  that  granularity.  Moving  directly 
across  the  figure  to  the  left,  it  is  observed  that  the  second  derivative 
method  used  about  25  pulses  per  scan  with  a  granularity  of  27  beam  loca¬ 
tions  in  the  scanning  window.  Returning  to  Figure  5-221,  the  mean  error 
of  the  second  derivative  algorithm  at  a  granularity  of  27  beam  locations 
is  approximately  1  mrad.  From  the  same  figure,  using  27  pulses  per  scan 
(granularity  of  45),  the  12  dB  estimator  had  a  mean  error  of  about  2.7 
mrad.  Further,  comparison  of  the  results  of  SDRV  at  a  granularity  of  27, 
Figures  5-197  through  5-260,  with  the  results  of  the  12  dB  estimator  at  a 
granularity  of  47,  shows  the  second  derivative  algorithm  to  be  clearly 
superior.  Since  the  mean,  median,  and  RADARCG  estimators  all  use  more 
pulses  than  the  12  dB  threshold,  SDRV  is  also  more  efficient  per  pulse  than 
those  methods. 

It  is  instructive  to  examine  the  error  of  the  estimators  as  a  func¬ 
tion  of  range.  Since  the  12  dB  estimator  is  representative  of  the  thresh¬ 
olding  estimators,  only  the  data  for  it  was  plotted  against  the  output  of 
the  RADARCG  and  SDRV  methods  in  the  accompanying  figures.  The  data  was 
obtained  by  using  the  results  of  the  eleven  flights  with  noise  and  turbu¬ 
lence  representing  each  antenna  offset  location  at  all  21  scanning  granu- 
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larities,  regardless  of  scan  SNR.  The  data  is  plotted  in  Figures  6-4  to 
6-11,  error  as  a  function  of  scan  number.  The  first  scan  occurred  with 
the  target  4.0  nmi  from  the  antenna,  and  the  103rd  scan  occurred  with  the 
target  over  the  release  point.  Clearly  shown  is  the  consistency  of  the 
12  dB  thresholding  method  and  the  SDRV.  method  with  about  half  the  error. 
Also  clearly  shown  is  the  marked  lack  of  robustness  of  RADARCG.  This  is 
most  pronounced  at  far  range,  and  decreases  to  the  level  of  SDRV  at  close 
range.  Note  that  RADARCG  appears  to  be  more  accurate  but  less  robust  than 
SDRV  at  close  range  from  the  plots  in  milliradians,  Figures  6-4  to  6-7, 
but  was  actually  approaching  the  quality  of  SDRV  as  seen  in  the  plots  with 
error  in  meters,  Figures  6-8  to  6-11.  Finally,  the  plots  show  the  error 
to  reduce  as  the  target  approaches,  an  expected  result. 

It  is  appropriate  to  conclude  this  work  with  a  comparison  of  the 
estimators  on  a  pulse-by-pulse  basis.  Since  a  typical  radar  system  has  a 
limited  pulse  budget,  it  is  prudent  to  employ  the  centroid  estimator  which 
uses  the  pulses  expended  most  efficiently.  The  following  figures  are 
plots  of  the  average  (mean)  error  of  an  estimator  per  pulse  as  a  function 
of  the  average  number  of  pulses  used  in  the  scan.  Each  data  point  is  the 
result  of  11  flights  of  which  the  total  error  was  divided  by  the  total 
number  of  pulses  used  to  obtain  a  single  data  point.  Figures  6-12  and 
6-13,  the  error-pulse  ratio  in  milliradians  and  meters,  respectively,  used 
all  scans  regardless  of  SNR.  The  second  derivative  estimator  is  clearly 
shown  to  use  the  available  pulse  budget  most  efficiently  by  having  the 
least  error  per  pulse  for  any  amount  of  pulses  budgetted  in  that  scan. 

Note  that  the  SDRV  method  is  most  efficient  with  pulse  integration,  and 
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Figure  6-5.  Mean  square  error  in  milliradians  as  a  function  of  scan  number 
for  flights  with  noise  and  turbulence,  all  scans  used. 
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Figure  6-6.  Standard  deviation  of  error  in  mi  1 1 i radians  as  a  function  of 
scan  number  for  flights  with  noise  and  turbulence,  all  scans 
used. 


Figure  6-7.  Variance  of  error  in  mi  11 i radians  as  a  function  of  scan  number 
for  flights  with  noise  ami  turbulence,  all  scans  used. 
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Figure  6-8.  Mean  error  in  meters  as  a  function  of  scan  number  for  flights 
with  noise  and  turbulence,  all  scans  used. 
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Figure  6-9.  Mean  square  error  in  meters  as  a  function  of  scan  number  for 
flights  with  noise  and  turbulence,  all  scans  used. 


Standard  deviation  of  error  in  meters  as  a  function  of  scan 
number  for  flights  with  noise  and  turbulence,  all  scans  used 
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Variance  of  error  in  meters  as  a  function  of  scan  number  for 
flights  with  noise  and  turbulence,  all  scans  used. 
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Linear  error-pulse  ratio  for  flights  with  noise  and  turbulence, 
all  scans  used. 
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that  the  other  estimators  approach  the  quality  of  SDRV  at  low  pulse  usage 
only  with  finer  granularities. 

For  comparison,  in  periods  of  high  SNR,  RADARCG  is  of  the  same 
quality  as  SDRV  with  pulse  integration,  as  shown  in  Figures  0-14  and  6-15. 
The  thresholding  methods  are  of  equal  quality.  As  is  to  be  expected,  SDRV 
is  the  most  efficient  estimator  in  low  SNR,  Figures  6-16  and  6-17,  with 
RADARCG  least  efficient.  Here  the  effect  of  pulse  integration  in  the  SDRV 
algorithm  is  most  dramatic. 

In  conclusion,  the  second  derivative  algorithm  is  shown  to  have  the 
least  error  per  pulse  used  in  comparison  to  the  other  estimators  employed 
in  this  simulation.  It  is  the  most  robust  of  those  compared  due  to  the 
unique  scan  rejection  capability  inherent  in  the  architecture  of  the  algo¬ 
rithm. 


Figure  6-14.  Angular  error-pulse  ratio  for  flights  with  noise  and  turbulence 
all  scans  with  a  SNR  13  dB  or  greater  used. 
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Ingular  error-pulse  ratio  for  flights  with  noise  and  turbulence, 
ill  scans  with  a  SNR  at  or  below  10  dB  used. 
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Figure  6-17.  Linear  error-pulse  ratio  for  flights  with  noise  and  turbulence, 
all  scans  with  a  SNR  at  or  below  10  dB  used. 
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This  appendix  contains  a  complete  listing  of  the  FORTRAN  program 
used  to  compute  the  flight  of  the  target  and  the  output  of  the  centroid 
estimators.  The  program  is  extensively  documented  and  a  reading  of  the 
comment  cards  reveals  much  about  the  computations  that  are  performed. 

The  program  requires  15  minutes  to  compile  and  execute  one  flight 
with  a  granularity  of  29  beam  pointing  locations  on  an  IBM  3031  computer 

The  plots  presented  in  the  text  are  special  purpose  programs  devel 
oped  to  be  compatible  with  the  plotting  facilities  at  Auburn  University. 


C  NOTE:  THIS  IS  INS  3ATCALS  5  AD  Ac  3I3CLATI0N  AS  HOD'FIED  A?.-.  193  1 
C 

c 

c 

c 

c«. . . 

c 


C  PA3A.1ETE3  LIST 

C 

C  A 

C  Ag-D'JSHY  7A3IA3L2,  COSTAINS  THE  S’JN  0?  3A.NDC3  FLOATING  POINT  SOSoisS 
C  A  Z  =  3  A I  N  P30G:  ANTENNA  A l IN  DTK  TO  EACH  SC  ATT  E  EE  3 

C  S  03  HOOT  I  N  ES  :  3  AW  3ETU3MS  IN  A2I30T.H 

C  A  2  D 1 3  =  All  NOTH  3INENSI0.N  ON  THE  SCANNING  VINDD7. 

C  A  23  =  All  NOTH  ANGIE  TO  I A3 GET  AXIS 

C  A  2  5  W  =  HAL?  0?  THE  A2I NOTH  32A3WIDTH  EETWEEN  NULLS. 

C  A2S?C=  All  NOTH  SPACING  0?  3EAN3  IN  32A3WI0THS,  3  E  3  TO  I  N  C  ?  E  N  E  S  T 

C  THE  ANTENNA  3E AN  POINTING  POSITION  (IE  TOO  WANT  THE  3EAN  TO 

C  NOTE  TWC  POSITIONS  AT  A  SHOT,  SET  A2S?C  =  2). 

C  A 2 T H S S  =  AZINOTH  TH3ESKOLO  D-T23NINE0  E30N  SEAN,  NEMAN,  0?  TWELVE 

C  03  S03300TENE3. 

C 

C  B 

C 

c  c 

C  C3LCCA,C3LOCE=>  CENTS?.  3  E  AN  LOCATION  IN  A  E  I  NOTH ,  E  LE  7  ATI  0  N  IN  ?  A  DA  S 
C  COO 3 01  NATES.  EJOAL  TO  PkE.PEL  WHEN  L0C:<=1. 

C  C3I*  YAW  ANGLE, 

c 

c  0 

C  DG  A I  N  =  ANTENNA  GAIN 

C  DELTA  *  3  A  X I N  U  3  ANTENNA  POINTING  E3E03. 

C 


C  E35CT*  3  ACES  CAT  T  E3  E  0  E-FIZLO. 

C  E  L  *  N  A  I N  P30G :  ANTENNA  EL  S7  AT  TON  TO  EACH  SCAT  TE  3  E  E 

C  TH3 E S HOLDING  3U3H0UTINSS:  SAW  PETOr.NS  EN  EL37 AEION 

C  EL3W=HAL?  0?  AN  ELEVATION  3VFN. 

C  ELD  1 3*  ELE7ATI0N  OI3ENSION  ON  THE  SCANNING  WINDOW. 

c  elo=elevai:om  io  taeget  axis 

C  ELSCA  N*ON/Of ?  INPUT  V  A3  I  A3  LE  C0NT30L  ?03  ELEVATION  SCAN  OPTION. 

C  ELSPC*  ELE  7  A  TI ON  SPACING  OF  3EA33  IN  3 EA 3«I OTH S . 

C  ELT3  ?S  =  ELE7ATI0N  TH32SHOLO  0ETE3  3  IN  EC  3T  1 E  A  N ,  3E0IAN,  C?  TWELVE 
C  OB  S U3 30 U TINES. 

C 

C  F 

C  FH  LH=  FI3ST  HIT,  LAST  HIT  THE23  HOLDING  PF.OCEOUEE 

C  FLTI3E*  S I  30  L ATI  ON  TI3E  IN  INTiSSi  SECONDS. 

C  P?.E2=F?E0O  2NCY. 

c  fxa*»:nal  antenna  x  cooscinate  location 

C 

C  G 

C  G  A  3  3  A  *  PITCH  ANGLE. 

C  3LSL?= SLIDESLOPE  IN  DEG5EES. 
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c  grana*  azisuth  null-to-null  scan  granularity,  detsshined  3y  addin. 

C  GB ANE*  AS  IN  GHANA,  3UT  IN  EL27ATIDN. 

C  GRNG* PROJECT  ION  OF  ?.NG  ON  X“  T  PLANS. 

C 

c  a 

C  HEIGHT' HEIGHT  OF  ANTENNA  CENT'S. 

C 

C  I 

C  I CODE'CODS  ?C H  CHOICE  OF  SCAT? ES-TN G  SODELS  IN  SCS. 

C  IMCXA* INTEGER  3Y  WHICH  ANTENNA  X  COORDINATE  IS  I NC REN 2N TED . 

C  I Y = R  A  N  D  C  N  NON  3  S3  GSNSEAT03  ENT3ANCE  NU33ER. 

C  IXA'INTI AL  ANTENNA  X  COORDINATE. 

C  IT'INTEGES  OUTPUT  0?  RANDOM  NUMBER  GESERATOS  AND  USED  AS  THE 
C  NEW  INPUT  10  THE  3  AS  D  0  BOOT  IN  E  THE  NEXT  TINE  IT  IS  CALLED. 

C 

C  J 

C 

C  K 

C  X*  'WATS  NUS3SR. 

C 

C  L 

C  LABEL'  A  DON. NT  7A3IA3LS  USED  TO  KEEP  T2ACK  OF  THE  FLOW  OF 
C  DATA  IN  THE  CSNTSOID  CALCULATIONS 

C  LAN3DA=WA7ELEHGTH. 

C 

c  a 

C  HE  A  NA=  SEAN  AZI3UTH  70LTAGE 

C  N  8  A  N  S*  NEAN  ELE  7  ATI  ON  70LTAGE 

C  SEDA, NEDS'  SEDIAN  70LTAGE  RETU3N  IN  AZIMUTH  AND  ELS7ATI0N 
C  DETE3SINED  BT  SEDIAN  SC3S0UTI NS  WHICH  IS  PASSED 

C  TO  FHLH  AS  THB3SH0LD. 

C  HTI»N07ING  TAaGET  INDICATION. 

C  HU-  30LL  ANGLE. 

C 

C  N 

C  NL0C*NUS3ES  OF  T3AN SS I S S IONS  IN  EACH  SCANNING  WINDOW. 

C  NOH*  NUN3ES  OF  SCA7T3SS3S  IN  EACH  TARGET  COSPLEX. 

C 

C  0 

C  0HEGA'3ADIAS  PRZQUSNCY  OF  THE  TA3G2T  NOTIONS,  IN  3AD/SSC0ND. 

C 

C  P 

C  PAHGLZ*  PRESENT  ANGLE  TO  CROSS-TRACK  CENTE3  AS  SEEN  P303  RUNWAY. 

C  PC A*  PHI  10  TARGET  CENTROID  IN  AZISUTH.  DETER HINED  9T  THRESHOLDING 
C  PROCEDtJ  3SS  AND  PASSED  TO  ERROR  SU3HOOTINE  FOR  COSPARISON 

C  AGAINST  ACTUAL  TARGET  LOCATION. 

C  PERIOD*  PERIOD  OF  THE  SELECTED  TARGET  NOTION,  IN  SECONDS. 

C  PSI'SPHERICAL  COOaDINATE.  (TARGET) 

C  PR  ATE*  PITCH  RATE  OF  THE  TARGET  SO  DSL,  IN  DEG RSES/SEC3 N D . 

C  PWRT-PEAK  TRANSaiT  POWER. 

C 

C  Q 

C 

C  R 

C  RCGAZ*  RADAR  CENTER  0?  G8A7TTT  IN  AZI30TH  AS  SEEN  ?T 
C  ROHWAT  AND  DETEESINED  ST  RADARC. 
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c  pcgel*  radar  csntsh  op  gravit*  in  elevation  is  seen  av 

C  RUNWAT  ISO  DETERMINED  31  RADASC. 

C  RECVOL-INPUT  OPTION  TO  DISPLIT  SCAM  RETURN  VOLTAGES. 

C  SNG*  BANGS  TO  EACH  SCATTE3EB  IN  TBS  TARGET  COMPLEX. 

C  RMGTGT-RANGE  TO  THE  T1RGET ,  IS  NAUTICAL  BILES. 

C  S3 ATE*8 01 L  RATE  OF  THE  TARGET  33 DEL ,  IN  DEGR EES /SECOND . 

C 

C  S  ‘ 

C  SCAT*  AN  A  SB  AT  DIMENSIONED  (ISCAT,»).  CONTAINS  COORDINATES  OF 
C  EACH  SCATTEBER  IN  TARGE?  COORDINATES  AND  THE  RCS  OF  EACH. 

C  SIGNLA.SIGNLE-  THE  SIGNAL  AHPLITODE  NHEN  THE  3EA3  IS  CENTERED  ON 
C  TARGET  IN  AZIMUTH  AND  ELEVATION,  USED  TO  DSTER3INE 

C  THE  SNR  FOR  THAT  SCAN. 

C  SNRAZ.3N3ZL*  THE  SIGNAL  TO  NOISE  (POWER)  RATIO  DETE33I NED  3T 
C  (SIGNAL  AHPLITUDE/SHS  NOISE  VOLTAGE) *SQRT (2*SNR) 

C  SPEED*  TARGET  SPEED  IN  MPH. 

C  SONS*  ARGUHENT  USED  TO  SOB  VOLTAGES  TOGETHER  IN  3ADARC. 

C  SUHSTH*  ARGUMENT  USED  TO  SUM  VOLTAGE  RETURN  TIMES  THE  ANGLE 
C  TO  THAT  RETURN  AS  SEEN  PROS  THE  ANTENNA. 

C 

C  T 

C  TCS*  THETA  TO  TARGET  CENTSOID  IN  ELEVATION,  DETERMINED  BI 
C  THRESHOLDING  PROCEDURES  AND  PASSED  TO  ERROR  SUBROUTINE 

C  FOR  COMPARISON  AGAINST  ACTUAL  TARGET  LOCATION  AS  SEES 

C  8T  THE  ANTENNA. 

C  THETA* SPHERICAL  COORDINATE. (TARGET) 

C  THRESH-  THRESHOLD  USED  IN  FHLH  (FIRST  HIT  LAST  HIT)  SUBROUTINE 
C  TIMS*  TIME  VARIABLE  USED  FOB  DEAD  TIME  BETWEEN  SCANS. 

C  TIME!-  INITIAL  TIME  OF  SCENABIO  AND  TI3E  VABIA3L2. 

C  TI2EP-  PULSE  TISE 

C  TPHSFT*  PHASZ  SHIFTEa  RESPONSE  TIME. 

C 

c  u 

c  OCESTA.OCERT E»  THE  ASOUST  OP  33AH  POINTING  ERROR  INTRODUCED  INTO 
C  THE  PLACE3ENT  OP  THE  ANTENNA  SEAS,  CONSTANT  THHOOGH- 

C  OUT  A  SCAN,  BUT  R  ANDOMLT  CHANGING  FROM  SCAN  TO  SCAN. 

C 

C  V 

C  7DI8CT-  DIRECT  VOLTAGE  RETURN. 

C  VLT*  VOLTAGE  MATRIX  DIMENSIONED  (JSCAN,II).  JSCAN-1  ARE  AZIMUTH 
C  VOLTAGES,  JSCAN-2  ARE  ELEVATION  VOLTAGZS,  II-SCAN  SIZE. 

C  VNRBS*  THE  RMS  NOISE  VOLTAGE  FOR  THAT  SCAN  COMPUTED  BT  SQUARING 
C  EACH  BOISE  VOLTAGE,  SUAHING  THE B,  AND  THEN  TAKING  THE  SQUARE 

C  BOOT  OF  THE  AVERAGE  NOISE  70LTAGZ  SQUARED. 

C 

c  ■ 

c 

C  X 

C  IA,TA, HEIGHT*  INITIAL  COOIDINATES  OF  THE  SADAB  ANTENNA. 

C  XD,TD,!0*  DOUBLE  PHSCZSZON  RADAR  COORDINATES  OP  TABGET  LOCATION. 

C  XOIF,TDIF,ZDIP*  RADAR  COORDINATES  OF  THE  SCATTEEERS  WITH  RESPECT 
C  TO  THE  LOCATION  OF  THE  ANTENNA. 

C  TO , TO , ZO*  TARGET  POSTION  ALONG  TRAJECTORT,  IN  RADAR  COORDINATES. 

C 

C  T 

C  TPL-  SEAL  RANDOM  RUBBER  BETWEEN  0  AND  1  PROS  OUTPUT  DF  THE 
C  OOTPOT  OF  THE  aiNDU  SUBROUTINE. 


o  o  o  n  n  o 


TS  ATE*  TAN  3112  Of  THE  TARGET  MODEL,  IS  DEGasZS/SECOND 

Z 


C 


o 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DIMENSION  7LT(2,49)  ,  SCI?  (3,9)  ,TOD*T  (3)  ,AZB1DB(49)  ,VSOI  S2  (2,  49) 
DOOBLS  PRECISION  EL  7  ( 49)  ,  A  ZI  ( 4  9j 

I STEGZB  OPTION  (5)  , 3L SCAN , PXA ,R ETVOL.COO NT A ,COUSTE, PSP 
DOOBLE  PRECISION  X DI P, 1 0 IP , ZDI E , X D , TO, 20 
BEAL  BOR,  LAMBDA,  BO,  L7 ,  LS  ,  7  3!)S  (  2, 49 )  ,  X 

INTEGER  AZD1 1, SLDIS , DIM, AGILE , 3DI 7 , ADI  7 , f LTI BE , PLS 1 2 , PLSf D, PLSS D 
COSMOS/ BPT/SB  ATS, PS  AT E, TRATE ,0  SEGA , SPEED , GLSL P, 8MGT3T, 

•AZ.GRNG 

COMHON/RTA/MU , GAMMA, CHI, XO, TO, ZO 
COMHOB/SC ATM/SCAT, NOB 
COHMON/B  A7  E/K, MUR, EPS  B 

COSMON  /ALL/3 L7,AZI, RANGE, OPTION, LA  BEL, ELDIS.AZ DIM, 2LSCAN, DELT1 
COMMON  /SAC2/P2L,PAZ, C3LCCE, CBLOCA , GR A NE, GBAS A , A ZBB, EL3W , LOCK 
COMMON  /EILEIT/XSEAN, XM2D, XB ADRC , X 1 2DB , XSD RV , S NR AZ, PLS  12 . 
•PLSfO.PLSSD, DMEAN,DMED,D3ADBC, D12DB,DSDS7 
COM MON/POO  31 7/CO  ON  TE, COUNT A, MODES , MODE  A, MS  DR7 , SNREL, 30DE1 2, 30DA 12 

DATA  PI.PI2,T?I.3100SG,DEGRAD,S0P2,C/3. 141S92, 1.570796, 

•6.  29  31  35, 57.  2  9578,.  0  )7  4533,  12.  566  37,  3.  0E*Q8/ 

DATA  ITABLE,J0ST1/1 ,0/ 


DATA  ISPOT  FOBS  AT 
COL.  7  ABIA3L2  NAME  OPTIONS 

1  TCODS  1  -  JETSOD  »  1  AS  TARGET  MODEL 

2  -  ONE  POIST  ISOTROPIC  SCATTER EB 

2  NOISE  0  -  NO  GAUSSIAN  NOISE 

1  -  ADD  GABSSIAN  NOISE  TO  RETURN 

3  I HPI  0  -  NO  SOLI, PITCH,  OH  TAN 

1  -  ROLL,  PITCH  AND  TA*  ARE  INCLUDED 

*  LOCK  9  -  ADD  UNCERTAINTY  OP  TARGET  LOCA¬ 

TION  INTO  ANTENNA  SCANNING. 

1  -  LOCN  ANTENNA  ON  TARGET 

5-7  IXA  INITIAL  INTEGER  STARTING  LOCATION  DP 

ANTENNA  X  COORDINATE,  MET  ESS. 

9-10  PXA  PINAL  INTEGER  ANTENNA  LOCATION,  NETESS 

11-13  INCXA  INTEGER  BY  WHICH  IA  IS  INCBEMSNTED 

PROS  IXA  TO  PXA.  IP  RUN  DSSS  ONLI 
ORE  ANTENNA  LOCATION,  SET  IXA*PX1, 
INCXA* ANY  NON-ZERO  P05ITITS  INTEGER. 

14  IPILSS  0  -  NO  OOTPOT  TO  DATA  PILES 

1  -  OOTPOT  TO  DATA  PILES 
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C  ••  NOTZ:  PROPER  JCL  ROST  32  INSERTED  •• 

C 


c 

1  5 

IP3IHT 

0  -BO  DETAILED  OOTPOT  OF  SOLOTIOi 

c 

1  - 

DETAILED  OOTPOT  OF 

SOLUTION 

c 

c 

1  6 

ELSCAB 

1  - 

ELE7ATIOH  SC1B  HOT 

OSED. 

c 

2  - 

ELEVATION  SCAB  OSED 

c 

c 

17 

R  ETVOL 

0  -* 

OFF 

c 

1  - 

DISPLAT  SCAB  70LTAG 

ES,  ELE7ATIDB 

c 

ABD  THEN  AZIR07H 

c 

19-19 

AZDIfl 

ODD 

INTEGERS  BETWEEN  6 

AND  SO 

c 

20-21 

2  LOIR 

ODD 

ISTEGEaS  BETWEEN  6 

ABD  50 

c 

22-24 

FLTIHE 

c 

c 

c 

c 

OB  THE 

BEIT  LIBE,  EMTEB 

0  03 

1  POa  TBR SSHOLDISG 

SSTHOD: 

c 

c 

1 

RE  IB 

0  - 

OFF 

c 

1  - 

CALCOLATE  SEAS 

c 

c 

2 

R  EDI AH 

0  - 

OFF 

c 

1  - 

CALCOLATE  HEDIAN 

c 

c 

3 

BADA3CG 

0  - 

OFF 

c 

1  - 

CALCOLATE  WEIGHTED 

ANGLES 

c 

c 

4 

TWEL7E  D 3 

0  - 

OFF 

c 

1  - 

OS 

c 

5 

SECOBD 

0  - 

OFF 

c 

0EBITATZ7E 

1  - 

ON 

c 

c  — • 

»•  «*  «< 

c 

CALL  DATS  (TODAX) 

■  BITS  (6,  1)  (TOOAr(I)  ,1-1,3) 

1  POaHATC  1  ',9A4) 

BEAD  (5,100)  ICODE,  BOISE,  IHPT,  LOCK,  IX A , PXA , I HCXA , IFILES, I  PS  I  NT, 
•ELSCA  N,  BET70L,  AZDIR,  SLDXB,  FLTIHE 
BEAD  (5,2)  (OPTION  (I)  ,  I- 1,5) 

2  FOai!A?(SI1) 

■  BITS  (6,#) 

4  FOB  II  AT  ('  THE  FOLLOWING  I BPOB  RATION  iAS  BEAD  IS') 

■  BITS  (6,5)  ICODE,  BOISE,  IBPT.LOCS 

5  POaSATC  ICODE-  ',11,  '  BOISE- •  , 1 1 ,  •  IBPT-',I1,'  LOCK-*  ,1  1) 

■  BITS  (6,6)  I XI, P X A, IN CXA, IFILES, IPS! NT, ELS CAN,  BET  TO  L 

6  FOSRAT  (•  IX  A- '  ,  13,  2X,  •  PXA-'  ,13,  *  TNCXA-',I3,'  IFILES-',I1, 

••  IPBIST-* ,11, •  ELSCA S— '  ,12 , •  8E?70L-',I1> 

■  BIT  S  (6, 7)  AZDIH,  EL  DIB 

7  FOBSAT  (•  AZDIB-',I2,2I, 'ELDIH-',I2) 

■  HITE  (S,9)  FLTIHE 

9  POBH  AT  (21,  '  FLT  IKE-'  ,14) 

■BITS  (5,  10)  (OPTIOB  (I)  ,1-1,3) 

10  FOSSA;  (’  REAS-  ',11,'  RED  I  IB-  ’,11,'  2ADABC-  ',11) 


uuuuuu  uou  u u u  uuu  u  u  u  ouo 
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NHITS  (6,11)  (OPTION  (I)  ,1*4,  5) 

11  POSH  AT ('  TS2LVE  03*  ',11,*  SECOND  DE3I7ATI7E*' ,11) 

KBIT!  (6,103) 

THE  FOLLOWING  LIST  OP  7ABIA3L2S  CAN  32  CHANGED  TO 
03S2BT2  THE  INFLUENCE  OP  PABTICULAB  FACTORS  ON  T3ACKING 
PSaPOBHANCE. 

BANOOa  GZN2BATOB  SEED  ' 

IX-6554  9 

ANTENNA: 

NTI*  3 

AZBN*. 916 

EL37-.  336 

TA— 500.  0 

HEIGHT*  3. 0 

C-3.0E*09 

FS2Q-9. 1 !♦ 03 

PNHT— 1 0. OS*  05 

LANBOA«3.0E*08/F3S3 

(C-2.0-PI/LANBDA 

AZSPC-1.0 

ELS?C“1. 0 

AaTENNA  ONCEBtAINTt  FACTOR  -  GI7ES  */-4  N  2H30B  AT  P.ASGS*620  3 

DELTA-ATAN2 (4..620.) 

IP (LOCK. EQ. 1 ) DELTA-0. 

GHANA*  (AZ8H*0EGa AD* DELTA)  *2.0/ (FLOAT  (AZDIB)-  1.) 

G3ANE-  (ELBN*0EG3 AD+DBLTA)  *2.  0/  (FLOAT  (ELDIH)-l.) 

TASGET 

SPEED-14  9.63636 
SPEED-SPEED*. 447 0399 
GLSLP-3.5 
f NGTGT-3. 72 
PERIOD*  3.  0 
PRATE-5. 0 
THAT?- 5. 0 
BRATS- 1 0. 0 

SCANNING  TINE: 

TIHEI-0.  0 
PIP- 6000 

TIBEP-  (1. 0/P2P)  -FLOAT  ((IT I) 

TPHSPT-1 0B-06 

TISE-1.0  -  (TIBEP*  (ELDISMELDI3-  1)/2)  ♦  (AZDIB *(AZDIN-  1)  /2)  *TPRSPT) 

CONVERSIONS. 

BASF  *9. 1 E*09 
GULP— 6LSL  P-DEGSAD 
RNGTGT-INGTGT* 1853. 2288 


onnn  noon  non  noonoooooo 
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CMS3A*2*PI/PE8IOD 
PRATI«P8A?E-0EGPAD 
R8ATE-8HATE-DS38AD 
T8ATE-T8AT  E'DSGSAO 
SN8SL- 1 S 75 
S  NH  AZ- 1 E75 
INIT-0 


SAIN  LOOPS  THIS  LOOP  CONSTRUCTION  CAUSES  THE  AZX33TH  (A)  AMD  THE 
ELEVATION  (E)  TO  CHANGE,  SCANNING  THE  ANTENNA.  NOTE  THAT  THE 
SIMULATION  SCANS  PI8ST  IN  AZIMUTH  8ITH  S-0  ( 3 E A M  ON  TAP GET  IS 
ELEVATION) ,  AND  THEN  SCANS  IN  ELEVATION  JIT3  A-0  (SEAM  ON 
TARGET  IS  AZIMUTH). 


DO  52  JXA-IXA,FXA,INCXA 
XA--FLOAT  (JXA-1) 

NSCANS-0 
SNHEL-999.99 
SNRAZ-999.99 
NSCANS-0 
TI8EI-0. 0 
JUST  I -0 

DO  13  JK-t,AZDI3 
VLT  (1.JK)  -0.0 
1  3  VRMS  (1.JX)  -0.0 
DO  14  JX-1.2LDIM 
VLT  (2,  JK)  -0.  0 
19  78 MS  (2.JK)  -0.0 

NRITS  (6,15)XA,  TA,  HEIGHT 

15  FORMAT ( ' 1 ' » '  ANTENNA  POSITION  NITH  RESPECT  TO  END  OF  SURRAT:', 

♦  •  XA-  «,F10.1,«  TA-  '.FIO.!,'  HEIGHT-  ’,F10.3) 

DO  50  SFIT-1,FLTIME 

IF  (MFLT.SQ.  1)  GO  TO  16 
TIHEI-TIHEI-TI3E 

16  CONTINUE 

INCREMENT  TIME  IF  ELEVATION  SCAN  IS  NOT  CALCULATED. 

IF  (ELSCAN.  RE. 2)  ELSCAS-1 

IF  (  (3PLT.GT.  I)  .AND.  ELSCAN.  ZQ.  1 )  TISSI-TI MEI  »TIMEP* 

•  (ELDIM*  (SLDI3-  1 )  /2)  ♦T?HSFT-({AZDI3”l)/2) 

DO  31  JSCAN-1, ELSCAN 

SET  A  a A  NOON  NU33E8  FOa  ANTENNA  POINTING  UNCERTAINTY.  IF  LOCKED 
ON,  CALLING  8 AN  DO  NON  RILL  MAINTAIN  THE  SAMS  SET  OF  SCAN  NOISE. 

CALL  8ANDO  (IX,  IT,  IFL) 

II-IT 

IF  (J3CAN.  EQ.  1)  A—  (AZDI3-  1) /2 

IF  (JSC AN. EC« 2.0 a. LOCK.  10.  1)  GO  TO  17 

THE  CSNTSa  SEA  3  POINTING  LOCATION  IS  GIVEN  AN  UNCE8TAI  NTT  OF  */-■» 
8ETE9S  AT  THE  iELEASE  POINT  BT  MOVING  THE  STARTING  LOCATION  OF 
THE  NULL-TO-NQLL  SCAN  RANDOHLI. 
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0CSSTA*2.* (IPL-0.5)  ‘DELTA/GRASA 
A»A*DCSRTA 

17  ir  (JSCAN.  EQ.  1)  E*0 

IP  (JSCAN.  EQ.  1)  DCISTE-Q. 

IP  (JSCAN.  EQ.2)  A*0 
I?  (JSCAN. 2Q. 2)  (JCERTA»0. 

TP  (JSC AN.  EQ.2)  £*-  (ELDIH-  1)  /2 
IP<JSCAN.EQ.1.0R.L0CX.SQ.1)G0  TO  19 
0CSPTE*2.  •(TPL-Q.S)  •OELTA/GRANE 
P*E»'JC2RT! 

19  IP  (JSCAN.  2Q.  I)  DIS-AZDIS 
IP  (JSCAN. 2  Q.  2 )  0I.1-ELDI3 
DO  30  II-1,DIN 

IF(JSCAN.  SO-  1.  AM  D.  II.  GT.  1}  TIB  EI-TIS2I*  TPHSPt 
If  (JSCA M. 23.2.  AND. II.  G?>  1)  TIME  1**1  RE I+TIS2P 
IP  (JSCAN.  EC.  I)  AZ  2  ADS  (II)  *GRANA*A 
19  CONTINUE 


SOBSOOTIMS  TRAJ:  THIS  SOB BOOT  IN  E  GEMESATES  THE  C0B8ENT  LOCATION 
OP  THE  TA.3S2T  ON  A  GIVEN  TSAJECTORr.  ALSO  GENERATED  ARE 
PERIODIC  VALUES  OP  ROLL/  PITCH,  AND  TAN,  CORRESPONDING 
TO  LOCAL  TORBGLENCE  AND  N IN D  SHEAR. 


CALL  TBAJ (TI.1EI ,  1 3P  T) 

IP  (TO.  LE.  90.0)  GO  TO  20 
GO  TO  21 

20  «aiTE(6, 130)  TI3EZ, NSCAN S 
GO  TO  51 

21  COSTINGS 


PIRD  THE  ANTENNA  LOCATION  IN  THE  TARGET  COORDINATE  STSTE3. 

THESE  ARE  NOT  SPHERICAL  PHI  AND  THETA,  3G7  RELATIVE  ANGLES 
PROS  TARGET  TO  ANTENNA.  THE  CENTER  OP  THE  ENTIRE  COORDINATE 
STSTEH  IS  THE  R0N1AI  TOOCHDOWN  POINT,  NOT  THE  ANTENNA  LOCATION. 


CALL  ROTAT2  (IA,TA, HEIGHT, X1,T1,Z1) 

PHI*- ATAN2  (XI,  II)  *91/2.0 

THETA* ATAN2  (SQST  < X  t  ••  2  ♦  T  1  •  « 2 )  ,Z1) 


THE  9ACKSCATTEH  CROSS  SECTION  OP  EACH  SCATTERER  IN  THE 
TARGET  CONPLZX  IS  NOR  DETERSINSD  POE  A  GIVEN  ANTENNA 
POSITION  IN  TARGET  COORDINATES. 


CALL  RCS(?H2TA,PHI,XC0DE) 
IP  (IPRINT.EQ.O)  GO  TO  22 
•TRITE  (6,132) 

RSITEI6,  12  1)  TIflEI 
RRITZ  (6,  122)  A,  S 
RUTS  (6,  1  09) 

KBITS  (6,  129)  Z0,X0,Z0 
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P30«ND  ‘3ADDEG 

PGiS31*Gi.1.3i»3ADDIG 

PCHI*CHI»8ADDEG 

WRITE  (6,  125)  PSO,  PGA3HA,  PCHX 

BRITS  (6,  104) 

PP«PHI*HADD2G 
TT«THETA*RADDEG 
BRITS  (6,115) 

WRITE  (6 , 1 1 6)  PP.TT 
WRITS  (6,  100) 

SCAT?«10.0«ALOG10  (SC1T  (1  ,4) ) 
SCATS8»10.0*ALOG10 (SCAT (2,4)  ) 
SCATLW*  10.  0*ALOG 10  (SCAT  (3,  4)  ) 
WRIT!  (6,  112) 

BRIT!  (8,1 13)  SCAT?, SCATS B, SC AT  LB 
22  CONTINO! 


FIND  TARGET  ORIGIN  IN  RADAR  COORDINATES 


XDIP-DBL!  (XO)  -D3LE  (XA) 

TDIF-DBLE  (YO)  -DaLE(IA) 

ZDIFaDBL!  (ZO)  -DBLE  (HEIGHT) 

GRNG»SQST  (SNGL  (XD1F)  ••2*'SHGL  (IDIF)  »»2) 

GNDRNG'GRNG 

RAHGE-SQSr  (SNGL(XDIF)  »*2  ♦S8GL  (TDIP)  **2»SMGL  (ZDIF)  »»2) 
HT«SNGL(ZDIF) 

AZ0-ATAN2  (SNGL  (XDI7)  ,  SNGL  (YDI?)  ) 

EL0- ATAN2  (SNGL  (ZDIF)  ,  GRNG) 

IF  (II.  HE.  (  (DIE*  1)  /2)  )  GO  TO  90 

PAZ*  AZO 

PEL-ELO 

C3LOCA- AZ3«G3ANA«0CE8TA 
CBLOC5-ELO*GRANA*OCZBTE 


COHPOTE  ELEVATION  AND  AZISOTH  ANGLES  TO  THE  3EAH  POINTING 
POSITION  WITHIN  IRE  SCANNING  WINDOW. 


)0  AZ1»AZ0*G8ANA*A 
BL1«SL0»GHANE»S 
IF  (IPRINT.EQ.O)  GO  TO  23 
WRITE  ( 5, 1  1  7) 

PEL0*EL0*RADDEG 

PAZ0»AZ0*RADD!G 

WRITS  (6,119)  PELS, PAZO, RANGE 

WRIT!  (6,  1  OB) 

23  CONTINUE 
TDIRCT-O.O 


SOB-LOOP:  THIS  LOOP  SOHS  THE  BACKSCATTEB  FRCS  EACH  OF  THE 
SCATTSRERS  IN  THE  TARGE?  CORPLEX. 
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00  23  I*1,NUfl 


FIND  SCATT  SEES  IN  RADAR  COORDINATES 


CALL  BOTAT  (SCAT  (I,  1)  ,SCA?(I,2)  ,SCAT(I,  3)  ,XD,XD,ZD) 
XDIP*XD-XA 

TDIP«TD-TA  t 

ZOI?*ZD-HEIGH? 

BNG  aSQBT  (SNGL  (XDI?)  «*2*SNGL  (TDIP)  *»2»SNGL  (SOI?)  »»2) 
GBNG*SQRT  (SNGL  (XDIP)  »«2»SNGL  (TDIP)  «*2) 


C  C0HP0T2  AZI NOTH  AND  ELEVATION  ANGLES  TO  SCATTEREB  IN  3 ADAS 

C  COOED! NAT  IS . 

C 

. . 

AZ»ATAN2(  SNGL (XOI P)  .SNGL(VDIF) ) 

£  Lx  AT  A  N2  (SNGL(ZDIP)  ,  G8NG) 

IP  (I  PRI  NT. SQ.  0)  GO  TO  2» 

WHITE  (6.  119) 

pel»el«raddeg 

paz«az»baodsg 

WHITS  (6, 120)  I.PEL.PAZ.RNG 
WHITE  (6,  1  04) 

2«  CONTINUE 
AZD-A21-AZ 
ELD-EL1-EL 

CALL  A NTSNA (AZD.ELD.DGAIN) 

. . . 

c 

C  SBSCT  IS  THE  COUPLER  E-PISLD  POH  THE  SCATTE2EH  OP  INTEHEST 

C  IN  THIS  PARTICULAR  SCAN.  NOTE  THAT  THIS  FIELD  IS  RIGHT 

C  OS  LEFT  CIRCULARLI  POLARIZED  BACKSCATTEB,  AND  BECAUSE  OP 

C  ANTENNA  POLARISATION  SESSITIVITI,  APPR0XI3ATELT  aALP  THE 

C  9ACNSCATTS3  POWER  IS  AVAILABLE  FOR  PROCESSING. 

C 

. . . 

RGT*  3NG*,2*4.0»PI 

£BSC7*LA3B0 A*SQBT  ( (SCAT ( I,  »)  «PWRT« 377.  ) /RGT) 
VDIRCT*VDIBClaDGAIS» (E8SCT«*2) 

25  CONTINUE 

CaaaaaaaaaaattataaMaaaaaaiaa*aaaaaaaaaaaaaaataaaaaaaaaaaa.aaa>aaaa 

C 

C  73HS  IS  THE  SOS  OP  THE  C03PLSZ  SACKSCATTEH  70LT AGES 

C  INDUCED  31  EACH  OP  THE  TARGET  C03PLEX  SCATTERERS. 

C 

Caa«aaaaaaaaaaaaaaaaaaaaaaaaaaa<aaaaaaaaaataaaaaaaaa<aaaaa.<aaaaaaa 

IP  (IPBINT.EQ.O)  GO  TO  26 
PD*  ( AB  S  (TDI8CT)  •*2)*10E*05 
PGN*  1  0.  Q*ALOG  1 0  (DGAIN ) 

WRITE (6, 127)  TDIRCT,PD 

PEL1*EL1 •RADD2G 

PAZ 1* AZ 1 *R ADDSG 

WRITE  (6, 126)  ?GN,?2L1,PAZ1 

WRITE  (6,  104) 


T 

l 


I 

1 


i 


o  n 


2S  COST  I N3S 

vtT(Jscis,ri)  *sroiac? 

IP  (JSCAN.  EQ.  1.  ASD.  II.  EQ.  ( (  A2DI3  ♦  1)  /2) )  SIGNLA *7DI SCT 
IP  (JSC AS. SQ. 2. ASD. II. EQ.  ( (ELDI3*  l)/2) ) SIGNLE«7DI2CT 
IP  (NOISS.  EQ.O)  50  TO  23 

c« .............. 

c 

C  ADD  GADSSIAS  NOISE  USING  "RASDD"  ASD  "GAUSS" 

C 

C  IS  DB  OP  BOISE  IS  ADDED  AS  FOLLOWS: 

C  IP  (SIS N AL  A3PLITUDS/RNS  NOISE  70LTAGE) *SQRT (2«SNR)  ,  THEN  BITS 

C  SISSAL  A  3?L ITU  DE*  3.  1  (AT  A  REIATIVZ  HAXISUS  AT  BASSE)  , 

C  BBS  SOI SE  7OLTAGE»S«0.  3399 

C 

S*0.  3398 
A  3*0. 0 

c«. 

C 

C  GAUSS 

C  COSPOTES  A  NORMALLY  DISTRIBUTED  RASDCH  SUS3ER  WITH  A  SI7ES 

C  SEAS  ASD  STANDARD  DEVIATION. 

C 

C  DESCRIPTION  OP  PARAMETERS 

C  IE  -IE  BOST  CONTAIN  AN  ODD  INTEGER  SDS3E2  7ITH  NINE  OR 

C  LESS  DIGITS  ON  THE  PIRST  ENTRT  TO  GAUSS.  THEREAFTER 

C  IT  RILL  CONTAIN  A  UNIPORHLY  DISTRI3UTED  INTEGER  3 A S DOS 

C  NOBBEH  GENERATED  BT  THE  SUBROUTINE  FOR  USE  ON  THE  BEET 

C  ENTRT  TO  THE  SUBROUTINE. 

C  S  -THE  DESIRED  STANDARD  D27IATI0M  OF  IRS  N0H3AL 

C  DISTRIBUTION. 

C  AH  -THE  DESIRED  BEAN  OP  THE  NOPSAL  DISTRIBUTION. 

C  7  -THE  7ALUS  OP  THE  COHPUTZD  N0R3AL  3AND0B  7ABIA3LE. 

C 

C  RE3ARKS 

C  THIS  ROUTINE  USES  B1SDU  BHICH  IS  STST33/360  SPECIFIC. 

C 

C  SUBROUTINES  REQUIRED 

C  RASDU 

C 

C  HETHOD 

C  USPS  12  UNI PORN  RASDOH  NU83SRS  TO  C03PUTS  NORSAt  R ANDO 3 

C  RUBBERS  BT  CENTRAL  LIBIT  TBE0RE3.  THE  RESULT  IS  THEN 

C  ADJUSTED  TO  NATCH  THE  GI7EN  SEAN  AND  STANDARD  DE7IATI0S. 

C  THE  ONIFOaa  3 A N DOB  SUBBERS  COHPUTED  WITHIN  THE  SUBROUTINE 

C  ARE  POUND  31  THE  POWER  RESIDUE  SETHOD. 

C 

C««. a......... 

AR*0. 0 

DO  27  IS«1,  12 

CALL  BAIDU (II, IT, TPL) 

IE*  IT 

27  AR-AR.TPL 

7*  (AR-6.  0)  «S*AH 
7NOISE  (JSC AN . II ) >7 

NO*  ADD  NOISE  TO  70LTAGE  RETURN 
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7L?  (J  SCAN,  II)  *VLT  (J  SCAN,  II)  *7 
IF  (IPSINT.EQ.1)  KBITS  (6, 101)71?  (JSCAN.II) 

IF  (NOISE.  2Q.  1)  GO  TO  29 

CTCLE  BANDO  IF  NOISE  IS  OFF  TO  flAINTAIN  SANE  ANTENNA  UNCERTAINTY 
FOB  ALL  FLIGHTS. 

DO  29  IN*1 , 12 

CALL  BANDO  (IX,  IY.YFL)  4 

29  IX»IT 

29  I?  (JSCAN.  2Q.  1)  A»A*AZSPC 
IF  (JSCAN.2Q.2) £*2*ELSPC 

IP  ELE7ATI0S  SCAN  IS  0S2D,  AD7ANCE  TINE  TO  SOYS  BEAN. 

IF  (JSCAN. SQ.  1.  AND.ELSCAN.  EQ.  2.  AND.  DI N.  20- 1ZDI3)  TIS2I*TINEI* 
•TIHSP»(  (2LDI3-1)  /1\  *TP  HS  FT«  ( (A  ZOI «-  1 )  /  2) 

30  CONTINOE 

31  CONTINOE 
¥SCANS*NSCAN3*1 

IF  (NOISE. EQ.O)  GO  TO  35 


CALCULATE  SIGNAL  TO  NOISE  BATIO 
AZINOTH 


7NS  Q-0.0 

DO  33  JJJ»1,AZDIB 
7 NOISE  (1  ,JJJ)  *7N0ISE(1  ,JJJ)  «*2 
33  7NSQ*?HSQ*YS0IS2 (1, JJJ) 

7N8SS»SQHT (7NSQ/AZDIH) 

SNBAZ-O. 5«  (SIGNLA/YNR.NS) »*2 
SNRAZ-10.  0*ALOG10  (SNSAZ) 

IF  (ELSCAN.  EQ.  1)  GO  TO  35 

ELE7ATI0N 

7KS  0-0.0 

DO  3#  JJJ*1,2LDia 
7  NOISE  (2  ,  J  J  J)  >7  NOISE  (2, JJJ)  **2 
39  7N5Q*  7NS 0*  7 HOIS E  (2,  JJ  J ) 

YNBHS-SQBT (YNSQ/ELDIH) 

SN8EL*0.  5*  (SIGN!  2/7  NS  NS)  **2 
SSBEL«10.0*ALOG10 (SNBEL) 

35  IF  (  (OPTION  (1)  .  EQ.  1.  OB.  OPTION  (2)  .  EQ.  1.03.  OPTION  (3)  .EQ.  1.03. 
•OPTION (»>. IQ. 1. OB. OPTION (5).EQ. 1.) 

•.  AND.  J0ST1.EQ.  1.AVD.IPBINT.SQ.  0 >  f  H IT2  (6,  1 0  5) 

DO  36  1*1 , AZDIN 

36  AZI  (I) -DBLB  (7LT  (1,1)) 

DO  37  I*1,ELDIS 

37  ELT  (I)  *03 LE  (TLT  (2 , 1)  ) 

IP  (BET70L.  EQ.  1 )  B  BITE  (6,  133)  (ELY  (III)  ,  1 1 1*  1 ,  SLOI.N) 

IF  (BET70L.SQ.  1 )  N3ITE  ( 6 , 1  0#) 

IF(BET70L.EQ.1)BBIT2(6,133)  (AZI  (III) , III* 1 , AZDIN) 
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I?  (OPTION  (  1)  .EQ.O. AND. OPTION  (2)  .  EQ.  0  .  A  ND.  OPTION  (3)  .ED.O.  AND. 
■OPTION (ft) .EQ.J. AND. OPTION  (5) .10. 3)00  TO  45 
IF(J0ST1.EQ.  I.AND.IP8INT.EQ.  1)  GO  TO  39 
IP  (JUST1.EQ.  1)  30  TO  44 
JOST 1  *  1 
WRITE  (6, 33) 

38  FORMAT  <//,  54X,  •  CENTROID  MEASURES  ENTS' , //, 

**  MOTE:  A  MINDS  SIGN  ON  AZIMUTH  E8B08  MEANS  THE 

•  •CALCULATED  TA8GEI  POSITION  IS  T(J  THE  LEPT  OP  THE  ACTUAL', 

••  LOCATION.',/, 

*'  A  HINDS  SIGN  OS  ELE7ATI0S  ZBBOS  MEANS  THE  ’, 

••CALCULATED  TA83ET  POSITION  IS  3EL0W  THE  ACTUAL  LOCATION.', 

•/,'  THE  FOLLOWING  APPLIES  TO  120B  AND  SDaV:',/, 

•61, 'THE  TWO  NUMBSaS  PHECSDING  THE  .METHOD  NAME  A2E  THE  MODE  ', 

•  ’ P LA 3 S  IS  ELEVATION  AND  AZIMUTH,  8ES  PECTI7  ELI .  •  ,  /,  81 , 

• • P08  I2DB,  MODE-O  MEANS  ALL  8ETU8NS  IN  THAT  SCAN  WE3E  A  BO  7  E  THE  •, 
• • 1 2D3  TH8ESHOLO.  MODE-1  MEANS  AT  LEAST  ONE  3STUSN  WAS  \B07S  THE', 
«/,6I, ‘THRESHOLD.  P08  SDP7 ,  30DE=0  MEANS  THE  SCAN  WAS  REJECTED  ', 
••DUE  TO  NOISE  OR  SHAPE  CRITERION.  MODS-2  MEANS  ?H2',/,6X, 
•■CENTROID  WAS  CALCULATED  BT  FINDING  THE  MAXIMUM  SLOPES  ON  ', 

••THE  asrURN. ',/, 

*6X, 'THE  TWO  TWO-DIGIT  NUMBERS  POLLOWING  THE  METHOD  NAME  ARE  THE', 
••  NUMBER  OP  PULSES  USED  IN  THE  ELEVATION  AND  AZIMUTH  SCANS,',/, 

•6X , ' 3ESPECTI7  ELT • ' ) 

9  «aiTE(6,40) 

4  FORMAT  (/,  1  3  1  (*_')  ,/,48X,' 1' ,33X,*  |  •,/) 

WHITE  (6,4  1) 

4  1  FORMAT  (•  ♦•  ,2  IX,  •  ACTUAL*  ,20X,'  |  •  ,  12X,  •  CALCULATED*  , 

•  111,*  I  • ,221, 'ERROR* ,/,48X,‘| *,33X, »| •) 

WRITE  (6, *2) 

42  FORMAT (•♦• ,130 (•_•),/, IX,' aANGE  IN  METERS  SNR  IN  D3  ELEVA 
•IX, ' AZI3UTH  | ' ,4X, 'METHOD  ELEVATION  AZIMUTH  |',3X,'EL27A 
•7X,  •  )  •  ,9X,  '  AZIMUTH'  ,/,  1  X  ,'  ANTENNA  P.D  NWA I  •  ,  5X  ,  •  EL'  ,  4X  , 

•  •AZ  DEGRESS  DEGREES  | ' , 1ST, • DEGR  EES  DEGREES  I  MILLI2ADIANS  ', 

•  •  METERS  |  MILLI2  ADI ASS  METERS') 

WRITE  (6,43) 

43  FORMAT  (•  ♦ '  ,  1  30  (•  ')  ) 

44  IF  (OPTION  (1)  .EgTl)  CALL  MEAN 
IF  (OPTION  (2) .EQ.  1)  CALL  MEDIAN 
IF  (OPTION  (3). 20.  1)  CALL  SADASC 
IF  (OPTION  (4).  EQ.  1)  CALL  T12DB 

IF  (OPTION  (5)  .Eg.  1)  CALL  SDR7(ISIT, NOISE) 

ISIT-1 

IF  (I PRINT.  EQ.  1)  WRITE  (6,  102) 

C 

C  WRITE  TO  DATASETS 

C 

'45  IF  (I  FILES.  EQ.  0)  GO  TO  50 

W8ITE(10,106)XA,RANGS,X8EAN,XMSD,XRADBC,X12D3,ISDR7,SNRAZ, 

•PLS 1 2 , PLSFD, PLSSD 

WRITS  (1 1 ,106)  X  A,  RANGE,  OH  EA  N ,  DN  ID ,  DR  A  DRC ,  D 1  2D3  ,  D  SD  R  7 ,  SNP  A I , 

•PLS12, PLSFD, PLSSD 

50  CONTINUE 

51  CONTINUE 

52  CONTINUE 
STOP 

150  FORMAT  (4 1 1  ,313 , 4 1) ,  2 1 2,  1 13,  F5.  2,  F  3.  1) 
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13  1  P0RMAT(5X,'  RETOEN  VOLTAGE  WITH  NOISE*  ’,215.3) 

132  format  (•  i  <) 

133  FORMA?  ( 1  -  * ) 

134  FORE  AT  ( '  0 ' ) 

105  FOR  SAT  (#3X.*  |  '  ,33X,  '  J  •  ,2#X,'  P  ,/,43X,'  |  ',3  3X,  •  |  •  ,2<*X,  •  p  ) 

136  FORMAT (9E 16.3, 312) 

112  FORMAT  ('O'  ,  1(3X,  ’RCS  OF  FOS  EL  AG  Z,  E  IGH?  WING, AND  LEFT  WING  (D3SM)'! 

113  FORMAT(*0','  F3SELAGE*' ,?10.4,  •  R  H  W I NG  *  •  ,  ?  1 3  .  4 ,  •  LHi  I  S3*  <  ,  F 1  3.  A ) 

115  FORK  AT  ('  O'  ,101,  *  ANTENNA  ASPECT  ANGLES, WITH  RESPECT  TO  T3  E  ’, 

• ‘TARGET’) 

116  FORMAT  ( •  •  ,  15X, 'PHI*' ,P1  2.  9,  IT, 'THETA*  • ,  PI  2.  9) 

117  FORMAT ( ' O'  ,13X, 1  I  An GET  CENTER  LOCATION  IN  RA3AR  COO R 01 N AT ES • ) 

119  FORMAT ( '  ', 15X, 'ELEVATION*', riO. 6, IX, 'AZIMOTH*' ,F  10.6, IX, 

•  '  3ANG2*‘  ,  F  1  0.  4 ) 

119  FORMAT  0‘ , 1  OX ,‘ SCATTERSR  LOCATION  IN  RADAR  COORDINATES’) 

120  FORMAT  ( '  • ,  1  OX .  '  SCA  T«  •  ,  1 2 ,  1 X  ,  <  ELEV AT  10  N-  • ,  F 1  0.  6 ,  1 X  , 

«  •  AZIMUTH*’  ,  P  10.6,  IX,  '  RANGE**  ,  F  10.  4) 

121  FORMAT  (’ O’ ,  10X,  •  SCANNING  TIME*  •  ,F  1  2.  6,  '  MILLISECONDS’) 

122  FORMAT  (•  O’  ,  1  OX,  •  AZ  6  EL  REAM  CORRECTION  NOSBERS  FROM  EOR  ESI3ET:  • 

•  •  A*  '  ,  F5.  2,  '  E* '  ,  F5.  2) 

124  FORMAT ('O', ST, 'CURE  ENT  LOCATION  OF  TARGET  CENTS  :  X* * , F 1 2 . 7 ,  IX , 

•  'T=‘  ,P12.7,  IX.'Z*',  FI  2.7) 

125  FORMAT < » Q • ,5X, 'CORRENT  VALOIS  FOR  TARGET  ANGOLAS  MOT! ON ( DEGS EES) ; 
••  POLL-’ ,F6.3, IX,  'FITCa»‘ ,P5. 3, IX, 'TAR* • , F6. 3) 

126  FORMAT  (’ O' ,5X,  '  GAIN-' ,  F12.  8, ‘OB  AT  EL*  ' ,  F  I  0.  6,1'  6  AZ*',?10.6, 

•2X. ' (DEGREES) • ) 

127  FORMAT  (' O' ,5X,  'RETURN  VOLT  AGE*  •  ,  F  1 6 .  1  3, 2X ,  '  DI RECT  POWER-’, 

•E16.3, IX, 'MICROWATTS’) 

130  FORMAT (' O' ,' FREEZE  COMMA  NO  POINT;  SIMULATION  ENDS.  ELAPSED  TIME' 
•'*  ',F12. 3,  IX, 'SECONDS,  ISCA  NS* ',  I  10  ) 

132  FORMAT  (•  O' ,  100  (•••)  ) 

133  FORMAT  (7F16. 3) 

END 
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subroutine  ramdu(ix,iy,tfl) 

c* ****************************** **************** ***•*••*•********•**••• 

c 

C  PURPOSE 

C  C03P0TES  UNIFORMLY  DISTRIBUTED  RANDC3  SEAL  NUMBERS  BETWEEN  0 

C  ANO  1.0  AND  RANDOM  INTEGERS  3ETSEES  ZERO  AND  2**31-  EACH  ENTRY 

C  OSES  AS  INPUT  AN  INTEGER  SAN 00.1  NUMBER  AND  PRODUCES  A  SEN 

C  INTEGER  AND  REAL  RANDOM  .113323. 

C 

C  DESCRIPTION  OF  PARA3ETE2S 

C  IX  -FOR  THE  FIRST  ENTRY  THIS  EOS?  CONTAIN  ANT  ODD  INTEGER 

C  N033ER  WITH  NINE  CR  LESS  DIGITS.  AFTER  THE  FIRST  ENTRY, 

C  IT  SHOULD  3E  THE  PREVIOUS  VALUE  OF  IT  COS  POTED  3T  THIS 

C  SUBROUTINE. 

C  If  -A  RESULTANT  INTEGER  R  AN  DOS  NU33ER  REQUIRED  FOR  THE  NEXT 

C  ENTRY  TO  THE  ROUTINE.  THE  RANGE  OP  THIS  NUS3ES  IS  BETWEEN 

C  ZESO  AND  2»» 3 1 . 

C  TFL-THE  RESULTANT  UNI  FOSS LT  DISTRIBUTED,  FLOATING  POINT, 

C  RASD03  NUS3ER  IN  THE  RANGE  0  TO  1.0. 

C 

C  R  ESA  3  AS 

C  THIS  ROUTINE  IS  SPECIFIC  TO  SYSTEM/360  AND  WILL  PRODUCE 

C  2**29  TER.1S  3EF08E  REPEATING-  THE  REFERENCE  3EL0 3  DISCUSSES 

C  SEEDS  (65533  HERE),  RON  PR03LESS,  AND  PROBLSBS  CONCERNING 

C  EANDOH  DIGITS  USING  THIS  GENERATION  SCHEIE.  TACLASEN  AND 

C  SARSA3LIA,  JAC3  12,  PP.93-39,  DISCUSS  CONGRUSSTED  GENERATION 

C  SETHODS  AND  TESTS.  THE  USE  OP  TWO  GENERATORS  OF  THE  RANDO 

C  TYPE,  ONE  FILLING  A  TABLE  AND  ONE  PICXING  PROS  THE  TA3LE, 

C  IS  OF  BENEFIT  IS  SOHE  CASES.  65549  HAS  BEEN  SUGGESTED  AS  A 

C  SEED  WHICH  HAS  BETTER  STATISTICAL  PROPERTIES  FOR  HIGH  0? DEE 

C  3ITS  OF  GENERATED  DEVIATE.  SEEDS  SHOULD  3E  CHOSEN  IN  ACC08- 

C  DANCE  WITH  THE  DISCUSSION  GI7EN  IN  THE  RSFE2ENCE  BELOW. 

C  ALSO,  IT  SHOULD  BE  NOTED  THAT  IF  FLOATING  POINT  3  AN  DOR  NU33SRS 

C  ARE  DESIRED,  AS  ARE  1 V1IL ABLE  FRO  H  RANDU,  THE  RANDOH  CHA3AC- 

C  TERISTICS  OF  THE  FLOATING  POINT  DEVIATES  HAVE  HIGH  PROBABILITY 

C  OF  HAVING  A  TRAILING  LOW  ORDER  ZERO  BIT  IN  THEIR  FRACTIONAL 

C  PART. 

C 

C  HETHOD 

C  POWER  RESIDUE  3STH0D  DISCUSSED  IN  133  3ANUAL  020-901  1, 

C  RANDOH  NU  33  S3  GENERATION  AND  TESTING. 

C 

C***************************************** **•****•*•*••*••••****• ****** 

IY*IX*6553 3 
IF  (I  Y)  5,6,6 

5  IT*IY*2147»63647* 1 

6  TFL*I  T 

TP  L*TPL*. 4656613  S-4 

RETURN 

END 
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SU3R0UTI MS  ROTAT<XSCAT,rSCAT ,ZSCAT,XD,  TO, ID) 

. . . . . 

c 

C  THIS  SUBROUTINE  COMPUTES  THE  RADAR  COORDINATES  OP  THE 
C  POINT  (XSCAT.T SCAT, ZSCAT)  GIVEN  IN  TARGET  COORDINATES 

C 
C 

C  ID,TD,ZD  ARE  RADAR  COORDINATES  OP  POIST. 

C  X SC  AT, ISC  AT , ZSCAT  ASS  TARGET  COORDINATES  OP  POINT. 

C  10, TO, 20  ARE  TRANSLATION  COORDINATES. 

C  CHI  IS  TAN  ANGLE  -  POSITIVE  ROTATION  ABOUT  THE  Z-AXIS. 

C  3AM3A  IS  PITCH  ANGLE  -  POSITIVE  ROTATION  ABOUT  THE  I-AXIS. 

C  HU  IS  ROLL  ANGLE  -  POSITIVE  ROTATION  ABOUT  THE  X-AXIS. 

C  THE  ROTATIONS  ARE  ALHATS  PERFORMED  IS  THE  OHDER  -  I A N-PI TCH-ROLL. 
C 

. . . 

REAL  MU 

DOUBLE  PRECISION  XD, T D,Z D,  D3LS 
COMMON/ RTA/MU, GAMMA, CHI, 10,10,20 
1  CG»COS  (GAMMA) 

CC*COS  (CHI) 

CM-COS  (MU) 

SG-SIN (GAMMA) 

SC-SIN  (CHI) 

S  M*  S I  M  (MU) 

SMSG»SN*SG 

CMSG*CH*SG 

XD-DBLE (CG  *CC*X  SC  AT* (SMS G*CC-CM*SC ) *TSCAT* 

♦  (CHSG*CC*SM*3C) *2SCAT)  »DELE(X0) 

ID-DBLE  (C3  *SC*X  SC  AT*  (SMSG*SC  *CM»CC)  *TSCAT» 

*  (CMSG»SC-SM*CC)  *ZSC  AT)  *DBLE ( TO ) 

ZD*0BLZ  (-5G*XSC  AT*  5M»CG*  TSCAT*C M*CG* ZSCAT)  *DBLE(Z0) 

RETURN 

END 
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SU3HOOTINS  30?  A? 2  (X  A,  t  A,  HZ  10  H?  ,  X  1 ,  Y  1 ,  Z  1) 

c* •••••» 

c 

c  this  soshootine  computes  the  Tasgs?  cocpdinatss  op  the 

C  POIS?  (XA.XA, HEIGHT)  0ESC5IBS3  IS  31315  COCHDIS ATSS 

C 

C  Z  I , T  1 , 2 1  APS  TA3GE?  C003DINATES  0?  POINT. 

C  XA,TA, HEIGH?  A  HE  2AD15  C008CINATES  OP  POINT. 

C  XO  »  TO , Z  0  ASS  TEAM  STATION  COORDINATES. 

C  CHI  IS  T AH  ANGIE  -  POSITIVE  SOTATION  A300T  THE  Z-AXIS. 

C  GANNA  IS  PITCH  ANGLE  -  POSITIVE  HOT  AT  ION  ABOUT  THE  Y-AXIS. 

C  10  IS  SOIL  ANGL2  -  POSITIVE  ROTATION  ABOUT  THE  T-AXIS. 

C  THE  aOTATIONS  ASS  ALHAXS  PS8POR.1ED  IN  THE  05DS3  -  X  AH-PITCH—90LL. 
C 

. . . . . 

HEAL  10 

CONNON/8TA/N 0,3 ASHA, CHI,  X0,T0,Z0 

I2-XA-X0 

Y2-TA-T0 

Z2»HEI3HT-Z0 

CG*COS  [GANNA) 

CC*COS  (CHI) 

CN-COS  (NO) 

SG-SIN  (GANNA) 

SC*SIN  (CHI) 

SN*SIN(NO) 

SGCC»SG*CC 

SGSC”SG*SC 

CGZ2»CG*Z2 

X  1»CG  *CC  *X2  *CG  *  SC*  I2-SG*  Z2 

T  1* (SN»SGCC-Ca«5C) *X2»  (SN*SGSC*CH»CC)*T2*S 2*CGZ2 
Z 1» (C  3 •SGCC ♦  SN* SC)  »X2  *■  (C«*SGSC-S3»CC)  •Y2»CH»CGZ2 
HETU8N 
END 


I 


I 

I 

S0B?00TI>!E  T3AJ  (TI3E.IRPY) 

'.SAL  3U 

COUHOH/a  PY/33  ATS,  P3 ATE, SPATE. OREGA. SPEED, GL3LP.H, 
•A2.3SMG 

coafioa/  3ta/3U,ga.vi.\,  cai.ao.ro,  20 

I  PI»3. 141593 

xa«o. o 

ro*cos  (3LSL?)  •  (a-sPssc'Tias) 

20*  TO  *T  A  M  (GLS  L  P) 

IP  (ISPT.S2.0)  00  TO  20 
T*B9*Sia (OR£GA*TiaE) 
au*  (33  ATE/C3EG  A)  *TE8S 
GA3RA*(PBA7E/0SEGA)  *TE8S 
cal*  (T  BAT  2/032G  A)  *TEas 
SET03H 

20  ao*o.o 
CHI*0.  0 
GA3SA*0.  0 
BEToaa 
HMD 


i 

I 

I 

i 

i 

I 

I 

[ 

I 


o  n  n  n  n 
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o 

c 

c 

c 

c 

c 

c 

c 


SOSHOOTIMS  PCS  (THETA,  PHI, ICODE) 

DIH2MSI0H  SCAT<3,») 

COHaOH/SA7E/K,308,ZPS3 
cohsom/scath/scat. moh 
REAL  a 

DATA  PI,  PI2S ,  PI2/3.  1 «  1 59  3,  2.  167 90  1 ,  1. 57079  4/ 


ICS  COMPLEX  CODE 
Is  JETRCDil 

2:  1  POINT  SCAT?  ZaSS  30  CEL 


DO  2  I*  1.3 
2  SCAT  (I,»)  -12-10 

GO  TO  (100,200)  ,  I  CO  0  E 


JET30D 1 1  ECS  TARGET  30DSL 


103  SOB- 3 

ALPHA-A3C0S (SIM (THETA)  *COS  (PHI)  ) 

DELT  A*ASC0S  (SIM  (THETA)  »S  IS  (PHI)  ) 

SCAT  (l.l)-O.O 
SCAT  (1,2)  —3.0 
SCAT  (1,3)  -0.0 
SCAT  (2.  1)  —3-0 
SCAT  (2,2)  >0.0 
SCAT  (2,3)  >0.0 
SCAT  (3,  1)  -*3.0 
SCAT  (3.2)  >0.0 
SCAT  (3,3)  -0.0 
DXP0S«2. 0 
DTPOS-IO.O 
DZPUS-3. 0 
DXL8MG>«.  0 
DTLBSC-6.0 
DZLHNG*2. 0 
DX8f*G>*.0 
DTHIMG-6.3 
DZIB  HG>2.0 

P0SX>C0S ( (K*OXf OS/2.0) >C0S  (ALPHA) ) 
P05T>C0S  (  (K*0TP05/2. 0 )  >COS  (DELTA)  ) 
PUSZ»COS< (K«OZf 3S/2.0)  >COS  (THETA)  ) 

MMGt X>C03  (  (H*0  X8M3G  /2  •  0)  >COS  (ALPHA)) 

i>car>cos  (  (k»dt8mmg/2.0)  >cos  (delta)  ) 

MMGBZ>COS  ( (K*D  ZH1MG/2  . 0)  >C0S  (THETA)  ) 
>IGLX>C03  (  (A ’OILS MG/2. 0)  >COS  (ALPHA)  ) 
»MGLT«COS  (  (K  *D  IL8MG/2. 3)  >COS  (DELTA)  ) 
VIGLZ>COS  (  (K  *D  ZL  MMG/2 •  0)  >COS  (THETA)  ) 
ArOSBL>10.  0>  (THE!  A- PI/ 2.  0)  «  2*  1.0 
A M3GEL- 1  00. 0>  (THETA- PI/2.0)  ••2*1.0 
IP(ASS  (PHI). LE. PI/2.0)  GO  TO  10 
APOSAZ*(7S.O/PI2S)*(ASS(PBX)  -PH  ••2*3.0 
GO  TO  15 


noono  non  non  non 
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10  APCSAZ*  (75.  D/-PI2S)  *?HI**2*a.O 
15  CONTINOI 

AaWGAZ*  1.  0-SIN  (PHI) 

ALWGAZ»1.0»SIN  (PHI) 

COHPOT!  P0S2UGE  2CS 

SCAT  (1,4)  *A?aS2L»AP’JSAZ*A3S  (P'JSX*P0ST*P5SZ) 
COHPOT!  as  WING  BCS 

SCAT  (2,4)  *ANNG2i»  Aa  WG  AZ*  A3S  (  UN  Si  1  *  W  KG  BT*  1 NGSZ) 
CONPOTE  LH  wing  acs 

SCAT  (3, 4)  *A JNG2L*  AN  WG  A  Z»  A3S  (  WN  GLZ  f  V  NG  L 1  *  i  NGL  Z) 
H  STO  B  N 


COHPOTATION  OP  1  POINT  SCATTEBSB  BCS 


200  NON*  1 

SCAT  (1,1)  *0.0 
SCAT  (1,2)  *0.  0 
SCAT  ( 1.  J)  -0.0 
SCAT  (1,4)  -1.3 
•ETOSN 
END 
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SOBS 0071  NS  ANTS'!*  (A2, EL, GAIN) 

REAL  K 

INTEGER  ELA,  ELA2.ELB,  EL32 

coascvRAVE/K.soa, s?sa  . 

DATA  ?  1/3.  1  N 1  5  9  3/ 

DATA  SLA , ELE/35 , 35/ 

DATA  D1, D2.SCAL2/0.03, 0.029, 10.05*05/ 

PHI-PI/2. -AZ 

THETA-PI/2. O-SL 

BI-COS (T8STA) *D1-K 

B 2-5  IN  (THETA)  *COS(PHI)  -D2*K 

R 1-0. 

DO  100  1-1, 2LA 
100  81-8UCOS (31«PL0AT(I) ) 

R 1*2. *81*1.0 
82-0. 

DO  200  1-1, SL3 
200  R2«R2*C05 (82*PLOAT(I) ) 

82-2.  -82*  1 . 0 
8 182-81*82 

GAIN«81H2*R1 82/SC  AL  E 
82T08N 

END 


o  n  a  n  n  n  o 
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S03BOCJTINE  SEAN 


THIS  S0830UTINE  CALCOLATES  THS  SEAN  3ETU3H  VOLTAGE  IS 
AZIStJTH  AS 0  ELEVATION  (SEPABATSLY)  WHICH  IS  T3ES  USED 
AS  THE  THBESHOLD  IS  PHLH  TO  LOCATE  THE  TAHGET  CENTSOID. 


P  !A  L  *  9  EL  (U  9)  ,  AZ  (A?)  ,  S  OH  ,BEA  S3  ,  SEA  SA 

C0H90M  /ALL/  EL,  AZ  ,  SA  S  GE  ,OPT  IOS  ,  L  A  BE  L,  ELD  IS,  AZOIS,  ELSCAH,  DELTA 
INTEGZH  OPTION  (5)  ,  EL D IB ,  AZDIB, ELSCAH 
IP(ELSCAS.EQ. 1) GO  TO  W 
S0H*0.0 
HSANE-0.  3 
DO  2  I* 1, ELDIS 
2  SUH«StJ3*EL(I) 

BEASS»SOH/rLOAT(ZLDIB) 

#  san*o.o 

DO  5  1*1 , AZDIB 

s  sas*sus»  az  (I) 

REANA*SDS /FLOAT ( AZDIS) 

LABEL* 1 

CALL  FHL3  (HEME,  SEA  VA) 

HET03S 

ESD 


n  n  o  o  n  n 
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SU3SOOTIN2  22SIAN 


THIS  SOSPOUTINS  D5T233ISE3  ? HE  BEDIAM  VALGS  OP  SSTtJHN 
VOLTAGE.  AND  PASSES  THIS  VALUE  AS  THaSSHOLD  TO  P3LH. 


C0330N  /ALL/  EL,  AZ.3A8GE,  OPTION,  LABEL,  ELDI3,  AZDI3, EL  SCAN,  DELTA 
C0330N  /8ADB/  PEL , P AZ , C8LOCE.C3LOCA , S3 A NE, GS AH  A, AZ 83 , ELBV, LOCK 
SEAL -9  EL  (9  9)  ,  AZ  (99)  ,  BEDE.  SEDA  .VOLT  (2,  49)  /98-Q.  000/ 

I8TEGE3  OPTION (5) .EL0Xa.AZSI3.C3A3GE.LI3I? (2) .SLSCAN 
L*  2 

:?(ELSCAH.E3-t)L*1 

oo  4  x* i , Azoxn 

4  VOLT  (1,K)»AZ  (R) 

DO  9  X-1 , ELDI3 

5  VOLT  (2.K)  -EL  (K ) 

LI  31?  (1)  -AZDIN-t 
LI  SIT  (2)  -ELOIS-1 

t  CHANGE-0 
00  2  I-1.L 
LL-LI3IT  (I) 

DO  3  J» 1 , LL 

IP(V01T  (t.J)  .LT.VOLT(t.J*n)  GO  TO  3 

VBOLD-VOLT  (I.JJ 

VOL?  (I.J)  -VOLT  (I,J*1) 

VOl?(I.J»1)-VHOLO 
CHANGE -CHANGE* 1 
3  CONTINUE 
2  CO NT  1 8 US 

IP  (CHANGE.  ME. 0)  GO  TO  1 
3I0AZ«(AZDIH*1)  /2 
HIDEL-  (ELDI3»1)  /2 
320 A  * VOLT ( 1 . 91 S AZ) 

HSDE-VOLT  (2.3I0EL) 

LABEL- 2 

CALL  P9L3  (3E3E.8EDA) 

82TU  88 
END 


o  o  n  n  o  o  n 


S03R0UTI N 2  T12D3 


THIS  SUBROUTINE  CONPUTES  *  THRESHOLD  T2D8  DOAN  PROS  THE  COSPUTSD 
RETURN  70LTAGE  WHEN  THE  ANT2NNA  3SAR  IS  OS  THE  TARSZT.  SOTE  THAT 
IP  -  1 2*20 *LOG (7TH RES/ VCZNTER) ,  THEN  VTHP SS*0. 25 11 336 4 *7CEN TEH. 


COSHOS/A LL/S L,AZ, RANGE, OPT  ION,  LABEL,  ELDI S,  AZDIH ,  ELSCA N,  DELTA 
COHBON  /? ADR/  P EL ,P A2 , CBLOCS ,C 3L0C A, GR ANE, 38  AN  A, AZB* , ELBB , LOCK 
COSHON/POO PI7/C0UNTS, COUNT  A, HO  DEE, BODE A, SS DB 7, SNBEL, NODE  12, SODA  12 
INTEGER  OPTION  (5)  ,  ELDI  S ,  AZDIH ,  PHA,  LHA,  PHS.  LHE,  AZS,  ELS,  COONTE, 

•  COUNT  A, ELS CAN 

DOUBLE  PRECISION  EL  (43)  ,AZ  (4 9 )  , 7CSNTR ,  ELTH  RS ,  AZTHRS 
HODS  12*7 
SODA  12*1 

HIDDLE*  (AZDIH*  1)  /2 
TCENTR-AZ(HIODLE) 

AZTHBS»0.2S1 18364«VCESIR 

DO  I  1*1, AZDIH 

IP  (AZTHRS.  GE-AZ  (I))  GO  TO  2 

1  CONTINUE 
SODA  12*0 

2  HIDDLE*  (ZLDIH*  1)  /2 
7CENTR-EL (BIDDLE) 

ELTH RS*0. 251 ia964*7CS*TR 
IP (ELSCAN.EO. 1) ELTHBS*0.0 
DO  3  1*1, ELDIS 
IP  (ELTHRS.GE.SL(I) )  GO  TO  4 

3  CONTINUE 
RODE  12*0 

4  LABEL** 

CALL  PHLH(ELTSBS,A2TH85) 

RETURN 

END 


non  nno  non 
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SUBROUTINE  FHLH  ( ELTHRS  ,  A  ZT  HRS) 

Caaaaa**aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa«a.aaaa«*aaaaaaaaaaaaaaaaaaaaaaaa 

C 

C  THIS  SUBROUTINE  USES  TH2  THRESHOLDS  CALCULATED  3T  BEAN 

C  AND/OR  12DIAN  TO  DETERBINE  THE  TARGET  CENTROID.  AN  EDGE 

C  OR  THE  TARGET  IS  DEPINED  W BEN  T WO  CONSECUTIVE  VOLTAGE 

C  RETURNS  EXCEED  THE  CALCULATED  THRESHOLD.  THE  CENTROID  IS 

C  THEN  BIDWAT  BETWEEN  THE  EDGES.'  A  HIT  IS  A  7QLTA3E  EXCEEDING 

C  THE  THRESHOLD. 

C 

C»*» . . 

DOUBLE  PRECISION  VOLT  (49 )  ,  EL  <4  9)  ,  AZ  (49  )  ,  ELTH  3S  ,  AZTH3 S .  THR  ESH 
COHN  ON  /ALL/  EL,  AZ, R A NGE, 0 PT ION , L ABEL, ELDI 9, A ZDI ■ , E LSCAN , DELTA 
COSHON  /RADH/  PEL ,P A2 , CB LOCE ,C 3 LCCA, GR A N2, Ga AN  A, A2 3W , EL SW , LOCK 
COHHON/ROORI7/COUNTE,  COUNT  A,  NO  DEE,  NODE  A,  NS  DR  V ,  SN  REL  ,  .10  DE  1  2  ,  HODA  1  2 
INTEGER  OPTION  (5)  , ELDI H ,  A  2DIH  ,  FHA ,  LHA ,  PHE ,  LH E ,  AES ,  ELS ,  COUNTE , 
•COUNT  A , ELS C  AN 
COUNTA  *  1 
COUNTE*  I 
HODEA*  3 
BODES*  3 
TCE*0.  0 

CALCULATE  PI3ST  HIT  IN  A2IBUTH 
DO  1  A2S* 1 , A2DIB 

IP  ( (A2  (AZS).GS.  A2TH3S)  .AND.  (A2  (AZS*  1)  .  5E.A2THHS)  )  GO  TO  2 
t  CONTINUE 
PCA-3E5 
HODEA-0 
GO  TO  S 

2  PHA*A2S 
COONTA-COUNTA*  A2S» 1 

CALCULATE  LAST  HIT  IS  AZ  IB  UTS 

DO  3  AZS*1,A2DIH 

IP  (  (AZ  (AZDIS-AZSH)  .  GE.  A  ZT  HRS)  .  A  ND.  (AZ  (AZDIN-AZS)  .  GE.  AZTH3S)  ) 

•GO  TO  4 

3  CONTINUE 
WRITE  (6,  12) 

4  LHA* A ZDIB-AZS* 1 

PCA*  (FLOAT  (LHA+PHA)  /2.-PL0AT  (  (AZDI  N»  1 ) /2)  )  «S  RAN  A  fCBLU  CA 
COUNT  A*COONT  A*  AZSH 

CALCULATE  PIRST  HIT  IN  ELEVATION 

5  IP  (ELSCAN.  ED.  1)  GO  TO  10 
DO  6  ZLS*1,ELDI3 

IP((EL(ZLS).G3.ELTHRS)  .AND.  (EL  (3LS»1).GZ.SLTHRS)  )  GO  TO  7 
S  CONTINUE 
TCS-9SS 
HODEE-O 
GO  TO  10 
7  PHE-ELS 

COUNT  E*COO  NT  Z*  ELS  H 
C 


SIT08M 

,  END 

I 

i 

I 

1 

I 

I 


I 
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SOSRCOTI NE  R  AD  ARC 

C« ............ 

r 

C  THIS  SUB  ROOTI  HE  CALCULATES  THE  A3P  UTUOE-N  EISETED  RADAR 

C  PETURN  TO  DETERMINE  TH2  RADAR  CENTER  OF  SR  AV  IT  T  (RAOAPC)  . 

C  P0R3ULA  USED  IS  (SUM  Of  VOLT  AG  ES  TI.1ES  ANGLE  AT  THAT 

C  701? AGS)  DIVIDED  31  (SOM  OF  THE  VOLTAGES). 

C 

C««« ....... 

COflHON  /ALL/  EL,  AZ ,  R  ANGE.OPT  10  N  ,  LA  BEL,  ELDI3,  A  ZDI  3,  EL5CA  N,  DELTA 
C0H30M  /RADR/  PEL  ,  P  A2  ,  CBLOCE  ,  C  BLOC  A ,  SR  A  N  E,  S3  A  N  A  .  A  23  J  ,  EL  3  «  ,  LDC.< 
DOUBLE  PRECISION  EL  (49 )  ,  A7,  (4 9)  . VOLT  (49 ) 

INTEGER  EL 3 13,  A2 013,  OPTION  <5)  ,DI3,ELSCAN 
PCGEL*0. 0 

IP  (2LSCAN.  EQ.  1)  GO  TO  3 
SUH7TH *0. 0 

saaw*o.o 

j»  i 

DiH'ELora 
CENTER*  (DIH*  1)  /2. 

GRAN-GRANE 
PINGLE’CBLOCE 
DO  2  1*1, DI3 
2  VOLT  (I)  =EL  (I) 

I  DO  10  I* 1 , DI 3 

S0B3TH*SU39TH» VOLT  (I)  *  (  (I-CENS  i ) »GR A N+ ? ANGLE) 

10  SUaW*SUMN*SSOL  (VOLT  (I)  ) 

IF(J.S3.2)GO  TO  20 
RCGEL»SU3* TH/S08R 
3  SUHRTH-0.0 
SUBB-O. 0 
J*  2 

Dia»A  zdi a 
CEMTEB*  (DIH»  1)  /2. 

GRAN*GRANA 
PANGLE*C3L0CA 
DO  11  1*1, DIS 

II  70 LT  (I)»A2  (I) 

GO  TO  1 

20  RCGA2*SUaiTH/S0as 
LABEL*  3 

CALL  ERROR (RCSEL, RCGAZ) 

RETORN 

END 


o  n  n  noon 
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S03ROUTINS  E3303 (ICE, PCA) 


THIS  S 

03300TI: 

(2  CALCOLAIES  THE  23S03 

3  IT  i  E  EM 

THE  LOCATION 

0?  THE 

TARGET 

CENTROID  DETERMINED  BY 

THE  ESTI 

HATOP  S3B300Ti:iES 

TO  THE 

ACTOAL 

LOCATION  0?  THE  TARGET 

IS  5  SIN 

BY  THE  RADAS. 

CO.M'OM  /ALL/  EL, AZ, R ANG2, OPT  ID N, LABEL, ELD 1 *, A  ZD  I N , EL5CA N , DELIA 
C32B0N  /SAD  3/  PEL , ? AZ , C3L0CE , C3L0CA , 33 ANE, 3? AS  A , A 23 « , EL  3 M , LOCK 
COSSON  /FI  LEI  1/  XSEAN  ,  l.NED,X3AD3C,  !(1  2DB,  :(SD3V  ,  S  N  3  AZ  ,  ?LS  1  2, 

•  PLSFD,  PLSS  D,  D3EAN,  DU  ED,  DRADRC,  D  1  2DB,  DSDR7 
CO  MR ON/? 00  RI 7 /CO  ONTE ,  CO'JNTA  , NODES , NODE  A,  SS  037  ,  S  Jf  S  EL,  NODE  1  2 ,  30  DA  1  2 
conscn  /aiA/  fja,3ANSA,cHi,xo,YO,zo 
DOOBLE  PRECISION  EL  (<4  9  )  ,  A  2  (4  9) 

INTEGE3  OPTION  (5)  ,  ELDI.1,  AZOI.H  ,  COONT*  ,  COGSTS,  ELSC A N , ?L S 1 2 , PL S FD, 
•PLSSD 

RADDEG=57. 29578 

C  2330R  CAL CO LAI ION 

C 

ER3EL=TCE-PEL 
ER3AZ«PCA-PAZ 
SLa«EANG2»SIN (2332LI 
AZ3*R A  NGE  *  SI S  (ES3AZ) 

TCED«TCS*3AD0EG 
PCAD«?CA»HADD2G 
P2LD*PSL»RADDEG 
?AZD-*PAZ*3ADDEG 
IF  (RCDEA.EC.O)  ES3AZ=1  E71 

C0N7ERTING  TO  P.ILLI2  ADIASS 

ERRS  L=23  3SL» 1  S3 
EHRAZ«EBRAZ*1E3 

IP  (LABEL.  EC.  1)  raZAN*2FSAZ 
IP  (LABEL. Eg. 2)  I3ED-ER3AZ 
IP  (LABEL.  EQ.  3)  X  SADRC*  ERR  AZ 
IP  (LABEL.  EQ.  4)  X12D3«ER3AZ 
IP (LAB2L.ZQ.5)  XSDRV«E3BAZ 
IP  (LA3EL.  Eg.  4)  PLS 12»CO0N?A 
IP  (LA  a  EL.  Eg.  5)  PLSSD*CO  3NTA 
IF  (LA3EL.  EQ.  1)  03  E  AN  *  A  Z  3 
IF  (LABEL.  Eg.  2)  OBED-AZS 
IP  (LABEL-  Eg.  3)  DaAD8C*AZ3 
IP  (LABEL.  EC. 4)  DUDB-AZ.R 
IP  (LABEL. EC. S)  D  SD  3  7*  A  ZB 

OOTPOT  ERS03  DATA 

IP  (LABEL. Eg.  3)  30DEE-3 
IP(LABSL.Eg.3) 300EA-5 

MBIT  E  (5,  1)  SANG  2,  rO,SNSEL,SNRAZ,P2LD,?AZD 
IP  (LA  3  EL.  E g.  1 )  MSI  TE  (6 , 21 
IP  (LABEL.  Eg.  2)  MRITE  (6, 3) 

I P  ( LA  3  EL.  2  C.  3)  MSITE(S,4( 
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IP  (LA  3  21.  EC*  3)  3  SITS  (6,6)  BO  DE  12 ,30DA  t  2,  COO  STS  ,  C  COSTA 
IF(LA3SL.£C.S)  3aZTE(6,5)  BODEE  ,  BODE  A  ,  COO  STS  ,CCO  NT  A 
I?  (ELSCAN.  EQ.  1 )  SO  TO  II 
IP  (BODEE.  HE.  0)  aaiTS  (5, 7}  TCSD,  ESSEL,  El.3 
I?  (BODES.  E  Q.  0)  3SITE  (6,3) 

M  IF  (BOD  EA.  NS.  0)  V  SITE  (6, 3)  PCAD  ,  ER3  AE  ,  A  EB 
I?  (SODEA.2Q.O)  3  31 TS  (6  ,  10) 

1  PORBAT  <2F9.2,2X.  PS.  1 ,  I  X,  PS  .  1  ,  1  X, PS,  2 , 3X,  PS.  2 , 2X  ,  •  J  • ) 

2  POSBAT  ('  *'  ,39X,  •  SEAR' ,23X,  •  | ’,  23X,  •  |  •) 

3  FOR  BAT  ( '  ♦  '  , 4  3X ,  *  BEDIA  S '  ,2  J  X  ,  '  |  ' ,  23X  ,  •  |  •  ) 

3  POSBAT  ('♦•  ,49X,  •  SADAHCG*  ,22X,  •  J  ’,  24X, ')  ') 

5  POSBAT  <•♦•  ,39X,2H,  U,  'SDRV*  ,1X,2I2,20X,  •  |  ’,237,  <  |  <) 

6  POSBAT  ('  ♦ '  ,  4  3X  ,  21  1,  I  X, '  12D3'  ,  1X,2T2,20X,  *  |  * , 23X ,  '  (  ' ) 

7  POSBAT  (’  *•  ,6  3X,  F6.2,  15X,  G10.  3,  1X.P7.  2) 

3  POSBAT(,»’,73X,?6.2,29X,G10. 3.2X.P7. 2) 

9  POSBAT  (’  ♦•  ,  65x  ,  '*2EJ  ECTED*  1  ,2X  ,  ••REJECT  «') 

10  POSBAT  (’  ♦  • ,  1  03X ,  *  *8  2J  ECTED*  *  ,  3X,  ’•REJECT*’ ) 

RETU3S 

END 


n  n  o  o  o 
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SUBROUTINE  SD97  (INIT,  NOISE) 

CONSON  /A  LL/E  L, AZ , RANGE, OPTION , L A3  EL, E  LD IS, AZ  013 , EL  SC AN , DELTA 
COflHON/RADR/PEL.P AZ, C B LOCE .C3L0CA , GR A NE , OR  AN  A , AZEW, EL  8W, LOCK 
COflflON/FOO  SI7/COO  NT£,CCUNTA,N0DEZ,H0D£A,.'15DRV,SNPEL,30DE12,30DA12 
9  SAL *9  EL  (4  9)  ,  AZ  (4  9)  ,  7REP.7LT  (» >  /»  *0.  D  DO/,  70LT  (2,  U9)  /9B*0 .3 DO/, 
•7END,S0S 

INTEGER  E LD IS, AA0I3, OPTION  (5)  ,  A  ,  3,  C,  D.  01 3,  EL  3C  AN  ,  P  ULS  ES,  Q  , 
•COOMTS.COUNTA, COUNT.  S0DE2,  .NODE A 
SEAL  NOSLIfl.N 

c... ......................... ....... .................................... 

C 

C  THIS  SUBROUTINE  CALCULATES  AN  ESTIMATE  0?  THE  T  AEG  ET  CENTROID 
C  POSITION  3ASE0  CN  THE  SSAPS  07  THE  SET 03 S. 

C  THE  CE  NT  90  ID  IS  CALCDLAT2D  AS  FOLLOWS: 

C  HALT  THE  SCAN  IS  BROKEN  INTO  POOR  'WINDOWS*.  THE  RETURNS  PROS  THE 

C  BEAS  POINTING  LOCATIONS  INSIDE  THE  'WINDOWS  ARE  A7E9AGED  TOGETHE? 

C  TO  OBTAIN  POOH  BEAN  7AL0ES.  THESE  PODfi  HEAN  VALOES  ARE  EXPECTED  TO 

C  TAKE  ON  ONE  HALF  OP  A  BELL  SHAPE.  AN  EDGE  IS  PLACED  AT  THE  POINT 

C  OP  KAXIHUS  SLOPE,  Oa  WHERE  THE  SECOND  DE3I7ITI7E  CHANGES  SIGN.  IF 

C  NO  CHANGE  IN  SIGN  OCCURRED,  OR  IP  THE  KA'X  SLOPE  3  AS  BELOW  THE 

C  NOISE  CRITERION,  THE  WIN  DO  SS  ARE  SHIFTED  AND  THE  PROCESS  REPEATED. 

C  THE  WINDOWS  ABE  ALLOWED  TO  SHIFT  TWO  WINDOW  WIDT3S  PROS  THE  END  OP 

C  THE  SCAN,  oa  TO  THE  FEDGE.  A  SCAN  3EJECTED  DOE  TO  SHAPE  OP  NOISE 

C  IS  TAGGED  f!ODE»Q.  IP  A  CENTROID  IS  CALCULATED,  THE  SCAN  IS  TAGGED 

C  HODS-2. 

C 

c«. 

NOSLIN’O.  3898*  I.  5 

IP  (INI?. EC.  I)  GO  10  30 

7REP-0. 

SUH*0. 

N=0. 

30  VEND-O. 

NODEE-O 

NODE  A«0 

COONTS*0 

COUNTA—O 

TCE-0.0 

PCA-O.  0 

DO  1  I* 1 , AZDIS 

1  70LT  (2,1)  »AZ  (I) 

IP  (ELSCAN.  E3.  1)  30  TO  12 
DO  2  I  *  1 , E  LDI 3 

2  70LT  (1,1)  *EL  (I) 

DIS-ELDIS 

GRAN*GRANE 

BW-EL33 

J*1 

3  1-1 
INC*  1 
PEDGE-O. 0 
SEDGE-0. 0 


BEGIN  HA  IN  LCCP 


FIND  SIZE  Of  WINDOW  AH  D  DETEHNIN!  NECSSSARr  ?  A3  AS  ETERS 


I 


TO  DETERMINE  TH2  PBOPER  SIZE  WINDOW,  WE  FIRST  NEED  TO  CALCULATE 
THE  OPTIMUM  WINDOW  SIZ2  WITH  A  "LOCK2D  OH"  ANTENN  k  BEAN,  THE 
C 2 HTTP  3EAM  POIHTIHG  LOCATION  32THG  NORB1L  OH  THE  TARGET. 

ir  AH  UNLOCKED  ANTENNA  BEAN  IS  USED,  THE  OPTIHU3  WINDOW  SIZE  IS 
TH2  SANE  IS  IN  THE  LOCKED  OH  CASE.  WE  THEN  DETEPSINE  THE  MAXIMUM 
HUN322  OP  BEAN  LOCATIONS  WITH  THE  WIDER  GRANULARITY  TO  NAINTAIN 
THE  OPTINUN  WINDOW  SIZE. 


HIDPT-  (DIN*  I)  /2 

A 4TH- FLO A?  (NIDPT)  /4  .  0 

A4THBI»AINT(A4TH) 

RE  SAIN-  (A4TH-A4TH3I)  *4.0 
PULSES-INT  (A4THHI) 

LIMIT«DIS-6*PULS2S 
COUNT-4  *  POLS  ES- 1 
Q-0 

IP  (LOCK.  EQ.  1)  GO  TO  » 

GHANOP-2. *8W*0.0174S33/( FLOAT (DIN)- I.) 

W IDT HO— PULSE  S-GRANOP 

PACTOR-3IDT  HO/GR  AN-  AI HT  (  WIDT  HO/GRAN) 

IP  (FACTOR. GT. 0.7  5) PULSES-INT (WIDTHO/GRAH) ♦ I 
IP  (PACTOR.LE.  0.7S)  PULSZS-INT  (W  IDTHO./GBA  H) 

IP  (PULSES.  EQ.O)  PULSES-  1 

LIBIT«OIN-«*PULSES 

COUNT-4-PUL3ES-1 


BEGIN  INNE3  LOC? 


CALCULATE  70LTAGE  WINDOWS 


4  IF  (Q.  EC.  LIMIT.  OR.  (ISC.  EQ.-1.  AND.  I.LE.  INT  (FEDGE)  >  )  GO  TO  II 
COUNT-COUNT* 1 

Q-0*  1 

DO  5  JJ-1,4 

5  7LT  (JJ)  -0. 0D0 
DO  7  J  J  J  —  1 .4 

(C-I*  (PULSES*(JJJ-1)  ) 

IP  (JJJ.EQ.  1  .AND. INC. 2Q.  1)  TEND- 70 L?  (J,K) 
l- 1*  PU LS IS  *JJJ-1 

IP  (JJJ.EQ. 4. AND. INC.  EQ.  -  1)  7END-70LT  (J,  K) 

IP  (JJJ.EQ. 2)  HIT-PLOAT  (L)  *0.5 
DO  6  II-K.L 

6  7LT(JJJ)  »7LT(JJJ)  »70LT(J,II) 

7  7LT  (JJJ) -7LT  (JJJ) /FLOAT  (PULSES) 

IP  (INC.EQ.-I)  GO  TO  9 

A-  1 
B-2 

C-3 
D-  4 


ooo  oonno 
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GO  TO  9 
9  A*» 

3-3 
C-2 
D*  I 


C 

C  FI  HO  DBS  17  ATIV  ES 

C 

. . ••*«•****•«••»«•»••••*•••*••••*••**•••*••••••••••••*••.•.••••..• 

c 

C  THE  riSST  DERIVATIVE  OS  SLOPE  3ETNEEN  THO  ADJACENT  WINDON 

C  VOLTAGES  IS  THE  INNER  SINUS  THE  OOTEH  70LTA3E,  DIVIDED  ST  A  UNIT 

C  ANGLE-  THE  SECOND  DERIVATIVE  OR  RATS  OP  CHANGE  OF  SLOPE  IS  THE 

C  INNER  NINOS  THE  OOTEH  FIRST  DERIVATIVE.  SINCE  THE  SAXIHOS  SLOPE 

C  OCCOSS  9HEN  THE  SECOND  DERIVATIVE  IS  ZERO,  SB  LOOK  FOR  THE  CHANGE 

C  IN  SIGN  Or  THE  SECOND  DERIVATIVE. 

C 

C  IF  THE  CHANGE  IN  SIGN  IS  NOT  FOOND,  OB  IF  THE  SAIIHOH  SLOPE  IS 

C  LESS  THAN  THE  NOISE  CRITERION,  THE  SINDONS  ARE  ADVANCED  OS  THE 

C  SCAN  REJECTED  ACCOBDINGLY. 

C 

. . . . . . . 

C 

C  REJECT  SCAN  IF  SEAN  VOLTAGE  3  IS  BELON  THE  VOLTAGE  BEFERENCE. 

C 

C,. 

9  IF  (VLT  (8) -LT.VREP)  GO  TO  10 
FDAB-SNGL  (VLT  (8)  -  VLT  (A) ) 

FDBC-SNGL  (VLT  (C)  -  VLT  (3)) 

SDB*  FDBC-FOAB 

IF  (SOB. LI. 0.0)  GO  TO  10 

FDCD-SNGL  (VLT  (D)  -VLT(C) ) 

SDCFDCD-FDBC 
IF  (SDC.  GT.  0)  GO  TO  10 

IF (NOISE- EQ. 1 .AND.FCBC.LT. N0SLI3)G0  TO  10 
IF  (POLSZS.  NS-  1.AND.PDCD.  LT»0.0)GO  TO  10 
IF  (INC. EQ.-1 ) GO  TO  IN 
FEDCE-HIT 


SETUP  FOR  SECOND  EDGE 


.■•■.•••..•.■•.•■••.••.a**..*.****...**.**'*-. 

1*01 . 1-4 • POLSZS* 1 
0*0 
INC--1 

C00NT*C00NT*4*  POL  SES- 1 
GO  TO  » 

ADVANCE  BINDONS 

10  I*I*INC 

SOB-SOa.TSHO 
N*N*  1. 

TBEF-SUH/N 
GO  TO  • 


non  no  on non  non non 
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IP  THE  PR0G3ES  REACHED  THIS  POINT,  T8ES  NO  SET  OP 
RETURNS  PISSED  IKE  CRITERION  PCR  THIS  SCAN. 


1  1  IP  <J.  S3.  2)  GO  TO  13 
COUNTE-COUNT 
5100!E»0 


SETUP  POR  EZI.1UT3 


12  di»-ezdi3 

GREN“GRAN A 
3i“E  Z3 i 
J«2 

GO  TO  3 

13  80DEE-0 
COUNTE “COUNT 
GO  TO  IS 


flSD  TARGE?  CENTER 


1*  SEOCE“HIT 

IE  (J.  EQ.  2)  GO  TO  IS 

TCE“ ( (5EDGE*PEDGE) /2. -PLCET  (BIDPT) ) •GRENZ*C3LOCE 
RODES-2 
COUNTS “COUNT 
GO  TO  12 

15  PCE“  (  (SEDGE *FEDGE)  /2. -PLCET  (3IDPT)  )  •GBESE»C3I.OCE 
RODEE-2 

COO RTE “COUNT 

16  LEBEL-5 

CELL  EBSOa  (TCE.PCE) 

IET0IN 

END 


END 


DATE 

FILMED 


DTIC 


